Survminer，生存軟件包：“ survdiff”對數秩檢驗與“常規”對數秩檢驗有什么區別？

Question

我在R中使用survminer生存軟件包。為了計算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" 。