在lme4中指定随机效果

Question

我想请您帮助我指定我在lme4中正在研究的模型的随机效应。我有现场调查的数据。该研究的目的是将机翼的大小（响应变量）与栖息地（探索性变量，具有2个级别的类别变量）相关联。我们通过在一个简单的栖息地对50种3种物种的个体进行采样，并在一个区域的复杂栖息地对3种物种进行了配对设计（再次n = 50）。我们在20个地区复制了该内容。我对物种的特定影响和社区影响都感兴趣。 对于运行每个物种的分析，此模型正确吗？

modela<-lmer(Wingsize~Habitat+(1|Region))

对于社区范围的影响，此模型正确吗？

modelb<-lmer(Wingsize~Habitat+(1|Region/Habitat/Species))

示意图（来自on-hold CrossValidated question：

Short Dataset

Answer 1

我不会同时将Habitat用作固定效果和小组效果，因为这样，您在小组级别和固定效果之间就具有“完美”的关联。相反，我只会将Habitat用作固定效果：

modelb<-lmer(Wingsize~Habitat+(1|Region/Species))