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?
La solution
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