International
Tables for Crystallography Volume F Crystallography of biological macromolecules Edited by M. G. Rossmann and E. Arnold © International Union of Crystallography 2006 |
International Tables for Crystallography (2006). Vol. F. ch. 11.4, p. 233
Section 11.4.8. Multiplicative corrections (scaling)
a
UT Southwestern Medical Center at Dallas, 5323 Harry Hines Boulevard, Dallas, TX 75390-9038, USA, and bDepartment of Molecular Physiology and Biological Physics, University of Virginia, 1300 Jefferson Park Avenue, Charlottesville, VA 22908, USA |
Proper error estimation requires the use of Bayesian reasoning and a multi-component error model (Schwarzenbach et al., 1989; Evans, 1993
). In SCALEPACK, the estimated error of the measurement is enlarged by a fraction of the expected, rather than the observed, intensity. This algorithm reduces the bias towards reflections with an integrated intensity below the average.
The scaling model allows for a large number of diverse components to contribute to the multiplicative correction factor s for each observation, where
are a priori unknown coefficients of the scaling components and
represent different functional dependence of the scale factor for each observation. The simplest scaling model has a separate scale factor for each group (batch) of data, for example, one scale factor per image. In such a case,
where j is the batch index for a particular reflection. For resolution-dependent decay, represented by one temperature factor per batch of data,
where n is the number of batches needed to make
represent the logarithm of the overall ith batch scale factor and
represents the temperature factor of batch i. S is the scattering vector for each reflection.
Coefficients of crystal absorption are much more complex. Scaling coefficients are associated with spherical harmonics (Katayama, 1986; Blessing, 1995
) as a function of the direction of vector S, expressed in the coordinate system of the rotating crystal. Each spherical harmonic index lm has two (or, in the case of
, one)
coefficients. One of these spherical harmonic functions is given by
where Ψ and Φ are polar-coordinate angles of vector S in the crystal coordinate system, and
is a Legendre polynomial (Press et al., 1989
). The other spherical harmonic of index lm has a sine instead of a cosine as the last term.
The multiplicative factor is applied to each observation and its σ to obtain the corrected intensity and associated σ. The averaged intensity over all symmetry-related reflections
is obtained by solving the two following equations:
where
and
are the user-specified error scale factor and estimated error, respectively, and
Thus,
During parameter refinement, the scale (and B, if requested) for all scaled batches are refined against the difference between the
's and
's for individual measurements, summed over all reflections (Fox & Holmes, 1966
; Arndt & Wonacott, 1977
; Stuart & Walker, 1979
; Leslie & Tsukihara, 1980
; Rossmann & Erickson, 1983
; Walker & Stuart, 1983
; Rossmann, 1984
; Schutt & Evans, 1985
; Stuart, 1987
; Takusagawa, 1987
; Tanaka et al., 1990
).
's are recalculated in each cycle of refinement. There is full flexibility in the treatment of anomalous pairs. They can be assumed to be equivalent (or not) and they may be merged (or not). This approach allows the crystallographer to choose the best scaling and merging strategy.
A polarization correction may be applied during DENZO calculations. Sometimes the exact value of polarization is not known. This error may be corrected during the scaling procedure. This feature can be used to refine the polarization at synchrotron beamlines. Very high resolution data should be used for this purpose.
References















