使用lmer的Glmer警告消息

Question

我正在努力适应以下模式：

cmod_lme4C_L <- glmer(yield ~ Location + treatment + (1|block),data=df,
                      family=gaussian(link = "identity"))

运行此代码后，我收到警告消息：

Warning message:
In glmer(yield ~ Location + treatment +  :
  calling glmer() with family=gaussian (identity link) as a shortcut to lmer() is deprecated; please call lmer() directly

有人能帮我理解这条信息吗？

看起来它建议使用lmer，但我不确定lmer与本例中的glmer有何相似之处。

Answer 1

这一警告在它的信息中非常准确。

当用gaussian(link = "identity")拟合混合效应模型时，它等价于拟合具有正态随机效应的线性混合效应模型。

glmer只需更改对lmer(yield ~ Location + treatment + (1|block),data=df)的调用并发出警告。

这个警告已经存在很长时间了，我敢打赌在不久的将来它不会被反对，但是出于所有的意图和目的，你应该使用lmer(...)而不是glmer(..., family = gaussian(link = 'identity'))。