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