将因子嵌套在R中的交互项中

Question

我在SAS中开发了一个混合效果模型，我试图在R中重新创建。因变量Y是3个固定效果的函数，A，B和C，以及随机效果D.有3个级别对于AxB的每个组合的D，使其嵌套在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;

我正在尝试使用lmer在R中复制相同内容，但无法在AxB交互中嵌套随机因子D.有人可以帮我这么做吗？