[SciPy-user] looking for a negative binomial distribution

[Chris Fonnesbeck]

I notice in the scipy dev wiki that the negative binomial random number
generator is broken. In particular, it appears to round the first parameter,
which is incorrect -- any real number is valid, not just integers.

As a result, the method bombs out when passing it a parameter value of n<1, 
since it gets rounded to zero. This is a bit worrying for a common
distribution like the NB, and it makes me wonder about other random number
generators in numpy/scipy. 

Does anyone else know of a stable random number generating library for python?