如何提取R语言Kaplan-Meier生存图Y轴随访末期累积风险值

你拟合的survmod1对象已经包含了所有累积风险的计算结果，无需从图中手动提取，可直接通过如下代码获取两组随访结束时的累积风险百分比：

提取代码 #

# 方法1：直接查看随访最大时间点的模型结果摘要
# 输出结果中的1 - survival值乘以100就是你要的累积风险百分比
summary(survmod1, times = max(frmgham_recoded$timedth, na.rm = TRUE))

# 方法2：直接整理成结构化的数据框输出
library(dplyr)

# 将模型结果转为数据框，匹配分组信息
surv_res <- data.frame(
  time = survmod1$time,
  survival = survmod1$surv,
  group = rep(c("非吸烟者", "吸烟者"), survmod1$strata)
)

# 筛选每个组随访终点的累积风险值
surv_res %>%
  group_by(group) %>%
  filter(time == max(time)) %>%
  mutate(累积风险百分比 = round((1 - survival)*100, 2)) %>%
  select(group, time, 累积风险百分比)

结果说明 #

• 你绘图时使用fun="event"参数，本质就是绘制1 - 生存率的曲线，和上述代码计算的累积风险完全一致。
• 输出的结果就是两组随访期结束时的累积风险百分比，和图中Y轴的数值完全对应。

内容的提问来源于stack exchange，提问作者Natasha H