lmer（）でsim（）を使用して

Question

同じ予測因子を使用して2つのマルチレベルロジスティック回帰を実行しましたが、2つの異なる応答で次を実行しました。

fruitMLM <- lmer(InsuffFruit ~ Income + HDI + Income:HDI + (1 + Income | Country),family=binomial(link="logit"))  
fuelMLM <- lmer(Pollution ~ Income + HDI + Income:HDI + (1 + Income | Country),family=binomial(link="logit"))

Income 値は離散です c(-2,-1,0,1,2), HDI 0から1の間で連続しています、 Country カテゴリーであり、応答は両方とも1/0です。

信頼バンドをプロットするには、ARMパッケージからSIM（）関数を使用してシミュレーションを実行します。

sim(fruitMLM,100)  
sim(fuelMLM,100)

最初のものは正常に計算されます。 2つ目は次のエラーを返します。

Error in mvnorm(n.sims, bhat[j,], V.beta) :  
  'Sigma' is not positive definite

私は実際に8つの異なる応答でこれを行っています。そのうちの6人はうまく機能し、そのうち2人はこのエラーを返しました。

誰かがこれを修正する方法を知っていますか？

Answer 1

友人が私に答えを与えた。

プロジェクトの参照では、Copy localをfalseに設定する必要があります。

trueの場合、Visual Studioツールは成果物からのDLLとローカルに見つかったDLLの間でめちゃくちゃになります。

問題を解決しただけでなく、セーフコントロールエントリも正しく設定されています。