Utilizzando sim () con lmer ()

Question

Ho eseguito due regressioni logistiche multilivello utilizzando gli stessi predittori, ma su due risposte diverse:

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 è discreta con valori c(-2,-1,0,1,2), HDI è continua tra 0 e 1, Country è categorica, e le risposte sono entrambi 1/0.

Per tracciare bande di confidenza ho eseguito una simulazione utilizzando la funzione di sim () dal pacchetto braccio:

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

La prima calcola bene. Il secondo restituisce il seguente errore:

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

Io in realtà sto facendo questo con 8 risposte diverse. Sei di loro ha funzionato bene e due di loro sono tornati questo errore.

Qualcuno sa come a rettificare questo?

Answer 1

Si può vedere se il messaggio di errore è corretto, vale a dire, controllare "a mano" se la matrice di covarianza dal modello non aver in realtà è positiva semi-definita? Non sono sicuro, ma immagino ci sono alcuni controlli si potrebbe effettuare - forse c'è qualcosa qui: http://en.wikipedia.org/wiki/Positive-definite_matrix