使用lmer

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')) < / p>