我有一个使用GLMM的全局模型,该模型具有5种具有交互作用的固定效果,以及两种随机效果。
x ~ a*b + a*c + a*d + a*e + (1|f) + (1|g)
我正在使用信息理论方法,并且在此全局模型上使用了dredge()
中的MuMIn
函数。
空模型的排名很低(95),我认为这意味着固定效应对系统很重要,因为有94种可能的模型组合要比随机模型好。
但是,一位同事问,是否将随机效应包括在null模型中,否则,许多模型当然比null更好,因为它们考虑了随机效应的变化,但是null模型没有考虑。
但是,我认为模型仍然会具有这些随机效应,只是没有固定效应。这是正确的吗?
是
[这是一个实验性的例证,dredge()
仅删除固定效果,而将空模型保留为(拦截+所有随机效果)
图书馆(LME4)图书馆(MuMIn)m1 = lmer(Reaction〜Days +(Days | Subject),sleepstudy,na.action = na.fail,REML = FALSE)dd =挖泥机(m1)固定术语是“(拦截)”##全局模型调用:lmer(公式=反应〜天+(天|主题),数据=睡眠研究,## 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