lmerTestを使用して次のモデルを実行しましたlme4。
model2 = lmer(log(RT)~Group*A*B*C+(1|item)+(1+A+B+C|subject),data=dt)
コマンドを入力するlmerTest と、次のエラーが表示されます。summary()
> summary(model1)
Error in `colnames<-`(`*tmp*`, value = c("Estimate", "Std. Error", "df", :
length of 'dimnames' [2] not equal to array extent
これはすでに他のユーザーの問題であり、1 人のユーザーが実行中の問題を回避できたことがわかりましたlsmeans()。lsmeans を試したところ、次のエラーが発生しました。
Error in asMethod(object) : not a positive definite matrix.
共分散行列を調べると、NA は見当たりませんでした。グループ因子のコントラストを逆にするだけで、このモデルを実行できることに注意してください。なぜそうなのか理解に苦しむ。
lmerTest ではなく lme4 を使用して同じモデルを実行すると、summary() のすべての出力を取得できますが、p 値は取得できません (予想どおり)。pvals.fnc は lme4 で廃止されました。私はまだ代替手段を見つけていません。さらに、lmerTest をうまく使用できた他のモデルと同じ方法で model2 の p 値を推定するとよいでしょう。
この時点で私が何をすべきか知っている人はいますか？どんな助けでも大歓迎です！