### a fit Shared
data(readmission)
fit.shared <- frailtyPenal(Surv(time,event)~dukes+cluster(id)+
strata(sex),n.knots=10,kappa1=10000,kappa2=10000,data=readmission,
Frailty=TRUE)
### calling survival
survival(20,fit.shared)
Run the code above in your browser using DataLab