I have fit a random intercept model in glmer and stan_glmer with the following model formula where my event is 0/1 and x1 and x2 are continuous covariates:
model <- stan_glmer(binary event ~ x1 + x2 +(1 | group), family="binomial"
The predicted probabilities (i.e. the fitted values from the model) are not normally distributed. Does this mean the model is wrong?