我正在尝试在下面拟合此模型:
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
有何相似。
警告在其消息中非常准确。
[用gaussian(link = "identity")
拟合混合效果模型时,等效于用正常随机效应拟合线性混合效果模型。
glmer
只是将调用更改为lmer(yield ~ Location + treatment + (1|block),data=df)
并给出警告。
警告已经存在很长时间了,我敢打赌,它不会在不久的将来被弃用,但是出于所有意图和目的,您应该使用lmer(...)
而不是glmer(..., family = gaussian(link = 'identity'))