negative binomial distribution: follow R argument convention
- Not yet tested