Stata中的st命令的R命令

Question

Stata中的stphtest的等效命令是什麼？

如果沒有一對一的命令，您如何測試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）。