[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