[R-sig-ME] A three-level GLMM with binomial link in R
bbo|ker @end|ng |rom gm@||@com
Fri Mar 5 22:01:39 CET 2021
What kind of estimation issues? As far as just data set size is
concerned, I would consider this a "medium-to-large" problem for glmer,
but not something I would expect to cause problems.
* MASS::glmmPQL is built on top of lme(), so might work well for you
* glmmTMB provides a more-or-less drop-in replacement for glmer you
* GLMMadaptive::mixed_model is another possibility
* if you really need speed, Doug Bates's MixedModels.jl in Julia
usually beats everything else listed here
But ... we could probably give much better advice if you say
something about the specific kinds of "estimation issues" you're having.
On 3/5/21 3:27 PM, Hedyeh Ahmadi wrote:
> Hi All,
> I was wondering what would be a powerful package in R to run GLMM with logit link that can handle a data set with N=22945 and 3 nested random intercepts. So far, I have tried glmer() from lme4 and it's giving me a lot of estimation issues. Any other package I should try?
> I am asking for another package as I am having the same issue with lmer() for similar LMM with continuous outcome, while lme() from nlme package runs the models with no problem.
> Thank you in advance.
> Hedyeh Ahmadi, Ph.D.
> Keck School of Medicine
> Department of Preventive Medicine
> University of Southern California
> Postdoctoral Scholar
> Institute for Interdisciplinary Salivary Bioscience Research (IISBR)
> University of California, Irvine
> [[alternative HTML version deleted]]
> R-sig-mixed-models using r-project.org mailing list
More information about the R-sig-mixed-models