En utilisant sim () avec lmer ()

Question

J'ai couru deux régressions logistiques multi-niveaux en utilisant les mêmes facteurs prédictifs, mais deux réponses différentes:

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 est discret avec des valeurs c(-2,-1,0,1,2), HDI est continue entre 0 et 1, Country est catégorique, et les réponses sont à la fois 1/0.

Pour tracer des bandes de confiance que je lance une simulation en utilisant la fonction sim () du paquet de bras:

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

Le premier calcule bien. Le second renvoie l'erreur suivante:

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

Je suis en train de faire cela avec 8 réponses différentes. Six d'entre eux a bien fonctionné et deux d'entre eux sont retournés cette erreur.

Quelqu'un sait comment corriger cela?

Answer 1

Pouvez-vous voir si le message d'erreur est correcte, à savoir, vérifier « à la main » si la matrice de covariance du modèle défaillant est en fait semi-définie positive? Je ne sais pas, mais je suppose qu'il ya des chèques que vous pourriez effectuer - peut-être il y a quelque chose ici: http://en.wikipedia.org/wiki/Positive-definite_matrix