title: Fitting For Discrete Data: Negative Binomial, Poisson, Geometric Distribution
body: <p>In scipy there is no support for fitting discrete distributions using data. I know there are a lot of subject about this.</p>
  <p>For example if i have an array like below:</p>
  <blockquote>
  <p>x = [2,3,4,5,6,7,0,1,1,0,1,8,10,9,1,1,1,0,0]</p>
  </blockquote>
  <p>I couldn't apply for this array:</p>
  <pre><code>from scipy.stats import nbinom
  param = nbinom.fit(x)
  </code></pre>
  <p>But i would like to ask you up to date, is there any way to fit for these three discrete distributions and then choose the best fit for the discrete dataset?</p>