R语言survival包log-rank检验中t0时刻死亡事件的处理方式

核心问题确认 #

• R对t0时刻死亡事件的默认处理：是的，R的survival包在处理生存数据时，会自动对t0时刻的事件（进入时间=终止时间且发生事件）添加极小的偏移量，这类观测不会被排除，会被正常纳入分析计算。这是包内生存分析逻辑的默认设定，支持"瞬间事件"的统计处理。
• 修改默认行为的方法：若要强制排除这类t0时刻的事件，可手动筛选数据后再分析：

  # 假设数据框为df，进入时间start、终止时间time、事件变量status
  df_filtered <- df[df$time > df$start, ]
  # 基于筛选后的数据执行log-rank检验
  survdiff(Surv(start, time, status) ~ group, data = df_filtered)
  

R与Stata结果差异的原因确认 #

Stata提示的341 observations end on or before enter()，说明这些观测在Stata中被默认排除；而R默认保留了这类观测并做偏移处理，直接导致两者的分析样本量不一致，这是卡方统计量出现差异的核心原因。此外需核对Stata提示的1条缺失事件时间观测，确认R中是否做了完全一致的缺失值处理（R默认自动删除含缺失值的观测，需确保该环节逻辑匹配）。

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