我已经在网上搜索了几个小时，但没有一个在线帖子是我要找的。我的问题在 SAS Proc 混合过程中很容易实现，但我不确定如何在 lme 和/或 lmer 包中实现。假设，我有一个模型， ，其中是固定的，但和是随机的。我的 R 代码是$y = \mu + \alpha + \beta +\alpha\beta + e$$\alpha$$\beta$$\alpha\beta$

 f1 = lme(y ~ factor(a), data = mydata,
     random = list(factor(b) = ~ 1, factor(a):factor(b) = ~ 1))

错误：意外=在：

 f1 = lme(y ~ factor(a), data = mydata,
          random = list(factor(a) =  

有人可以告诉我如何在 lme 中指定这些随机效果吗？提前谢谢了