fit_mle function

Fit basic Poisson and Negative Binomial models using glmmTMB