Thanks for your advices, Terry and Konrad, using the linear fit as initial condition for the pawerlow fit works pretty well for my data. (I already had the two calculations but performed them vice versa ... :-) Anyway, I had the impression that the leastSquaresFit in Scientific Python is an implementation of the Levenberg Marquardt algorithm as it is presented in the Numerical Recipes. Accoring to reviews, this algorithm is not famous for its stability (e.g. http://www.stanford.edu/class/cme302/wnnr/nr.html). Better implementations are out there (e.g. http://www.ics.forth.gr/~lourakis/levmar/). Are there any plans to improve the SciPy algorithm? Would it be a welcome contribution to SciPy to work this part out?
- harold - -- http://mail.python.org/mailman/listinfo/python-list
