排除在lme4

Question

我正在混合模型，使用lme4，其中包括三個因素（S，M，R）與2,3和4級中的一個因子水平相互作用的計算。但我失蹤了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)