[R-sig-ME] lme4/glmer convergence warnings

Corey Sparks corey.sparks at utsa.edu
Fri Apr 11 03:01:58 CEST 2014


Hello all,
I’ve been seeing the aforementioned convergence errors for weeks in a course i’m teaching using lme4, so I decided to follow Ben’s advice on reporting the :

I’m fitting a binomial GLMM for small area estimation model building here:
glmer(I(bmi>30)~povz+vachousz+baccz+blackz+hispanicz+factor(region)+(1|state)+(1|cofips), family="binomial", data=merged, weights=cntywt/mean(cntywt))

n=~240000, n_cofips=217, n_state=46
I get the warning:
In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
  Model failed to converge with max|grad| = 0.235915 (tol = 0.001)

and the gradients:

relgrad <- with(fit.1 at optinfo$derivs,solve(Hessian,gradient))
> max(abs(relgrad))
[1] 0.0001631008

and when I use refit(), I get:
fit.1<-refit(fit.1)
> relgrad <- with(fit.1 at optinfo$derivs,solve(Hessian,gradient))
> max(abs(relgrad))
[1] 2.369877e-06

For another model on a LMM, I get:

fit.mix<-lmer(bmiz~agez+lths+coll+black+hispanic+other+(1|cofips), brfss_11)

In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
  Model failed to converge with max|grad| = 0.00550412 (tol = 0.002)

relgrad <- with(fit.mix at optinfo$derivs,solve(Hessian,gradient))
> max(abs(relgrad))
[1] 8.099289e-08

Hope this helps


Corey Sparks
Assistant Professor
Department of Demography
The University of Texas at San Antonio
501 West Cesar E Chavez Blvd
corey.sparks 'at' utsa.edu
coreysparks.weebly.com
210 458 3166



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