Dear subscribers:
Please spread the word to those who may be interested in our next annual
short course at BNL on X-ray absorption fine structure.
The dates are: November 2-4. It will be virtual. The level will be
intermediate to advance (last year it was a beginner level). As in the
past, lectures, discussions and software demonstration sessions will be on
Wednesday and Thursday, and the full day data analysis practicum will be on
Friday.
More information about the course and the link to the application form will
appear soon on BNL web site, and on the web site of the Synchrotron
Catalysis Consortium (SCC).
Regards,
Anatoly
-------------------------------------------------------------------
Anatoly I. Frenkel
Professor
Department of Materials Science and Chemical Engineering
Stony Brook University
Stony Brook, NY 11794
Email: anatoly.frenkel(a)stonybrook.edu
<https://mail.bnl.gov/OWA/redir.aspx?C=F6wCGt4adziN6IvpdRQ-rWdgy3Oaoif_x-0yb…>
http://you.stonybrook.edu/frenkel
Co-Director and Spokesperson, Synchrotron Catalysis Consortium (SCC) at BNL
http://you.stonybrook.edu/scc2
Joint Appointment:
Senior Chemist
Chemistry Division, Brookhaven National Laboratory
Upton, NY 11973. Ph: 631-344-3013. Group: 631-344-3494
Email: frenkel(a)bnl.gov
https://www.bnl.gov/staff/frenkel
-------------------------------------------------------------------
Hello all,
I have been a frequent user of Athena for many years, mostly for interpreting P K-edge XANES spectra. Until last week I thought that the R factor in Athena was always defined as:
sum( [data_i - fit_i]^2 )
-------------------------------
sum( data_i^2 ]
This is also the definition given in the online manual, and it has been stated by me and by other colleagues in a number of papers dealing with P K-edge XANES. But well, this is not true when dealing with normalized XANES spectra! I realized this when I played around with a number of my old LC fits in Excel. While the chi-square value (or maybe more precisely, the sum of squared residuals) was reproduced perfectly, I always got "R factors" (i.e. with the above definition) between 2 and 3 times lower than what Athena gives. After that I consulted the Demeter programming documentation (https://bruceravel.github.io/demeter/pods/Demeter/LCF.pm.html) to find that, for normalized mu(E), "Demeter thus scales the R-factor to make it somewhat closer to 10^-2". However, the equation stated on this page actually reproduces the R factor even more poorly, and therefore I won't reiterate it here. After inspecting the Perl code, and trying out different alternatives in Excel, I now believe that the following equation provides a more accurate definition of the R factor (correct me if I'm wrong!):
sum( [data_i - fit_i]^2 )
-------------------------------
sum( [data_i - avg data]^2 )
where "avg data" is the arithmetic mean of the data in the LC fitting range. It would be great if others could confirm this. As far as I understand, this won't affect the interpretations that any of us have made over the years, it only affects the understanding of what the R factor actually is...
Kind regards, Jon Petter
Jon Petter Gustafsson, Professor in Soil Chemistry
Department of Soil and Environment
Swedish University of Agricultural Sciences (SLU)
Box 7014
750 07 Uppsala, Sweden
Phone: 018-671284; e-mail: jon-petter.gustafsson(a)slu.se<mailto:jon-petter.gustafsson@slu.se>
---
N?r du skickar e-post till SLU s? inneb?r detta att SLU behandlar dina personuppgifter. F?r att l?sa mer om hur detta g?r till, klicka h?r <https://www.slu.se/om-slu/kontakta-slu/personuppgifter/>
E-mailing SLU will result in SLU processing your personal data. For more information on how this is done, click here <https://www.slu.se/en/about-slu/contact-slu/personal-data/>
