使用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}组合。