mailr20264 - in /branches/relax_disp: docs/latex/dispersion.tex user_functions/relax_disp.py


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by edward on July 09, 2013 - 14:35:
Author: bugman
Date: Tue Jul  9 14:35:06 2013
New Revision: 20264

URL: http://svn.gna.org/viewcvs/relax?rev=20264&view=rev
Log:
Fix of the two-point dispersion error formula in the docs.

This includes the relax_disp.select_model user function docstring and the 
relax manual.


Modified:
    branches/relax_disp/docs/latex/dispersion.tex
    branches/relax_disp/user_functions/relax_disp.py

Modified: branches/relax_disp/docs/latex/dispersion.tex
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/docs/latex/dispersion.tex?rev=20264&r1=20263&r2=20264&view=diff
==============================================================================
--- branches/relax_disp/docs/latex/dispersion.tex (original)
+++ branches/relax_disp/docs/latex/dispersion.tex Tue Jul  9 14:35:06 2013
@@ -176,7 +176,7 @@
 
 Errors are calculated using the formula
 \begin{equation} \label{eq: dispersion error}
-    \sigma_{\textrm{R}_2}^2 = \frac{\left( 
\frac{\sigma_{I_1}}{I_1(\omega_1)} \right)^2  +  \left( 
\frac{\sigma_{I_0}}{I_0} \right)^2 }{T_\textrm{relax}} .
+    \sigma_{\textrm{R}_2} = \frac{1}{T_\textrm{relax}} \sqrt{ \left( 
\frac{\sigma_{I_1}}{I_1(\omega_1)} \right)^2  +  \left( 
\frac{\sigma_{I_0}}{I_0} \right)^2 } .
 \end{equation}
 
 In a number of publications, the error formula from \citet{IshimaTorchia05} 
has been used.  This is the collapse of Equation~\ref{eq: dispersion error} 
by setting $\sigma_{I_0}$ to zero:

Modified: branches/relax_disp/user_functions/relax_disp.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/relax_disp/user_functions/relax_disp.py?rev=20264&r1=20263&r2=20264&view=diff
==============================================================================
--- branches/relax_disp/user_functions/relax_disp.py (original)
+++ branches/relax_disp/user_functions/relax_disp.py Tue Jul  9 14:35:06 2013
@@ -494,12 +494,7 @@
     R1rho(nu1) = ------- * ln | ------- | ,
                  relax_T      \   I0    /\
 """)
-uf.desc[-1].add_paragraph("where I1 is the peak intensity in a spectrum for 
a given spin-lock field strength nu1.  Errors are calculated using the 
formula:")
-uf.desc[-1].add_verbatim("""\
-                    / (sigma_I1 / I1(nu1))^2 + (sigma_I0 / I0)^2 \ 
-    sigma_R2 = sqrt | ------------------------------------------ | ,
-                    \                  relax_T                   /\
-""")
+uf.desc[-1].add_paragraph("where I1 is the peak intensity in a spectrum for 
a given spin-lock field strength nu1.")
 uf.desc[-1].add_paragraph("For the variable relaxation time period type 
experiments, the R2eff/R1rho values are determined by fitting to the simple 
two parameter exponential as in a R1 or R2 analysis.  Both R2eff/R1rho and 
the initial peak intensity I0 are optimised using the minimise user function 
for each exponential curve separately.  Monte Carlo simulations are used to 
obtain the parameter errors.")
 # The no exchange model.
 uf.desc.append(Desc_container("The model for no chemical exchange 
relaxation"))




Related Messages


Powered by MHonArc, Updated Tue Jul 09 14:40:02 2013