logistic_mixed_model function

Fit a logistic mixed-effects regression model