[Ifeffit] High SO2

Scott Calvin SCalvin at slc.edu
Fri Apr 10 08:26:55 CDT 2009

Hi Abhijeet,

Well, I do know that many people do it that way, so it's not "wrong."  
But there are a number of problems with it:

* You're locking in errors, rather than using additional data to  
improve your fit. After all, the fit to the first shell has  
uncertainty, given by the error bars. But you're choosing the center  
point of the range and locking it in. By doing that, you're then  
distorting the fit to the next shell.

* You're ignoring "leakage" from the second shell when you fit the  
first shell. The second-shell paths have some width (look at them in  
Artemis to see this). Inevitably they have some tail that extends into  
the first shell. A fit to the first shell therefore has a systematic  
bias. If all you're doing is fitting one shell, then it may be small  
enough to live with it. But if you're adding a second shell anyway,  
why not let it inform the fit?

* It becomes difficult to keep track of degrees of freedom. By doing  
fits in layers, you of course aren't actually reducing the number of  
parameters you're fitting, but ifeffit doesn't know that. The reduced  
chi-square statistic is no longer reduced by the right number of  
degrees of freedom, so you lose the ability to use it to compare fits.  
And if you decide during the fitting of the second shell to change the  
k-range, degrees of freedom can get quite confusing.

Maybe someone who uses this strategy should speak up for it; I'd like  
to understand what the advantages are.

Just to be clear, what I'm suggesting instead is that you fit the  
first shell alone first, but when you add the second shell you  
continue to allow the first shell parameters to vary.

In any case, either way that you do it, you should be using the S02  
from the first shell for the subsequent shells. Otherwise you're not  
taking advantage of your understanding of the physical meaning of the  

--Scott Calvin
Sarah Lawrence College

On Apr 10, 2009, at 5:52 AM, abhijeet gaur wrote:

> Hello Scott Sir,
>                        I got that ifeffit basically works on the  
> principle of making the  theoretical data look similar to  
> experimental data. But what I know about fitting is that we have to  
> first fit the first shell and then keeping the first shell parameter  
> constant , we go for the next shell. Is it not right. Please tell me
>                                       Abhijeet Gaur
> _______________________________________________
> Ifeffit mailing list
> Ifeffit at millenia.cars.aps.anl.gov
> http://millenia.cars.aps.anl.gov/mailman/listinfo/ifeffit

More information about the Ifeffit mailing list