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 using Penalized likelihood
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