data(kidney)
##-- Shared frailty model --##
modSha <- frailtyPenal(Surv(time,status)~age+sex+cluster(id),
Frailty=TRUE,n.knots=8,kappa1=10000,data=kidney,hazard="Splines")
##-- Cox proportional hazard model --##
modCox <- frailtyPenal(Surv(time,status)~age+sex,n.knots=8,
kappa1=10000,data=kidney,Frailty=FALSE,hazard="Splines")
#-- confidence interval at 95
summary(modSha)
summary(modCox)
#-- confidence interval at 99
summary(modSha,level=0.99)
summary(modCox,level=0.99)Run the code above in your browser using DataLab