if (FALSE) {
#-- a fit Shared
data(readmission)
fit.shared <- frailtyPenal(Surv(time,event)~dukes+cluster(id)+
strata(sex),n.knots=10,kappa=c(10000,10000),data=readmission)
#-- calling survival
hazard(20,fit.shared)
}
Run the code above in your browser using DataLab