I used following command:
> frailtyPenal(Surv(time, status) ~var1 + cluster(family), Frailty=TRUE
> ,n.knots=8, kappa1=1500,
+ cross.validation=FALSE)
And got this error :
Error in Surv(time, status) : Time variable is not numeric
In addition: Warning message:
In is.na(time) : is.na() applied to non-(list or vector) of type 'closure'
-----------
It looks like you forgot data= in your model statement, which makes 'time'
refer to the function of that name.
Terry T.
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.