Dredge MuMIn：在 GLMM 上使用挖泥机时，空模型是否包含随机效应？答案

【问题标题】：Dredge MuMIn: when using a dredge on a GLMM, does the null model include random effects?Dredge MuMIn：在 GLMM 上使用挖泥机时，空模型是否包含随机效应？
【发布时间】：2020-06-08 04:55:01
【问题描述】：

我有一个使用 GLMM 的全局模型，该模型具有 5 个具有交互作用的固定效应，以及两个随机效应。

x ~ a*b + a*c + a*d + a*e + (1|f) + (1|g)

我正在使用信息论方法，并在此全局模型上使用了 MuMIn 中的 dredge() 函数。

空模型的排名非常低 (95)，我认为这意味着固定效应对系统很重要，因为有 94 种可能的模型组合优于随机模型。

但是，一位同事询问随机效应是否包含在空模型中，否则，当然许多模型都比空模型更好，因为它们考虑了随机效应的变化，但空模型不会。

但是，我假设模型仍然会有这些随机效应，只是没有固定效应。这是正确的吗？

【问题讨论】：

标签： r data-modeling information-theory mumin

【解决方案1】：

是的。

这是一个实验说明，dredge() 仅丢弃固定效果，将空模型保留为（截距 + 所有随机效果）

图书馆（lme4）图书馆（木敏） m1 = lmer(Reaction~Days+(Days|Subject),sleepstudy,na.action=na.fail,REML=FALSE) dd = 疏浚(m1) 固定术语是“（拦截）” ## 全局模型调用：lmer(formula = Reaction ~ Days + (Days | Subject), data = sleepstudy, ## REML = FALSE，na.action = na.fail) ## --- ##型号选择表 ## (Intrc) 天数 df logLik AICc 增量权重 ## 2 251.4 10.47 6 -875.970 1764.4 0.0 1 ## 1 257.8 5 -887.738 1785.8 21.4 0 ## AICc(x) 排名的模型 ## 随机项（所有模型）： ##‘天数|主题”

这里的第一点是打印输出的末尾明确显示了所有术语包含的随机效应。

手动拟合空模型并将AICc与上表进行比较：

m0 <- update(m1, . ~ . - Days)
formula(m0)
## Reaction ~ (Days | Subject)
AICc(m0)
##    [1] 1785.821

【讨论】：