https://lacucinadiadine.blogspot.com/2023/12/fitting-negative-binomial-in-python.html