为什么这个生存图不从 100% 开始
Why doesn't this survfit plot start at 100%
当我用两个不同的审查事件绘制数据的生存图时,总体图 (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)
这似乎是 ggfortify 中的一个错误。作为临时修复,您可以通过执行以下操作将 t = 0 时的生存百分比设置为 100%:
p <- autoplot(kaplanMeier, facets = TRUE)
p$layers[[1]]$data[1, c(5, 7, 8)] <- 1
p
当我用两个不同的审查事件绘制数据的生存图时,总体图 (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)
这似乎是 ggfortify 中的一个错误。作为临时修复,您可以通过执行以下操作将 t = 0 时的生存百分比设置为 100%:
p <- autoplot(kaplanMeier, facets = TRUE)
p$layers[[1]]$data[1, c(5, 7, 8)] <- 1
p