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

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

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

我使用的是信息论方法，并在此全局模型上使用了

MuMIn

中的

druge（）

函数

空模型排名很低（95），我认为这意味着固定效应对系统很重要，因为有94种可能的模型组合比随机模型更好

然而，一位同事问，零模型中是否包括随机效应，否则，当然许多模型比零模型好，因为它们考虑了随机效应的变化，但零模型不会

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

是的

这里有一个实验性的例子，

druge（）

只删除固定效果，将空模型保留为（截距+所有随机效果）

图书馆（lme4） 图书馆（MuMIn） m1=lmer（反应天数+（天数|受试者），睡眠研究，na.行动=na.失败，REML=FALSE） dd=挖泥船（m1） 固定术语为“（截取）” ##全球模型调用：lmer（公式=反应~天数+（天数|受试者），数据=睡眠研究， ##REML=FALSE，na.action=na.fail） ## --- ##选型表 ##（国际）日对数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

m0是

这里有一个实验性的例子，
druge（）
只删除固定效果，将空模型保留为（截距+所有随机效果）

图书馆（lme4）
图书馆（MuMIn）
m1=lmer（反应天数+（天数|受试者），睡眠研究，na.行动=na.失败，REML=FALSE）
dd=挖泥船（m1）
固定术语为“（截取）”
##全球模型调用：lmer（公式=反应~天数+（天数|受试者），数据=睡眠研究，
##REML=FALSE，na.action=na.fail）
## ---
##选型表
##（国际）日对数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

m0