title: trying to get reasonable values from scipy powerlaw fit
body: <p>I'm trying to fit some data from a simulation code I've been running in order to figure out a power law dependence. When I plot a linear fit, the data does not fit very well. </p>
  
  <p>Here's the python script I'm using to fit the data:</p>
  
  <pre><code>#!/usr/bin/env python
  from scipy import optimize
  import numpy
  
  xdata=[ 0.00010851,  0.00021701,  0.00043403,  0.00086806,  0.00173611, 0.00347222]
  ydata=[ 29.56241016,  29.82245508,  25.33930469,  19.97075977,  12.61276074, 7.12695312]
  
  fitfunc = lambda p, x: p[0] + p[1] * x ** (p[2])
  errfunc = lambda p, x, y: (y - fitfunc(p, x))
  
  out,success = optimize.leastsq(errfunc, [1,-1,-0.5],args=(xdata, ydata),maxfev=3000)
  
  print "%g + %g*x^%g"%(out[0],out[1],out[2])
  </code></pre>
  
  <p>the output I get is:
      -71205.3 + 71174.5*x^-9.79038e-05</p>
  
  <p>While on the plot the fit looks about as good as you'd expect from a leastsquares fit, the form of the output bothers me. I was hoping the constant would be close to where you'd expect the zero to be (around 30). And I was expecting to find a power dependence of a larger fraction than 10^-5. </p>
  
  <p>I've tried rescaling my data and playing with the parameters to optimize.leastsq with no luck.  Is what I'm trying to accomplish possible or does my data just not allow it? The calculation is expensive, so getting more data points is non-trivial.</p>
  
  <p>Thanks!</p>