# NOT RUN {
data(udca)
# values found in table 8.3 of the book
fit1 <- coxph(Surv(futime, status) ~ trt + log(bili) + stage +
cluster(id), data=udca1)
fit2 <- coxph(Surv(futime, status) ~ trt + log(bili) + stage +
cluster(id) + strata(endpoint), data=udca2)
# }
Run the code above in your browser using DataLab