title: Fitting Gamma distribution in Python
body: <p>I'm finding the parameters of Gamma distribution for a small sample. Later, I need to use the parameters to predict future data. However, the result shown incorrect answer.</p>
  
  <p>This is the result I got from Excel and it was the correct answer I'm looking for
  Alpha   0.458718895
  Beta    96.76626573</p>
  
  <pre><code>import scipy.stats as stats 
  data=[0.0621,0.046,0.0324,0.0279]
  fit_alpha, fit_loc, fit_beta=stats.gamma.fit(data,floc=0)
  print(fit_alpha, fit_loc, fit_beta)
  ll=[1,2,3,4,5,6,7,8,9,10]
  plop=stats.gamma.pdf(ll,fit_alpha, fit_loc, fit_beta)
  print(plop)
  </code></pre>
  
  <p>Expected result:
  6.29%   4.28%   3.40%   2.88%   2.53%   2.27%   2.06%   1.90%   1.76%   1.65%</p>