data(dataNested)
modNested<-frailtyPenal(Surv(t1,t2,event)~cluster(group)+subcluster(subgroup)+
cov1+cov2,data=dataNested,n.knots=8,kappa1=50000)
# It takes around 90 minutes to converge (depends on the processor) #
plot(modNested,conf.bands=FALSE)Run the code above in your browser using DataLab