clesleycode / regression Goto Github PK

In your regression notebook's Python code for curve_fit() you have fixed parameter values coded as the p0 starting point, and if p0 is not given then scipy uses default initial parameter values of all 1.0 to start the "descent" into error space. This can be suboptimal, and in more complex equations often results in the algorithm finding a local minimum. For this reason, the authors of scipy added a genetic algorithm for initial parameter estimation. The module is named scipy.optimize.differential_evolution.

I have used scipy's Differential Evolution genetic algorithm to determine initial parameters for fitting a double Lorentzian peak equation to Raman spectroscopy data for carbon nanotubes and found that the results were excellent. The GitHub project, with a test spectroscopy data file, is:

https://github.com/zunzun/RamanSpectroscopyFit

My background is in nuclear engineering and industrial radiation physics, and I love Python, so if you have any questions, please let me know.

James Phillips