scores_nbinom function

Calculating scores for the negative binomial distribution