[R-sig-ME] Convergence warnings in lme4

Ben Bolker bbolker at gmail.com
Tue Oct 21 15:39:49 CEST 2014


On 14-10-21 09:21 AM, adam.hayward at ed.ac.uk wrote:
> Dear list users,
> 
> This is a general question, rather than a specific example of my data. I
> have recently upgraded my computer and to R ver 3.1.1 from 3.0.2. Models
> (glmer) which previously ran perfectly are now producing warning
> messages regarding failure to converge, for example:
> 
> Warning messages:
> 1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
>   Model failed to converge with max|grad| = 0.0238312 (tol = 0.001,
> component 13)
> 2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
>   Model failed to converge: degenerate  Hessian with 1 negative eigenvalues
> 
> and this is true with a range of optimizers. Running identical models on
> my old and new versions of R produces identical output (estimates, AIC
> etc) but the new version of lme4 produces the warnings.
> 
> My questions are thus (1) has lme4 recently been updated with
> significant changes in what constitutes acceptable convergence, and (2)
> can my models now be trusted?
> 
> Many thanks for any comments or help,
> 
> Adam

  Presumably you're using version 1.1-7 (the latest version on CRAN); it
would be helpful to know that.  You can see changes to lme4 at
http://cran.r-project.org/web/packages/lme4/news.html . The significant
changes you're referring to started in version 1.1-3.

  We are still working to figure out how to get the right level of
sensitivity and specificity on these tests; we haven't come up with
anything definitive yet.  I've got an example of various attempts at
troubleshooting (restarting the fit from the same point, trying
different optimizers, etc.) posted at
http://rpubs.com/bbolker/lme4trouble1 .

The bottom line is that your models are no more nor less trustworthy
than they used to be.  If you have tried them with a variety of
optimizers, and if the results make sense, it's *probably* the case that
they're OK, and just a little unstable.

  Ben Bolker



More information about the R-sig-mixed-models mailing list