使用 lme4 进行随机效果 [关闭]

Question

我有 3 个随机变量，x,y z（全是随机效应）

x 嵌套在 y 中，但是 y 在 z 中交叉

我在 lme4 中使用了以下函数，但它不起作用。

<- lmer(A ~ 1 + (1 | x/y) + (1 | y*z) + (1|x/y*z), my data)

有人帮助我吗？非常感谢

Answer 1

恐怕这还很不清楚。更多的上下文将是有用的。我的猜测是你想要

A ~ 1 + (1|y)+ (1|z) + (1|y:z) + (1|y:x)

或等效

A ~ 1 + (1|y*z) + (1|y:x)

但几乎不可能确定。

• 前两个随机效应项给出了 y 间和 z 间的方差
• 第三项给出了 y 和 z 组合之间的方差 - 只有当您对每个 {y,z} 组合有多个观察值时，您才需要这个
• 最后一项给出了 x 嵌套在 y 中的效果。

表达式A ~ 1 + (1|y/x) + (1|z/y) 应该给您相同的结果，因为a/b 通常扩展为a + a:b（/ 的顺序很重要，但: 不重要），但不太清楚。

交叉随机效应通常由(1|y) + (1|z) 或(1|y*z)（扩展为(1|y) + (1|z) + (1|y:z)）表示。如果如上所述，每个 {y,z} 组合都有多个观察值。