r - 在 lme4 中排除一个因素水平的交互作用计算

Question

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

有没有办法告诉lme4忽略估计这个级别的参数？

而不是它给我的信息：Error in mer_finalize(ans) : Downdated X'X is not positive definite, 21.

Answer 1

手动构建交互变量，并删除级别。如果您给出一个可重现的示例，那就太好了，但是：

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