排除lme4中的一个因子级别交互计算

Question

我使用lme4制作混合模型，其中包括三个因子（S，M，R），2,3和4级。但是我错过了对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)