在lme4中指定隨機效果

Question

我想在您指定我在lme4中正在研究的模型的隨機效應方面尋求幫助。 我有現場調查的數據。 該研究的目的是將機翼的大小（響應變量）與棲息地（探索性變量，具有2個級別的類別變量）相關聯。 我們通過在一個簡單的棲息地對50種3種物種的個體進行采樣，並在一個區域的復雜棲息地對3種物種進行了配對設計（再次是n = 50）。 我們在20個地區復制了該內容。 我對物種的特定影響和社區影響都感興趣。 對於運行每個物種的分析，此模型正確嗎？

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

對於整個社區的影響，此模型是否正確？

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

示意圖（來自保留的CrossValidated問題 ：

短數據集

Answer 1

我不會將“ Habitat同時用作固定效果和組級別效果，因為這樣您在組級別和固定效果之間就具有“完美的”相關性。 相反，我只會將Habitat用作固定效果：

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