library(KMsurv)
data(bmt)
fit.km.bmt=survfit(Surv(t2,d3)~group,data=bmt)
plot(fit.km.bmt,xlab="Days Post Transplant",
ylab="Estimated Disease-free Survival",
col=c(1,2,4))
legend(locator(1),c("ALL","AML-low","AML-high"),
lty=1,col=c(1,2,4),bty="n")
Как добавить срединную информацию о времени и степени опасности на графике в км, как показано ниже, спасибо.
