R中Cox模型的似然比检验

Question

是否有人知道似然比测试,如lmtest包中的lrtest,适用于使用coxph生成的cox比例风险模型？lrtest似乎不适用于coxph模型.

谢谢

Answer 1

存在一个anova.coxphpkg:survival,它可以比较模型对象.

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测试.