在Stata中用于stphtest的R命令

Question

St中的stphtest的R中的等效命令是什么？

如果没有一对一命令，您如何在Cox比例 - 危害回归模型中测试比例风险的假设？

Answer 1

请参阅?survival:::cox.zph：

  

描述

     

测试Cox回归模型拟合（coxph）的比例风险假设。

Answer 2

只是为了扩大格伦的答案： 我知道R对于Stata用户来说有多么不同。 首先运行coxph模型：

model.coxph <- coxph(Surv(t1, t2, event) ~ var1 + var2, data = data)

> summary(model.coxph0)
                      coef  exp(coef)   se(coef)      z Pr(>|z|)   
var1             1.665e-01  1.181e+00  1.605e-01  1.038  0.29948   
var2            -1.358e-03  9.986e-01  6.852e-04 -1.982  0.04746 * 

在模型上运行cox.zph：

(viol.cox <- cox.zph(model.coxph0))

            rho chisq     p
var1     0.0486 1.095 0.295
var2    -0.0250 0.258 0.611
GLOBAL       NA 1.462 0.691

每个p值应该是> 0.05这样就不会违反这些假设。 您还可以在cox.zph上使用简单的plot绘制缩放的schoenfeld残差：

> plot(viol.cox)