Fit negative binomial mixed-effect regression model.