SAS で開発した混合効果モデルを R で再作成しようとしています。従属変数 Y は、3 つの固定効果 A、B、C と変量効果 D の関数でした。D には 3 つのレベルがありました。 AxB の各組み合わせで、AxB 双方向相互作用内にネストされます。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をネストできません。誰かが私にこれを手伝ってもらえますか？