I am a student in M1 in experimental economics. I am struggling with a glmer model because I don’t know if the R function clusters the standard errors and if including a random intercept already had a similar effect (of clustering the standard errors). More precisely, I introduced a random intercept for controlling the intra-individual variability in my participants decisions (they take 18 decisions), but I am afraid to finish with biased results. In fact, I see that my residuals are not normally distributed (see below) and I am hoping it’s just because of this …

Btw, I’ve found this previous Overflow topic : /questions/72971844/cluster-robust-standard-errors-for-lmer-and-glmer-in-stargazer-lme4-package but it wasn’t really informative (and it didn’t work).
And this topic that let me think I must cluster my SE: /questions/22179631/individual-random-effects-model-with-standard-errors-clustered-on-a-different-va
but the function is different and I know it allows to cluster the SE quite easily.

enter image description here

The model:

summary(m1<-glmer(sharing~iv_realf*condf+(1|ResponseID)+(1|replicates)+age+gender, family=binomial, data=subset(d3, ac==2)))

Thus, I am wondering if I still have to cluster the SE (maybe the glmer function automatically does it) and if the random intercept accounts for the intra-variability from my participants answers.

Thanks in advance and sorry if it’s not clear (I don’t know a lot about econometrics and coding)