我目前在 glmer 中拟合模型(尽管同样的问题适用于 lmer)。
但是我有两个(我认为是)嵌套的固定效果。我正在测试人们与他们的性别和体重相比花多长时间进食。显然,某人的体重取决于他们的性别(性别之间的体重分布不同),但我在我的模型中关心这一点。我想测试最重的女性是否比最轻的女性吃得更多,最重的男性是否比最轻的男性吃得更多,但我想要一个模型来看看是否存在人口范围效应。
这就是我想通过将两个固定效果也作为嵌套随机效果来做到这一点的方式:
model <- glmer(TimeEating ~ Sex + Weight + (1|Sex:Weight),
data, family = poisson(link = "log"))