如何使用R包对Cox回归模型执行vuong测试

Question

据说vuong检验可用于比较两个未嵌套的回归模型。但是我想知道如何为cox回归模型做到这一点。 谢谢

Answer 1

从技术角度来看，

好吧，vuong函数接受来自类glm，negbin或zeroinfl的回归模型对象。由于coxph是不同类的对象，因此您可能不能将其提供给vuong函数。您可以将coxph馈送到plrtest，这是针对非嵌套coxph模型的部分似然比测试：

require("survival")
pbc  <- subset(pbc, !is.na(trt))
mod1 <- coxph(Surv(time, status==2) ~ age, data=pbc, x=T)
mod2 <- coxph(Surv(time, status==2) ~ age + albumin + bili + edema + protime, data=pbc,  x=T)
mod3 <- coxph(Surv(time, status==2) ~ age + log(albumin) + log(bili) + edema + log(protime), data=pbc, x=T)
plrtest(mod3, mod2, nested=F) # non-nested models
plrtest(mod3, mod1, nested=T) # nested models

现在，从理论角度来看，vuong是一种基于似然比的测试。从统计的角度来看，您可能要确保可以将它们用于Coxph。有些posts声称可以，有些声称shouldn´t。