我有一个在 SAS 中开发的混合效应模型,我试图在 R 中重新创建。因变量 Y 是 3 个固定效应 A、B 和 C 以及随机效应 D 的函数。D 有 3 个级别AxB 的每个组合,使其嵌套在 AxB 2-way 交互中。在 SAS 中,代码如下:
proc mixed data = S;
class A B C D;
model Y = A B A*B C C*A C*B A*B*C / outp = S2 influence residual;
repeated C/ subject = D(A*B);
random D(A*B);
run;
我正在尝试在 R 中使用 复制相同的lmer
内容,但无法在 AxB 交互中嵌套随机因子 D。有人可以帮我这样做吗?