我正在尝试生成一条KM曲线,并希望审查事件像往常一样显示在图表中(如十字架或等效曲目)。
我正在使用上面的R版本,并更新了包“生存”,“KMsurv”和“rms”。
我尝试在某些网站上显示的代码没有成功。
来自http://www.visualcinnamon.com/2013/07/plotting-survival-analysis-results-in-r.html
的示例library(survival)
data(colon)
fit = survfit(Surv(time,status)~rx, data=colon)
plot(fit, xlab="Time", ylab="Survival Probability", main="Kaplan-Meier plot")
我确实生成了KM曲线,但是没有以图形方式表示删失事件。
我也试过(使用rms包):
fit <- npsurv(Surv(time,status)~rx, data=colon)
survplot(fit)
没有任何效果。我错过了什么?
PS:基础R中的mark.time工作,但我特别感兴趣的是累积概率图(AKA“倒置KM”),我不认为这可以在基地R中轻松完成。这对rms来说非常容易使用参数:'fun = function(x){1 - x}'