Использование sim () с lmer ()

Question

Я запустил две многоуровневые логистические регрессии, используя одни и те же предикторы, но по двум разным ответам:

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 () из пакета arm:

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