r-如何获得glmer()模型的null偏差

Question

有没有一种方法可以获得用glmer((拟合的广义线性混合模型的零偏差和df？summary((输出中没有包括这一点，就像它与glm((对象一样，这是有原因的吗？

Answer 1

您可以通过仅使用截距项重新拟合模型来计算零偏差，例如

gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
data = cbpp, family = binomial)
gm0 <- update(gm1, . ~ 1 + (1|herd))
deviance(gm1)  ## 73.47
deviance(gm0)  ## 92.42 (null deviance)

• 我不知道你说的"；空df"；GLMM；有效样本量的"分母自由度"测量非常适用于平衡方差分析，也很难为GLMM定义线性混合模型[通过包含/排除、Satterthwaite、Kenward Roger等]
• 我可以想到lme4不会自动为您进行此计算的几个原因：
  • 这可能是一次昂贵的重新安装(即使对于GLM，它也需要重新安装模型，请参阅此处glm中的代码(
  • 对于GLMM来说，什么是合适的空模型进行比较并不明显。你是否同时去除了随机和固定效应，并将模型简化为GLM？你是保留所有的随机效果，还是只保留截距级别的随机效果或其他一些混合效果，这取决于问题的上下文？让用户自己做这件事迫使他们做出这个选择

(也就是说，我不认为省略零偏差是一个明确的选择。(

如果您确实选择放弃所有随机效应(即，与上例中的deviance(glm(cbind(incidence, size - incidence) ~ period, data =cbpp, family = binomial))进行比较，您应该能够与glmer拟合进行有意义的比较，但也有一些微妙之处：您可能想阅读?deviance.merMod中GLMM的偏差和对数似然性的章节。