从 clogit 运行疏浚模型子集

【问题标题】：Running a Dredge Model Subset from clogit从 clogit 运行疏浚模型子集
【发布时间】：2014-05-15 20:45:34
【问题描述】：

我正在尝试在我的完整模型 clogit（来自包 survival）上运行挖泥机，但每次尝试时 R 都会崩溃。我只是收到一条错误消息，说“R 遇到了一个致命错误。会话已终止。”

library(survival)
FullModel <- clogit(OBSERVED ~ Canopy + distgr_sca + Near_Hwy3 + strata(PID1), Compiled)
library(MuMIn)
dredge(FullModel)

这是我的代码。关于为什么会崩溃的任何想法，或者我可以在我的模型上运行完整挖泥机的其他方法？

谢谢。

【问题讨论】：

标签： r survival-analysis model-comparison

【解决方案1】：

当模型公式中未指定 strata 或 strata() 是唯一的模型项时，

clogit 似乎会使 R 崩溃。您需要告诉dredge 在所有模型中保留strata(PID1)，并将子集限制为至少有 3 个变量：

dredge(fm0, eval = T, fixed = "strata(id)", m.min = 3)

可重现的例子：

library(survival)
library(MuMIn)

## from example(clogit)
resp <- levels(logan$occupation)
n <- nrow(logan)
indx <- rep(1:n, length(resp))
logan2 <- data.frame(logan[indx,], id = indx, tocc = factor(rep(resp, each=n)))
logan2$case <- (logan2$occupation == logan2$tocc)
fm <- clogit(case ~ tocc + tocc:education + strata(id), logan2, na.action = "na.fail")
##

dredge(fm, fixed = "strata(id)", m.min = 3)

【讨论】：