Utilizzando sim () con lmer ()
Domanda
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?
Soluzione
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