[R-sig-ME] Failure to converge with glmer on gamma data

moses selebatso selebatsom at yahoo.co.uk
Sun Mar 6 22:43:54 CET 2016


Hello
I m trying to determine the effects of season and habitat on forage protein content. With my replica nested within location. 
model<-glmer(Protein~Habitat*Season + (1|Location/Replica), family=Gamma)
I am getting a warning when running my model. I have tried read on previous post and most of the time I get lost because I fairly new to R. I understand there is a maximum "maxlgradl" of 0.001 that should be accepted in this kind of warnings. However, mine is a little higher than that. Can someone advise how I can proceed here and get an acceptable output.
I also tried using "....family =Gamma(log)"  instead of just "....family = Gamma" and it doesn't give me the warning. Can this be an acceptable option? 
Model and output below. I can provide data sample if needed.
> model<-glmer(Protein~Habitat*Season + (1|Location/Replica), family=Gamma)
Warning message:
In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
  Model failed to converge with max|grad| = 0.00123611 (tol = 0.001, component 1)
> model2<-glmer(Protein~Habitat+Season + (1|Location/Replica), family=Gamma)
Warning message:
In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
  Model failed to converge with max|grad| = 0.00254076 (tol = 0.001, component 1)


 Thank you

Moses SELEBATSO Home:    (+267) 318 5219 (H)  Mobile:  (+267) 716 39370  or  (+267) 738 39370
	[[alternative HTML version deleted]]



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