add some stats function for negative binomial distribution