r - 为什么这个生存图不从 100% 开始

Question

当我用两个不同的审查事件绘制数据的生存图时，总图 (s0) 不是从时间 = 0、pstate = 100% 开始，而是在第一个审查事件发生时跳到 100%。

在这里你可以看到一个例子，跳跃发生在时间 1，这是第一个 cencoring 事件。

library(survival)
library(ggfortify)
library(tidyverse)
set.seed(1337)

dummy_data = tibble(time = sample.int(100, 100, replace = TRUE),
                    event = sample.int(3, 100, replace = TRUE))%>%
  mutate(event = factor(event))

kaplanMeier <- survfit(Surv(time, event) ~ 1, data=dummy_data)
autoplot(kaplanMeier, facets = TRUE)

Answer 1

这似乎是 ggfortify 中的一个错误。作为临时修复，您可以通过执行以下操作将 t = 0 时的生存百分比设置为 100%:

p <- autoplot(kaplanMeier, facets = TRUE)

p$layers[[1]]$data[1, c(5, 7, 8)] <- 1
p