R survminer包surv_median函数95%CI估计方法及与SAS匹配问询

survminer包的surv_median()函数未独立实现置信区间计算逻辑，其输出的中位生存时间95%CI完全提取自底层依赖的survival包survfit()函数的计算结果。

survfit()默认采用Brookmeyer-Crowley法计算中位生存时间的95%置信区间，该方法基于Kaplan-Meier生存曲线的翻转置信带推导得到区间估计值，恰好也是SAS PROC LIFETEST过程计算中位生存时间95%CI的默认方法。

若需保证双方计算结果完全一致，请注意以下配置要求：

• 两边需使用完全一致的生存时间定义、删失标记规则、数据截断规则
• SAS端使用PROC LIFETEST时无需修改METHOD=、CL=等参数，保留默认配置即可
• 你可在R端直接调用survival::survfit(Surv(时间, 删失) ~ 分组变量, data = 数据集)打印结果，其输出的中位生存时间及95%CI与surv_median()的输出完全匹配

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