CPHASEMATCH (CCP4: Supported Program)

NAME

cphasematch - Calculate agreement between phase sets, allowing for possible origin differences

SYNOPSIS

cphasematch -mtzin filename -mtzout filename -colin-fo colpath -colin-hl-1 colpath -colin-phifom-1 colpath -colin-fc-1 colpath -colin-hl-2 colpath -colin-phifom-2 colpath -colin-fc-2 colpath -colout colpath -resolution-bins number-of-bins -no-origin-hand [-stdin]
[Keyworded input]

DESCRIPTION

'cphasematch' will calculate the agreement between two phase sets, represented either by a phase / figure of merit pair, or as a set of Hendrickson Lattman coefficients. By default, cphasematch will check for possible origin shifts between the two sets.

INPUT/OUTPUT FILES

-mtzin
Input MTZ file containing the phase sets to be compared. These can be represented as a phase / figure of merit pair, or as a set of Hendrickson Lattman coefficients. These columns can be created from a set of atomic coordinates (i.e. a PDB file) by running csfcalc.
-mtzout
Output MTZ file containing, in addition to the input columns, phase information for the second phase set (HL coefficients, phi/fom, Fc depending on what was supplied) adjusted for any origin shift which has been applied.

KEYWORDED INPUT

See Note on keyword input.

-colin-fo colpath

Specify the Fobs column in the input MTZ file, e.g. "*/*/[FP,SIGFP]"

-colin-hl-1 colpath

Specify the Hendrickson Lattman coefficients for the first phase set, e.g. "*/*/[nat.ABCD.A,nat.ABCD.B,nat.ABCD.C,nat.ABCD.D]"

-colin-phifom-1 colpath

As an alternative to -colin-hl-1, specify the phase/FOM for the first phase set, e.g. "*/*/[nat.Phi_fom.phi,nat.Phi_fom.fom]"

-colin-fc-1 colpath

As an alternative to -colin-hl-1 or -colin-phifom-1, specify an Fcalc/phi pair for the first phase set, e.g. "*/*/[nat.F_phi.F,nat.F_phi.phi]". The associated figure of merit is set to 0.9998.

-colin-hl-2 colpath

Specify the Hendrickson Lattman coefficients for the second phase set.

-colin-phifom-2 colpath

As an alternative to -colin-hl-2, specify the phase/FOM for the second phase set.

-colin-fc-2 colpath

As an alternative to -colin-hl-2 or -colin-phifom-2, specify an Fcalc/phi pair for the second phase set.

-colout colpath

Default: "phasematch"
Specify the root of the output MTZ columns.

-resolution-bins number-of-bins

Default: 12
Specify the number of resolution bins for output statistics.

-no-origin-hand

By default, cphasematch will check for possible origin shifts, and changes of hand. Use this option to turn off this behaviour.

Note on keyword input:

Keywords may appear on the command line, or by specifying the '-stdin' flag, on standard input. In the latter case, one keyword is given per line and the '-' is optional, and the rest of the line is the argument of that keyword if required, so quoting is not used in this case. End input with <Ctrl-D> (END is not recognised).

Program Output

The program outputs phase statistics as a function of resolution, and as a function of FOM. In general, figures of merit (FOM) should lie between 0 and 1, with higher values being better. Correlation coefficients should lie between 1 (perfect correlation) and -1 (perfect anti-correlation).

The following graphs are output:

Finally, the program gives overall average values for several quantities, calculated over the full dataset, including:

<dphi>
mean phase difference
w1<dphi>
mean phase difference weighted by the FOM of the first set
wFcorr
correlation coefficient between the structure factor amplitudes of the first and second sets, weighted by the cosine of the phase difference
Qfom1
Average ratio between the cosines of the phase difference and the phase error (latter from the FOM).

Problems

None known.

AUTHOR

Kevin Cowtan, York.

SEE ALSO

phistats