# NOT RUN {
# }
# NOT RUN {
#-- here is generated cluster (30 clusters)
readmissionNested <- transform(readmission,group=id%%30+1)
# Baseline hazard function approximated with splines with calendar-timescale
model.spli.AG <- 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)
summary(model.spli.AG)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab