Survminer，生存包："survdiff" log-rank 和 "regular" log-rank 测试有什么区别？

Question

我正在使用 R 中的 survminer 生存包。为了计算 p 值以比较生存曲线，我正在使用 surf_pvalue 函数。

使用此函数，您可以通过两种不同的方式确定对数秩 p 值。一个称为 "survdiff"，另一个称为“1”或 "LR"。有关详细信息，请参阅 this link。第二种说是一个"regular log-rank test, sensitive to detect late differences"，但是两者有什么区别呢。有人知道吗？

谢谢。

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 的 source code，我们会看到

[...]
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"。

Survminer，生存包："survdiff" log-rank 和 "regular" log-rank 测试有什么区别？

Survminer, survival package: What is the difference between the "survdiff" log-rank and "regular" log-rank tests?

r

statistics

survival-analysis

survival

statistical-test