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

Question

据说vuong test可以用来比较两个没有嵌套的回归模型。但我想知道如何为 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。