title: Fitting negative binomial in python
body: <p>In scipy there is no support for fitting a negative binomial distribution using data
  (maybe due to the fact that the negative binomial in scipy is only discrete).</p>
  
  <p>For a normal distribution I would just do:</p>
  
  <pre><code>from scipy.stats import norm
  param = norm.fit(samp)
  </code></pre>
  
  <p>Is there something similar 'ready to use' function in any other library?</p>