r - lme4で1つの因子レベルの交互作用計算を除外します

Question

lme42、3、4レベルの3つの要素（S、M、R）を含む、を使用して混合モデルを作成しています。しかし、S3：M2：R1の観測値が欠落しているため、3方向の交互作用を確認できません。

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)