survdiff log-rank 和 coxph log-rank 之间的区别

机器算法验证 r 生存 cox模型 对数秩检验
2022-04-08 08:08:51

我正在使用 R 中的生存包来分析临床数据。我正在分析两组不同的患者,当我计算 survdiff 以比较曲线时,我得到 p = 0.135,但是当我使用 coxph 和不同的协变量调整模型时,比如说临床癌症分期,我得到了一个整体 logrank 5 个协变量的得分为 0.0005793。我的问题是,我可以使用这个后期的 logrank p 值来说明用更多协变量调整模型曲线之间的差异是显着的吗?

这是数据

survdiff(formula = my.surv ~ final_table$G)

n=56, 14 observations deleted due to missingness.

                 N Observed Expected (O-E)^2/E (O-E)^2/V
final_table$G=1  4        2     1.43    0.2294     0.247
final_table$G=2 52       24    24.57    0.0133     0.247

 Chisq= 0.2  on 1 degrees of freedom, **p= 0.619** 

这是 coxph 结果

coxph(formula = Surv(final_table$Time_surv, final_table$Survival) ~ final_table$G + final_table$ST)

  n= 56, number of events= 26 
  (14 observations deleted due to missingness)

                   coef exp(coef)  se(coef)     z Pr(>|z|)
final_table$G2    2.094e-01 1.233e+00 7.532e-01 0.278    0.781
final_table$STII  1.883e+01 1.501e+08 5.739e+03 0.003    0.997
final_table$STIII 1.998e+01 4.773e+08 5.739e+03 0.003    0.997
final_table$STIV  2.089e+01 1.186e+09 5.739e+03 0.004    0.997

              exp(coef) exp(-coef) lower .95 upper .95
final_table$G2    1.233e+00  8.111e-01    0.2817     5.396
final_table$STII  1.501e+08  6.662e-09    0.0000       Inf
final_table$STIII 4.773e+08  2.095e-09    0.0000       Inf
final_table$STIV  1.186e+09  8.430e-10    0.0000       Inf

Concordance= 0.74  (se = 0.057 )
Rsquare= 0.37   (max possible= 0.957 )
Likelihood ratio test= 25.86  on 4 df,   p=3.381e-05
Wald test            = 4.02  on 4 df,   p=0.4033
Score (logrank) test = 19.67  on 4 df,   **p=0.0005793**

谢谢

感谢评论我做了这个分析,survdiff with roup and stage

survdiff(formula = Surv(final_table$Time_surv, final_table$Survival) ~ 
final_table$G + final_table$ST)

n=56, 14 observations deleted due to missingness.

                                 N Observed Expected (O-E)^2/E (O-E)^2/V
final_table$G=1, final_table$ST=III  3        2    1.149     0.630     0.668
final_table$G=1, final_table$ST=IV   1        0    0.279     0.279     0.285
final_table$G=2, final_table$ST=I   15        0    8.715     8.715    13.547
final_table$G=2, final_table$ST=II   2        1    1.816     0.367     0.402
final_table$G=2, final_table$ST=III 30       19   13.067     2.693     5.540
final_table$G=2, final_table$ST=IV   5        4    0.973     9.413     9.935

Chisq= 23.2  on 5 degrees of freedom, p= 0.000313

所以最终的值是完全显着的,但现在我得到了 6 条曲线,或多或少这就是我想要的,组和阶段如何影响生存。你怎么看?

1个回答

没有实际输出很难判断,但通常“整体对数秩分数”将检验所有系数均为 0 的零假设。因此,显着结果可能是由于您的一个或多个协变量与生存相关,而您的2组仍然相同(或者它们可能不同)。

最好用您的组变量(和协变量)拟合模型,并在没有组变量的情况下拟合另一个模型(但仍然使用相同的协变量)并比较 2 个拟合。