[SciPy-User] optimize.leastsq and improper input parameters

ms devicerandom@gmail....
Thu Jul 8 12:46:07 CDT 2010


Hi,

I am stuck with optimize.leastsq.
I am writing a quite complicated code that fits data to two variations 
of the same function, each of those can be with or without some 
parameters fixed. The concept is:
- write the two variations
- have a list of what parameters will be fixed and which not
- use reasonable starting points of the non-fixed params to fit, and 
keep the fixed params as fixed within the function
- then have a generalized routine that has as input the function to 
actually fit, among other things
- apply leastsq on that function

What comes out is similar, in structure, to this simplified example:
-----
import scipy as sp
import scipy.optimize as opt
import numpy as np


#Initial data
x = np.arange(0,10,1)
y = [i**2 for i in x]

#Three nice functions to fit
def xexponent1(param,i):
     exp = param[0]
     return i**exp

def xexponent2(param, i):
     exp_a=param[0]
     exp_b=param[1]
     return i**(exp_a) + i**(exp_b)

def line(param,i):
     A = param[0]
     b = param[1]
     return i*b + A


def f_to_minimize(pars,args):
     #Generalized function we use to minimize the fit
     #calculates function and squared residuals
     function,x,y = args[0],args[1],args[2]
     y_estimate = [function(pars,xi) for xi in x]

     #calculate squared residuals
     resid = [(i-j)**2 for i,j in zip(y,y_estimate)]
     return sum(resid)

def minimize(x,y,func,p0):
     #calls minimization
     args = [func,x,y]
     i = opt.leastsq(f_to_minimize, p0, args)
     print i

minimize(x,y,line,[1,2])
-----

Here you can play with p0 and the function to give to f_to_minimize as 
an argument.

What comes out is that if I give a p0 = [1] and I use a single-variable 
function, it works. As soon as I try a two-variable function (and thus I 
need two input parameters), I get:

massimo@boltzmann:~/work$ python test_norm_leastsq.py
Traceback (most recent call last):
   File "test_norm_leastsq.py", line 44, in <module>
     minimize(x,y,[1,2])
   File "test_norm_leastsq.py", line 41, in minimize
     i = opt.leastsq(f_to_minimize, p0, args)
   File "/usr/lib/python2.6/dist-packages/scipy/optimize/minpack.py", 
line 300, in leastsq
     raise errors[info][1], errors[info][0]
TypeError: Improper input parameters.

The funny thing is that it worked *before* I messed with the thing to 
simplify the function-choosing mechanism (before I had N different 
functions for each combination of fixed/nonfixed params, now I just have 
two and I fix stuff *inside* the function), and I can't see however how 
can this be different. Also, the example above leaves me perplexed -it 
seems leastsq simply doesn't want two-variable functions to be minimized 
in this case. Any hint?

Thanks a lot,
Massimo


More information about the SciPy-User mailing list