0

我已经使用相同的预测变量运行了两个多级逻辑回归，但有两个不同的响应：

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)

第一个计算得很好。第二个返回以下错误：

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

我实际上是用 8 种不同的反应来做这件事的。其中六个工作正常，其中两个返回此错误。

有谁知道如何纠正这个？

4

1 回答 1

1

您能否查看错误消息是否正确，即“手动”检查来自失败模型的协方差矩阵是否实际上是半正定的？我不确定，但我想你可以执行一些检查 - 也许这里有一些东西：http ://en.wikipedia.org/wiki/Positive-defined_matrix

于 2010-03-02T22:55:22.263 回答