使用与SIM卡11聚物（）（）

Question

我已经运行使用相同的预测器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。

要绘制置信带我运行使用SIM（）函数从臂封装的仿真：

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

第一个计算精细。第二个返回以下错误：

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

我其实有8个不同的响应这样做。其中六运行良好，其中两个返回此错误。

有谁知道如何纠正呢？

Answer 1

您可以看到，如果错误信息是正确的，即，勾选“手动”从故障模型的协方差矩阵是否实际上是半正定的？我不知道，但我想有一些检查，你可以执行 - 也许这里有什么东西： http://en.wikipedia.org/wiki/Positive-definite_matrix