r - 如何在 nlme 包中建模嵌套随机效应

Question

我在线性混合模型中有两个因素。因子 A 被视为固定效应，因子 B 被视为随机效应并嵌套在因子 A 中。谁能告诉我如何使用 nlme R 包做到这一点？

我知道 lme( response~ factorA, random=~1|factorA/factorB) 是一种建模方式。但是，此函数将因子 A 视为随机效应。

Answer 1

您的数据结构似乎使得对个人的观察嵌套在由 factorB 标识的组中。这些组进一步嵌套在由 factorA 标识的更大的分组中。您不希望此层次结构的最高级别拥有自己的随机截距项。相反，您只需使用作为固定效应包含的因素来模拟变化。美好的。然后剩下的就是允许截距在 factorB 上变化。这正是 DWin 所建议的。

lme(response ~ factorA, random=~1|factorB)

不过，这是否真的是您想要的并不完全清楚。从您写的问题中也不清楚您的数据的真实结构是什么。如果您更新问题，我将更新此答案以适合。