if (FALSE) {
data(drs)
# Clustered by subject
fit.drs <- fitfrail(Surv(time, status) ~ treated + cluster(subject_id),
drs, frailty="gamma")
fit.drs
# Variance estimates
vcov(fit.drs)
# Plot the estimated cumulative baseline hazard
plot(fit.drs, type="cumhaz")
}
Run the code above in your browser using DataLab