Curve fitting a sum of Gaussian's to 6 peaks
I have successfully used a genetic algorithm to search error space and find initial parameter estimates for scipy's Levenberg-Marquardt nonlinear solver.  The most recent versions of scipy now include the Differential Evolution genetic algorithm which is what I had been using.  It will take a bit of experimenting on a previously fit data set to tune things like population size or genetic crossover strategy, but then you should have a way for the computer to find initial parameter estimates in this way.  It takes a while to run, so let it go overnight, but it can be automated.
I was wondering if anyone could show me a way of fitting multiple gaussian curves to a data set containing 6 peaks(data comes from a diffraction pattern from a copper gold alloy crystal). the way i have at the moment involves using multiple gaussian equations added together meaning i have to give multiple guesses of values when curve fitting.