if (FALSE) {
data(drs)
fit.drs <- fitfrail(Surv(time, status) ~ treated + cluster(subject_id),
drs, frailty="gamma")
# Plot the parameter and log-likelihood trace
plot(fit.drs, type="trace")
# This may take a while to run.
# Use parameter B to specify the number of repetitions in the weighted bootstrap
plot(fit.drs, type="cumhaz", CI=0.95)
}
Run the code above in your browser using DataLab