exponential_mixed_model function

Fit an exponential mixed-effects regression model