Regarding glm.nb and the theta parametre.

Hi, could you edit to make the code copy able, please, as reproducible example, called a reprex?

I know it's not all that long, but peeps are much more likely to take a look if they can just cut and paste.