r - Survminer，生存包：“survdiff”对数秩和“常规”对数秩测试有什么区别？

Question

我正在使用survminerR 中的生存包。要计算 p 值以比较生存曲线，我正在使用该surf_pvalue函数。

使用此函数，您可以通过两种不同的方式确定对数秩 p 值。一个称为“survdiff”，另一个称为“1”或“LR”。有关更多信息，请参阅此链接。第二个说它是“常规对数秩检验，对检测后期差异很敏感”，但两者之间有什么区别。有人知道吗？

谢谢你。

Answer 1

方法"survdiff", "log-rank", "LR","1"都给出相同的结果。

这是一个最小且可重现的示例

library(survival)
library(survminer)

fit <- surv_fit(Surv(time, status) ~ sex, data = colon)

methods <- c("survdiff", "log-rank", "LR", "1")
purrr::imap(setNames(methods, methods), ~surv_pvalue(fit, method = .x))
#$survdiff
#  variable      pval   method pval.txt
#1      sex 0.6107936 Log-rank p = 0.61
#
#$`log-rank`
#  variable      pval   method pval.txt
#1      sex 0.6107936 Log-rank p = 0.61
#
#$LR
#  variable      pval   method pval.txt
#1      sex 0.6107936 Log-rank p = 0.61
#
#$`1`
#  variable      pval   method pval.txt
#1      sex 0.6107936 Log-rank p = 0.61

---

如果我们看一下源代码，survminer我们会看到

[...]
allowed.methods <- c("survdiff", "log-rank", "LR", "1",
                       "n", "Gehan-Breslow", "GB",
                       "sqrtN", "Tarone-Ware", "TW",
                       "S1", "Peto-Peto", "PP",
                       "S2", "modified Peto-Peto", "mPP",
                       "FH_p=1_q=1", "Fleming-Harrington(p=1, q=1)", "FH")

method.names <- c(rep("survdiff", 4),
                rep(c("n", "sqrtN", "S1", "S2", "FH_p=1_q=1"), each = 3))
[...]

四种方法统一为method.names = "survdiff".