# 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 DataCamp Workspace