R语言survival包：骨髓移植数据2年生存率比较的截尾疑问

• 不需要手动截断数据集到2年
  你完全不用把生存时长上限设为730天（2年）再调整截尾状态。survival包的生存分析模型（比如Kaplan-Meier）本身就支持直接聚焦到指定时间点的生存率计算。对于那些随访超过2年的患者，在2年这个时间点他们的状态天然就是截尾——我们只关心到2年的生存情况，2年之后的随访信息不会参与2年生存率的计算，保留原始数据也不会给分析添乱。

• 原始数据的“长尾巴”不会影响亚组间的2年生存率比较
  当你的分析终点是2年生存率时，Kaplan-Meier曲线在2年之后的部分（也就是你看到的长尾巴）对2年时间点的生存率估计、以及亚组间的统计检验（比如Log-rank检验）没有任何干扰。模型计算时只会用到每个样本在2年以内的信息：

  • 2年内死亡的患者：计入事件，用来计算风险率
  • 存活超过2年的患者：在2年节点被标记为截尾，只贡献到2年的随访时长，2年之后的情况不纳入分析

• R语言中的正确处理示例
  直接用原始数据拟合模型，然后提取2年的生存率即可，代码示例如下：

  library(survival)
  
  # 假设数据集名为bone_marrow，time是生存天数，status是结局（1=死亡，0=截尾），subgroup是亚组变量
  km_fit <- survfit(Surv(time, status) ~ subgroup, data = bone_marrow)
  
  # 提取并查看730天（2年）的各亚组生存率
  summary(km_fit, time = 730)
  

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