R-factors

From CCP4 wiki
Jump to navigation Jump to search

Definitions[edit | edit source]

Data quality indicators[edit | edit source]

In the following, all sums over hkl extend only over unique reflections with more than one observation!

  • Rsym and Rmerge - the formula for both is:


[math]\displaystyle{ R = \frac{\sum_{hkl} \sum_{j} \vert I_{hkl,j}-\langle I_{hkl}\rangle\vert}{\sum_{hkl} \sum_{j}I_{hkl,j}} }[/math]


where [math]\displaystyle{ \langle I_{hkl}\rangle }[/math] is the average of symmetry- (or Friedel-) related observations of a unique reflection. The formula is due to Arndt, U.W., Crowther, R.A. & Mallet, J.F.W. A computer-linked cathode ray tube microdensitometer for X-ray crystallography. J. Phys. E:Sci. Instr. 1, 510−516 (1968). Any unique reflection with n=2 or more observations enters the sums.

It can be shown that this formula results in higher R-factors when the redundancy is higher (Diederichs and Karplus [1]). In other words, low-redundancy datasets appear better than high-redundancy ones, which obviously violates the intention of having an indicator of data quality!

  • Redundancy-independant version of the above:


[math]\displaystyle{ R_{meas} = \frac{\sum_{hkl} \sqrt \frac{n}{n-1} \sum_{j=1}^{n} \vert I_{hkl,j}-\langle I_{hkl}\rangle\vert}{\sum_{hkl} \sum_{j}I_{hkl,j}} }[/math]


which unfortunately results in higher (but more realistic) numerical values than Rsym / Rmerge (Diederichs and Karplus [1] , Weiss and Hilgenfeld [2]).

measuring precision of averaged intensities/amplitudes[edit | edit source]

for intensities use (Weiss [3])


[math]\displaystyle{ R_{p.i.m.} = \frac{\sum_{hkl} \sqrt \frac{1}{n-1} \sum_{j=1}^{n} \vert I_{hkl,j}-\langle I_{hkl}\rangle\vert}{\sum_{hkl} \sum_{j}I_{hkl,j}} }[/math]

Rmrgd-I (defined in Diederichs and Karplus [1]) only differs by a factor (FIXME: what is the factor? 0.5 or 1.4142 or ?) since it likewise takes the improvement in precision from multiplicity into account. Rsplit , which is what the X-FEL community uses, is the same as Rmrgd-I but that community seems not to be aware of this.

Similarly, one should use Rmrgd-F as a quality indicator for amplitudes [1], which may be calculated as:


[math]\displaystyle{ R_{mrgd-F} = \frac{\sum_{hkl} \sqrt \frac{1}{n-1} \sum_{j=1}^{n} \vert F_{hkl,j}-\langle F_{hkl}\rangle\vert}{\sum_{hkl} \sum_{j}F_{hkl,j}} }[/math]


with [math]\displaystyle{ \langle F_{hkl}\rangle }[/math] defined analogously as [math]\displaystyle{ \langle I_{hkl}\rangle }[/math].

In the sums above, the summation omits those reflections with just one observation.

measuring radiation damage[edit | edit source]

We can plot (Diederichs [4])


[math]\displaystyle{ R_{d} = \frac{\sum_{hkl} \sum_{|i-j|=d} \vert I_{hkl,i} - I_{hkl,j}\vert}{\sum_{hkl} \sum_{|i-j|=d} (I_{hkl,i} + I_{hkl,j})/2} }[/math]


which gives us the average R-factor of two reflections measured d frames apart. As long as the plot is parallel to the x axis there is no radiation damage. As soon as the plot starts to rise, we see that there's a systematical error contribution due to radiation damage.

Strong wiggles at very high d are irrelevant as only few reflections contribute.

To my knowledge, the only program that implements this currently (December 2008) is XDSSTAT.

Comparing two sets of structure factor amplitudes or intensities[edit | edit source]

The following is symmetric, and suitable for comparing two data sets, or two model amplitudes:


[math]\displaystyle{ R_{scale}=\frac{\sum_{hkl}\vert F_{hkl,i}-F_{hkl,j}\vert}{0.5\sum_{hkl} F_{hkl,i}+F_{hkl,j}} }[/math]


