Stata中stphtest的R命令

Question

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

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

Answer 1

只是为了扩展格伦的答案： 我知道 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)

Answer 2

见?survival:::cox.zph:

说明

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