# NOT RUN {
library(survival)
data(Breast)
dim(Breast)
km <- survfit(Surv(time, status) ~ treat, data = Breast)
km
plot(km, col = c("black", "red"), lwd = 2, xlim = c(0, 5), xaxt = "n", yaxt = "n")
legend("bottomleft", legend = c("Control", "Experimental"), col = 1:2,
lty = 1, lwd = 2, cex = 1.5)
axis(1, cex.axis = 1.3)
axis(2, las = 2, cex.axis = 1.3)
mtext("Distant-recurrence free survival", side = 2, line = 3.2, cex = 1.5)
mtext("Time (in years)", side = 1, line = 2.5, cex = 1.5)
# }
Run the code above in your browser using DataLab