r - R 中 Cox 模型的似然比检验

Question

有谁知道似然比检验，如 lmtest 包中的 lrtest ，它适用于使用 coxph 生成的 cox 比例风险模型？ lrtest 似乎不适用于 coxph 模型。

谢谢

Answer 1

pkg:survival 中有一个 anova.coxph，它允许比较模型对象。

fit <- coxph(Surv(futime, fustat) ~ resid.ds *rx + ecog.ps, data = ovarian) 
fit2 <- coxph(Surv(futime, fustat) ~ resid.ds +rx + ecog.ps, data=ovarian)
anova(fit2,fit)

Analysis of Deviance Table
 Cox model: response is  Surv(futime, fustat)
 Model 1: ~ resid.ds + rx + ecog.ps
 Model 2: ~ resid.ds * rx + ecog.ps
   loglik  Chisq Df P(>|Chi|) 
1 -31.970                    
2 -30.946 2.0469  1    0.1525

这是一个LR测试。