glm.nb_mverse function

Fit negative binomial regression models across the multiverse