# NOT RUN {
# }
# NOT RUN {
data(kidney)
##-- Shared frailty model --##
modSha <- frailtyPenal(Surv(time,status)~age+sex+cluster(id),
n.knots=8,kappa=10000,data=kidney,hazard="Splines")
##-- Cox proportional hazard model --##
modCox <- frailtyPenal(Surv(time,status)~age+sex,
n.knots=8,kappa=10000,data=kidney,hazard="Splines")
#-- 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)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab