https://pypi.pythonprogramminghelp.com/how-to-create-negative-binomial-regression/