[R-sig-ME] about the glmer function
lme4 maintainer
bbolker at gmail.com
Tue May 20 20:32:38 CEST 2014
[forwarding to r-sig-mixed-models]
This is an unfortunate limitation of the current version of lme4. You
could try lme4.0, which did allow AGHQ on vector-valued random effects,
or wait until we can find enough resources to (re)implement this
feature. (Or fit your model in MCMCglmm, or live with the Laplace
approximation, or ...)
https://github.com/lme4/lme4/issues/123
sincerely
Ben Bolker
On 14-05-20 12:29 AM, Anıl AKTAS SAMUR wrote:
> Dear LME4 Authors,
>
> My outcome variable is binomial, and I have 11 independent variables and
> a time variable. The time variable has different slopes, so I fixed it
> to |time-before| and |time-after|. I used the |lme4| package (the
> |glmer| function). I have a random intercept and two random slopes. I
> created my model like this:
>
> |m3.glmm <- glmer(y ~ timebefore + timeafter + x1 + x2 +...+ x11 + (1+timebefore+timeafter|id),
> data = data, family = binomial (link="logit"), nAGQ=3)|
>
> When I used this model, I had this error:
>
> |Error in updateGlmerDevfun(devfun, glmod$reTrms, nAGQ = nAGQ) :
>
> nAGQ > 1 is only available for models with a single, scalar random-effects term|
>
>
> Could you please explain how to fit (or code) this model?
>
> Bests,
>
>
> Anil
More information about the R-sig-mixed-models
mailing list