https://mingredbestcheapprice.blogspot.com/2024/02/fitting-negative-binomial-in-python.html