require(survival)
leukemia.surv <- survival::survfit(survival::Surv(time, status) ~ x, data = survival::aml)
survival:::plot.survfit(leukemia.surv, conf.int = TRUE, lty = 2:3, col = 1:2)
qkmplot(leukemia.surv, conf_int = TRUE)
Run the code above in your browser using DataLab