lmer 的 drop1 函数

【问题标题】：drop1 function for lmerlmer 的 drop1 函数
【发布时间】：2015-01-03 06:39:39
【问题描述】：

我构建了一个混合效应模型，其中包含三个固定效应和一个随机效应。

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

【解决方案1】：

如果您只使用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？这可能很棘手/需要一些黑客攻击......

【讨论】：