data(readmission)
# Shared frailty model
modSha <- frailtyPenal(Surv(time,event)~dukes+cluster(id),
n.knots=10,kappa1=10000,data=readmission,Frailty=TRUE,
hazard="Splines")
print(modSha)
# Cox proportional hazard model using Penalized likelihood
modCox <- frailtyPenal(Surv(time,event)~dukes,n.knots=10,
kappa1=10000,data=readmission,Frailty=FALSE,hazard="Splines")
print(modCox)
Run the code above in your browser using DataLab