[R-sig-ME] Bglmer convergence warnings (covariance structure)

Ben Bolker bbolker at gmail.com
Mon Jun 12 18:04:46 CEST 2017


 I think if you'd failed to specify the covariance matrix correctly that
you would simply have gotten an error.
  I'm doing some test runs on simulated data. I did get convergence
warnings, but the results look reasonable when compared across widely
different fitting platforms (lme4, blme, glmmTMB, brms):

http://bbolker.github.io/mixedmodels-misc/notes/bglmer_cmp.html

code (Rmd file etc.) is in the corresponding github repo

  cheers
   Ben Bolker



On 17-06-12 06:45 AM, D.J. Damen wrote:
> Dear all,
> 
> I am trying to fit a binomial model with three factors using the bglmer package in R. My model includes a complex 3x2x2 design; factor c.con.tr has three levels, factor c.type.tr two, and factor c.diff.tr has also two levels. When I try to fit a random-intercept only model, the model produces the following warning:
> 
> Warning message:
> In get("checkConv", lme4Namespace)(attr(opt, "derivs"), opt$par,  :
>   Model failed to converge with max|grad| = 2.02439 (tol = 0.001, component 1)
> 
> I figured that I did not specify my covariance structure correctly, but I am at a loss as to how I should change it. The structure of my model is as follows:
> 
> riobglmer <- bglmer(contrast~c.con.tr*c.type.tr*c.diff.tr+(1|id)+(1|item.new), data=mydata, family=binomial (link='logit'), fixef.prior= normal(cov = diag(9,12)))
> 
> To me, (cov = diag(9,12)), seems to be correct.. Am I missing something important? Suggestions and/or remarks are more than welcome.
> 
> Thank you very much in advance.
> 
> Best regards,
> 
> Debby Damen
> PhD Student
> 
> Department of Communication and Information Science
> 
> Tilburg University
> Warandelaan 2, room D410
> 5037 AB Tilburg
> 
> T. +31 13 466  8245
> M. d.j.damen at uvt.nl<mailto:d.j.damen at uvt.nl>
> 
> 
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>



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