[R-sig-ME] glmer and influence.me - complaining about nAGQ==0

Cátia Ferreira De Oliveira cm|o500 @end|ng |rom york@@c@uk
Sat Apr 24 20:29:02 CEST 2021


Dear Professor Bolker,

I am really sorry for bothering you but I have been quite stuck and have
posted my question both on reddit and cross validated but nothing has come
from it. Do you know if there is a good way of running dfbetas for glmer
with a gamma distribution?  The model runs without issues,

lmer_log11 <- glmer(logRT ~ Probability*Block*Session*testing + (1 +
Block * Probability * Session|Participant), data= Data, family =
Gamma(link = "log"),
control=glmerControl(optimizer="bobyqa",optCtrl=list(maxfun=1000000)),
nAGQ = 0)

yet when I run the influence function I get this error:

influentialcases <- influence(glmer_log11, "Participant")

Error in glmer(formula = RT ~ Probability * Block * Session * testing +  :
  should not specify both start$fixef and nAGQ==0

I have tried to remove the "nAGQ==0" but the model is now taking 10 days to
converge, a lot more than usual. Given this, I haven't been able to check
if the influential cases will run normally.

Do you have any idea why this is happening? I would very much
appreciate your help!

Thank you so much!

Best wishes,

	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list