Запуск подмножества модели Dredge из clogit

Я пытаюсь запустить драгу на моей полной модели 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
person user3495408    schedule 03.04.2014    source источник

Ответы (1)


arrow_upward
1
arrow_downward

clogit вызывает сбой R, если в формуле модели не указан слой или когда strata() является единственным термином модели. Вам нужно указать dredge сохранить strata(PID1) во всех моделях и ограничить подмножества минимум тремя переменными:

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)
person Kamil Bartoń    schedule 04.04.2014