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. ______________________________________________ R-help@r-project.org 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.