Question

我使用R中的 coxme（）函数运行混合效果模型。该模型分析了不同国家公司的产品成功事件。固定效应例如是GDP，人口，技术和文化变量。随机效应是不同的国家。

我知道使用 coxph（）可以使用 cox.zph（）命令测试比例风险。

我的问题：如何使用 coxme（）检查比例风险？

Answer 1

可以使用与标准 coxme 模型相同的 cox.zph() 函数检查随机效应 coxph() 模型中的固定效应是否存在比例风险 (PH)。根据手册，fit 的 cox.zph() 参数是“使用 coxph 或 coxme 函数拟合 Cox 回归模型的结果。”

随机效应“不检查比例风险，而是将它们视为模型中的固定偏移量。”

一个例子，借自this Cross-Validated question：

> library(survival)    
> library(coxme)    
> df <- stanford2
> df$cid <- round(df$id / 10) + 1 ## generates some clusters
> fit <- coxme(Surv(time, status) ~ age + t5 + (1 | cid),data=df)
> fit
Cox mixed-effects model fit by maximum likelihood
  Data: df
  events, n = 102, 157 (27 observations deleted due to missingness)
  Iterations= 2 12 
                    NULL Integrated    Fitted
Log-likelihood -451.0944  -446.8618 -446.8261

                  Chisq   df        p  AIC   BIC
Integrated loglik  8.47 3.00 0.037317 2.47 -5.41
 Penalized loglik  8.54 2.04 0.014582 4.46 -0.88

Model:  Surv(time, status) ~ age + t5 + (1 | cid) 
Fixed coefficients
          coef exp(coef)   se(coef)    z      p
age 0.02960206  1.030045 0.01135724 2.61 0.0091
t5  0.17056610  1.185976 0.18330590 0.93 0.3500

Random effects
 Group Variable  Std Dev      Variance    
 cid   Intercept 0.0199835996 0.0003993443
> cox.zph(fit)
       chisq   df    p
age    0.831 3.00 0.84
t5     2.062 2.04 0.36
GLOBAL 2.767 5.04 0.74

这是通过 survival_3.1-11 和 coxme_2.2-16 完成的。

R：coxme中的比例风险假设（）

1 个答案: