0

同じ予測子を使用して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"))

Incomec(-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つはこのエラーを返しました。

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

4

1 に答える 1

1

エラーメッセージが正しいかどうかを確認できますか。つまり、失敗したモデルの共分散行列が実際に正の半確定であるかどうかを「手動で」確認できますか？よくわかりませんが、実行できるチェックがいくつかあると思います。おそらくここに何かがあります：http: //en.wikipedia.org/wiki/Positive-definite_matrix

于 2010-03-02T22:55:22.263 に答える