функция drop1 для lmer

Я построил модель смешанных эффектов с тремя фиксированными эффектами и одним случайным эффектом.

mdl1 <- lmer(yld.res ~ veg + rep + rip + (1|state),
       REML=FALSE,data=data2)

Я хочу получить самую экономную модель из вышеописанной модели. Для этого я хочу отбросить одну независимую переменную за раз и посмотреть, улучшит ли это соответствие модели (посмотрев на значение AICc). Но когда я использую drop1, это дает мне следующую ошибку:

drop1(mdl1, test="F")

Error in match.arg(test) : 'arg' should be one of “none”, “Chisq”, “user”

Я не совсем уверен, как это сделать, и был бы очень признателен за любую помощь.


r lme4
person user53020    schedule 06.11.2014    source источник

Ответы (1)


arrow_upward
3
arrow_downward

Если вы просто используете drop1() со значением по умолчанию test="none", это даст вам значения AIC, соответствующие модели с каждым фиксированным эффектом, отбрасываемым по очереди.

Вот немного глупый пример (вероятно, нет смысла тестировать модель с квадратичным, но без линейного члена):

library('lme4')
fm1 <- lmer(Reaction ~ Days + I(Days^2) + (Days | Subject), sleepstudy)
drop1(fm1)
## Single term deletions
## 
## Model:
## Reaction ~ Days + I(Days^2) + (Days | Subject)
##           Df    AIC
## <none>       1764.3
## Days       1 1769.2
## I(Days^2)  1 1763.9

Насколько сильно вам нужно AICc, а не AIC? Это может быть сложно/требуется некоторый взлом...

person Ben Bolker    schedule 06.11.2014