Ho costruito ed eseguito un modello di regressione logistica a effetti misti nel pacchetto lme4 per r per stimare la probabilità di occupazione di pesci in luoghi diversi (celle/habitat). Il quadro di dati è costituito da 1.207.140 osservazioni di 68 singoli pesci. Per ogni individuo (al giorno per ~ 1 anno) descrive il numero di occorrenze in ciascuna posizione unica rispetto al numero totale di occorrenze in tutte le località.Informazioni sui messaggi di avviso per il modello misto in r lme4
Qui è il modello base:
m.base = glmer(cbind(N,t.move-N) ~ jdate + snSurface.Area + Restoration..P.A. +
Release.Location+ Sex + (1|Station) + (0 + jdate|ID), data=allfishdat, family=binomial)
where N=# unique positions, t.move=total positions, jdate=julian date, Station=locations, ID=fish ID
ottengo il seguente messaggio di avviso:
Warning messages:
1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 3349.26 (tol = 0.001)
2: In if (resHess$code != 0) { :
the condition has length > 1 and only the first element will be used
3: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model is nearly unidentifiable: very large eigenvalue
- Rescale variables?;Model is nearly unidentifiable: large eigenvalue ratio
- Rescale variables?
ho fatto qualche ricerca per cercare di capire che cosa questi messaggi significato e le loro implicazioni sul modello, ma non hanno ancora capito gli avvertimenti.
Potete dire quello che hai guardato già, e che cosa fare e non capire? hai letto https://github.com/lme4/lme4/blob/master/README.md? http://stackoverflow.com/questions/23814130/glmer-model-from-early-2013-warning-message-about-convergence-when-re-running-i/23839952#23839952? http://stackoverflow.com/questions/21344555/convergence-error-for-development-version-of-lme4? –
Urto. Sarebbe utile se potessi dire quello che hai e non ho capito - sono felice di aiutarti ma non mi sento di ripetere informazioni che sono già là fuori ... –
Urto Bump ..... –