Usando sim () con LMER ()

Question

He corrido dos regresiones logísticas multinivel utilizando los mismos predictores, sino en dos respuestas diferentes:

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 es discreta con valores c(-2,-1,0,1,2), HDI es continua entre 0 y 1, Country es categórica, y las respuestas son ambos 1/0.

Para trazar bandas de confianza corro una simulación utilizando la función SIM () desde el paquete de brazo:

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

El primero calcula bien. El segundo devuelve el error siguiente:

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

De hecho, me estoy haciendo esto con 8 respuestas diferentes. Seis de ellos funcionaba bien, y dos de ellos regresó este error.

¿Alguien sabe cómo rectificar esto?

Answer 1

¿Puede usted ver si el mensaje de error es correcto, es decir, comprobar "a mano" si la matriz de covarianza del modelo no es en realidad semi-definida positiva? No estoy seguro, pero me imagino que hay algunas comprobaciones se podría realizar - tal vez hay algo aquí: http://en.wikipedia.org/wiki/Positive-definite_matrix