使用具有

我已經使用相同的預測運行兩個多級邏輯迴歸11聚物（）SIM（），但在兩個不同的反應：使用具有

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。

要繪製置信帶我運行使用SIM（）函數從臂封裝的仿真：

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

第一個計算精細。第二個返回以下錯誤：

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

我實際上正在做這8個不同的響應。其中六人工作正常，其中兩人返回此錯誤。

有誰知道如何解決這個問題？

來源

2010-02-25 Jared

什麼版本的'arm'？我知道這是早期版本中的問題（請參閱第6點：http://www.stat.columbia.edu/~cook/movabletype/archives/2007/09/recent_changes.html）。 – Shane 2010-02-25 22:44:25

如果您可以提供一些示例數據來重現問題，則調試問題也更容易。 – Shane 2010-02-26 12:36:00

Shane，對不起，我沒有意識到我收到了評論。我正在構建模型的人不提供一些模型的模擬。我用過的手臂可能已經過時。我想我最近再次嘗試了模擬，並安裝了新的手臂，並且工作正常。但我很可能很快會再次訪問，並會從那裏看到。感謝您的幫助。 – Jared 2010-04-27 18:11:05

您能否看到錯誤消息是否正確，即手動檢查來自故障模型的協方差矩陣是否實際上是半正定的？我不確定，但我想有一些檢查可以執行 - 也許這裏有一些東西：http://en.wikipedia.org/wiki/Positive-definite_matrix

來源

2010-03-02 22:55:22

感謝您的想法。我很快就會重新討論這個問題，所以我會嘗試一下。對不起，我沒有意識到我收到了回覆。 – Jared 2010-04-27 18:06:26

回答

相關問題