Kaplan-Maier曲线审查事件未显示检查员和事件时间之间的关系

Question

我们注意到R（survival package）和SAS（或survminer）如何显示KM曲线审查事件之间存在差异。当事件与审查事件同时发生时，图表中不会显示+标记。 survminer::ggsurvplot函数正确显示审查事件，如SAS。

我在这里遗漏了什么吗？如何让plot.survfit显示所有审查事件＆＃39;？

library(survival)
library(survminer)

df <- data.frame(time=c(4,6,8,11,15,15,17,18,19), 
                 cens=c(0,0,1, 0, 1, 0, 1, 0, 1)) 

mod <- survfit(Surv(time, !cens) ~ 1, data=df) 

# plotting with the survival package does not show a mark at time = 15 
plot(mod, mark=1, mark.time = TRUE, conf.int = FALSE)

# plotting with survminer or SAS display the mark at time 15.
ggsurvplot(mod, data = df, conf.int = FALSE)

Answer 1

找到了解决方案。

plot(mod, mark=1, conf.int = FALSE, mark.time=mod$time[mod$n.censor > 0])