# dontrun
# Load survival
# library(survival)
# Subset data
# bladder1 <- bladder[bladder$enum < 5, ]
# Estimate coxph model
# M1 <- coxph(Surv(stop, event) ~ (rx + size + number) * strata(enum) + cluster(id), bladder1)
# Survfit
# M1Fit <- survfit(M1)
# Plot strata in a grid
# ggfitStrata(M1Fit, byStrata = TRUE)
Run the code above in your browser using DataLab