https://camera-optical-zoom.blogspot.com/2022/10/fitting-negative-binomial-in-python.html