[R-sig-ME] Number of Iterations lme4

Jakob Gährken jakob.gaehrken at uni-kassel.de
Thu May 7 15:02:09 CEST 2015


Hello,
I was asked to forward my question to this mailing list. I would be glad 
to get some hints regarding my questions below:

I am an quit inexperienced R user but I try hard to learn fast. My 
problem is that while I try to run lme4 with the model:
model1 = glmer(Survival ~ Kohorte+KJ+Saison+KAE+(1|ID), family = 
binomial(link = "logit"), 
data=phenos,control=glmerControl(optimizer="bobyqa",optCtrl=list(maxfun=1e6)))
I always get the warning :

1:In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
    Model failed to converge with max|grad| = 0.75662 (tol = 0.001, component 2)
2:In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
    Model failed to converge: degenerate  Hessian with 1 negative eigenvalues

I know that my data has the promblem that the binary response variable has a skewed distribution (1.933 "0" vs. 68.031 "1"). And that this could be the reason for the converging problem and that I might look for alternatives.(btw. do you consider some kind of minimum proportion for the "minor observation" in this case "0" (2.8%) when you fit a model e.g. for diseases with low incidence)
But beside of the reason for the converging problem I am wondering because I tried to change the number of iterations using  : control=glmerControl(optCtrl=list(maxfun=...)).The used time which leads to the error message and also the  ...with max|grad| = 0.75662...are  always the same - no matter what value for maxfun I choose. So I am not sure how to control the number of iterations that are really done by R and whether the amount changes with different defined maxfun values in my case.

Thank you for your time and looking forward hearing from you,

Jakob Gährken



	[[alternative HTML version deleted]]



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