# NOT RUN {
# }
# NOT RUN {
data(dataNested)
modClu <- frailtyPenal(Surv(t1,t2,event)~cluster(group)+
subcluster(subgroup)+cov1+cov2,data=dataNested,
n.knots=8,kappa=c(50000,50000),hazard="Splines")
print(modClu)
#-- here is generated cluster (30 clusters)
readmissionNested <- transform(readmission,group=id%%30+1)
modJointNested_Splines <- frailtyPenal(formula = Surv(t.start, t.stop, event)
~ subcluster(id) + cluster(group) + dukes +
terminal(death), formula.terminalEvent = ~dukes,
data = readmissionNested, recurrentAG = TRUE, n.knots = 8,
kappa = c(9.55e+9, 1.41e+12), initialize = TRUE)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab