comp从survMisc包中尝试。它扩展了survival包。它计算对数秩检验的统计量和 p 值，以及 Gehan-Breslow、Tarone-Ware、Peto-Peto 和 Fleming-Harrington 检验以及趋势检验（对于上述所有内容）。取自手册的示例如下：

data(larynx, package="KMsurv")
s4 <- survfit(Surv(time, delta) ~ stage, data=larynx)
comp(s4)
comp(s4)$tests$trendTests # outputs only the results for trend tests

如果将结果与

survdiff(Surv(time, delta) ~ stage, data=larynx)

对于“传统”对数秩检验（不是趋势检验），您会得到相同的结果。

你的问题不是很清楚，所以不确定这是否是你要找的。要测试比例风险假设，您可以对比例风险模型的 Schoenfeld 残差使用Grambsch-Therneau 检验。这基本上测试了（缩放）残差的斜率作为后续时间的函数。