[R-sig-ME] Need help on convergence issue when fitting zero-inflated Poisson with random coefficients using gamlss

Ben Bolker bbo|ker @end|ng |rom gm@||@com
Thu Jan 23 20:28:55 CET 2020

    Since the proximal problem looks like (guessing from the error
messages) a call to the nlminb optimizer within a call to lme is
running out of iterations before converging, I'd try to see if there
is a way to set the number of parameters higher.  Unfortunately, I
don't know how to do this off the top of my head (a brief review of
the gamlss documentation didn't get me anywhere, and I haven't had
time to dig through the code to see if/how this is possible) ...

   There are other packages that are capable of fitting mixed ZIP
models (e.g. https://journal.r-project.org/archive/2017/RJ-2017-066/index.html
 shows examples using glmmTMB, inla, MCMCglmm, gamlss, ...) in case
you can't solve the problem with gamlss (and you aren't stuck using
gamlss for other reasons)

   Ben Bolker

On Wed, Jan 22, 2020 at 2:01 PM Xia Li <odditylee using gmail.com> wrote:
> Hello,
> I have trouble getting convergence when fitting ZIP with (individual level)
> random coefficients, the function that I used:
> m1 <- gamlss(y ~ re(fixed = ~ treatment, random = ~ treatment|unit_id),
>              family = ZIP, data = dat)
> Specifically I wanted to include global or population level treatment
> effect, and individual random effect nested under treatment.
> I always got the following error messages:
> Error in lme.formula(fixed = fix.formula, data = Data, random = random, :
>    nlminb problem, convergence error code = 1
>    message = iteration limit reached without convergence (10)
> Is there anything that I can debug from? Changing the control parameters? I
> tried different algorithms like method = CG() but seems it did not help
> much.
> Looking for help. Thanks!
> --
> Best,
> Xia
>         [[alternative HTML version deleted]]
> _______________________________________________
> R-sig-mixed-models using 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