numpy/numpy — pull request #4754 by ChienliMa (2014-05-29 00:28:11+00:00)
Title: modified distribution.c

In documentation , p is the possibility of success.
so scale θ = p/(1 − p) rather than (1-p)/p
cite: http://en.wikipedia.org/wiki/Negative_binomial_distribution