glm.nb function

Fit a Negative Binomial Generalized Linear Model