for amplitudes, and analogously for intensities.

Model quality indicators[edit | edit source]

  • R and Rfree : the formula for both is


[math]\displaystyle{ R=\frac{\sum_{hkl}\vert F_{hkl}^{obs}-F_{hkl}^{calc}\vert}{\sum_{hkl} F_{hkl}^{obs}} }[/math]


where [math]\displaystyle{ F_{hkl}^{obs} }[/math] and [math]\displaystyle{ F_{hkl}^{calc} }[/math] have to be scaled w.r.t. each other. R and Rfree differ in the set of reflections they are calculated from: R is calculated for the working set, whereas Rfree is calculated for the test set.

Relation between R and Rfree as a function of resolution[edit | edit source]

       mode(R) = 0.091*ln(resolution) + 0.134
       mode(Rfree-R)   = 0.024*ln(resolution) + 0.020

References:

  • Tickle IJ, Laskowski RA and Moss DS. Rfree and the Rfree Ratio. I. Derivation of Expected Values of Cross-Validation Residuals Used in Macromolecular Least-Squares Refinement. Acta Cryst. (1998). D54, 547-557 [5]
  • Tickle IJ, Laskowski RA and Moss DS. Rfree and the Rfree ratio. II. Calculation of the expected values and variances of cross-validation statistics in macromolecular least-squares refinement. Acta Cryst. (2000). D56, 442-450 [6]
  • Urzhumtsev, Afonine & Adams (2009) Acta Cryst., D65, 1283-1291.

what kinds of problems exist with these indicators?[edit | edit source]

  • (Rsym / Rmerge ) should not be used to judge data quality, Rmeas should be used instead. The reason is that the former depend on multiplicity, whereas the latter doesn't.
  • R/Rfree and NCS: reflections in work and test set are not independent if chosen randomly. It is better to choose the test set reflections in thin resolution shells. Since the twin related reflections have the same sin(theta)/lambda values they will not be split over the working and reference sets. DATAMAN from the Uppsala Software Factory and XPREP (a program which may be obtained from Bruker) offer this option. The "RFREE SHELL" command in sftools is another way to select thin shells. A disadvantage is the the maps may not be quite as good as when the free R reflections are selected randomly. (FIXME: which Phenix program does this?). A paper investigating this thoroughly is Fabiola, F., A. Korostelev, et al. (2006). "Bias in cross-validated free R factors: mitigation of the effects of non-crystallographic symmetry." Acta Cryst. D 62: 227-38.
  • Sets of reflections used for calculating Rfree should be maintained throughout a project. This is nicely discussed at http://www.bmsc.washington.edu/people/merritt/xplor/rfree_example.html . Note that none of the programs mentioned for selecting thin shells will allow you to extend the set of shells to higher resolution if you want to preserve your existing R-free set.
     twinning  twinning not
     modelled  modelled
twin   0.41      0.49
normal 0.52      0.58

Another paper which investigates the properties of R-values in the presence of twinning is P. R. Evans and G. N. Murshudov (2013) "How good are my data and what is the resolution?" Acta Cryst. (2013). D69, 1204-1214. As the title indicates, this paper discusses at what resolution the data should be cut. One important finding is that a perfect model gives an R value of 42.0% (for a perfect twin, 29.1%) against pure noise. This tells us that a model that gives significantly lower Rfree in the (current) high resolution shell may benefit from including higher resolution data.

  • R-values and pseudo-translation: if you have pseudotranslation you should be aware that if you solve the structure by molecular replacement, starting R factors could be 70-80%.

Notes[edit | edit source]

  1. 1.0 1.1 1.2 1.3 K. Diederichs and P.A. Karplus (1997). Improved R-factors for diffraction data analysis in macromolecular crystallography. Nature Struct. Biol. 4, 269-275 [1]
  2. M.S. Weiss and R. Hilgenfeld (1997) On the use of the merging R-factor as a quality indicator for X-ray data. J. Appl. Crystallogr. 30, 203-205[2]
  3. M.S. Weiss. Global indicators of X-ray data quality. J. Appl. Cryst. (2001). 34, 130-135 [3]
  4. K. Diederichs (2006). Some aspects of quantitative analysis and correction of radiation damage. Acta Cryst D62, 96-101 [4]