https://camera-optical-zoom.blogspot.com/2023/12/python-negative-binomial-regression.html