CRITIQ: Determine the Most Variable Reflections
Authors: Richard Alden, Jim Stewart and Keith Watenpaugh
Contact: Keith Watenpaugh, Physical and Analytical Chemistry,
The Upjohn Company, Kalamazoo, MI49001, USA.
CRITIQ is used to examine multiply observed intensity data and to pick out
the most variable reflections. A bdf containing replicated observations is read
and a list of the reflections whose multiple observations exhibit the most
variability is produced. The N most variable reflections according to each of
four categories are listed. N is user specified. It is also possible to list
all singly-observed reflections within a specified range of
sinθ/λ.
Overview
Five different discrepancy lists are prepared; a separate compact list, in
order of decreasing discrepancy, for each of the four ranking categories (see
below), plus a merged list containing all ranked reflection intensities. The
merged list is listed in reflection sort order. To aid editing, the merged list
contains, in additon to the ranked observations, up to a specified maximum
number of unranked observations for each reflection.
The ranking categories are as follows:
1. |K*I-Y|
2. (1-Y/(K*I))2 (corresponds to the SCALE1 R-factor)
3. |K*I-Y| / σY
4. |K*I-Y| / σI
where I is the observed intensity(rel), Y is the mean intensity
of a particular reflection and K is the scale constant for a particular scale
group. For observations where I < 0.1, the value for category number 2 is
not computed (i.e., the observation is not considered for ranking purposes in
that category).
To clarify what is meant by "decreasing discrepancy" for a category, consider
the observation which has the largest absolute value for the expression, K*I-Y.
It will be ranked number one in category 1 above and will appear as such in the
listing for that category. Of course it may also appear in one or more of the
remaining three categories. This observation will also be listed in the merged
set of all discrepancies under the reflection to which it belongs. All other
observations for that reflection which have rankings in any of the four
categories will also appear, as well as a user-specified maximum of unranked
observations (if any), still for that same reflection.
File Assignments
Reads reflection data from the input archive bdf
Example
title ANALYZE NATIVE DATA AFTER INITIAL SCALING
CRITIQ 1 50
The listing will consist of the 50 largest discrepancies in each of the 4
categories with all scale group numbers and all data being involved. No singly
observed reflections will be listed and Bijvoet pairs will not be considered as
separate reflections. All shift numbers for which there are scale factors will
be read from logical record lrrefl: of the input bdf.