I want to extract standard deviation of residual from `glmer()` function in R .
So I wrote :
lmer_obj = glmer(Y ~ X1 + X2 + (1|Subj), data = D, family = binomial)
sigma(lmer_obj)
I noticed that the last command `sigma(lmer_obj)` returns always "1" irrespective of data ...