The model consists 4 fixed effects and a random effect. I'm keep getting a convergence error
Error: nlminb problem, convergence error code = 1 message = false convergence (8)
How do I solve this?
You can try different optimisers. nlme takes control argument that can used to set different optimisers:
control=lmeControl(opt = "optim")
Hopefully glmmPQL can pass that argument to nlme.
Thank you very much for your answer.
I did try this but glmmPQL didn't pass this argument, I called a different package but that overrides the package MASS which I need for GLMM.
You can test if nlme:lme converge with your data with some different optimers. If so, open option is to create wrapper function which calls nlme:lme with a different optimiser and pass other arguments through. Then you need to modify glmmPQL (change mcall[[1L]]
- Badges
- Science topic
- Similar topics
- Computer Science and Engineering
- Coding