R 排除lme4中的单因素水平交互计算

我正在使用

lme4

制作混合模型，其中包括三个因子（S、M、R），分别为2、3和4级。但是我缺少对S3:M2:R1的观察，所以我无法检查3向交互

有没有办法告诉

lme4

忽略此级别的参数估计

---

相反，它给我的信息是：

mer_finalize（ans）中的错误：降级的X'X不是正定的，21。

手动构造交互变量，并降低级别。如果您能给出一个可复制的示例，那就太好了，但是：

mydata <- transform(mydata,SMRint=interaction(S,M,R,drop=TRUE))
lmer(response~SMRint+(1|grpvar),data=mydata)

---

mydata或
interaction（S，M，R，drop=T）
谢谢！工作得很有魅力：）