next up previous contents
Next:  My native Patterson function Up:  Pathology of GraphEnt calculations, Previous:  The GraphEnt map changes   Contents


 All my anomalous Pattersons are ``consistent with a uniform map''

Why : Assume for a minute that all your reflections have the same F/$ \sigma$(F) = $ \mathcal {C}$. Then, because for the initial (uniform) map all Fourier coefficients are identically zero (with the exception of F000), the statistic $ \chi^{2}_{}$ is given by

$\displaystyle \chi^{2}_{}$ = $\displaystyle \sum_{\bf h}^{}$$\displaystyle {\frac{\mid F_{c, {\bf h}} - F_{o, {\bf h}} \mid ^2}{\sigma(F_{\bf h})^2}}$ = $\displaystyle \sum_{\bf h}^{}$$\displaystyle {\frac{\mid F_{o, {\bf h}} \mid ^2}{\sigma(F_{\bf h})^2}}$ = $\displaystyle \sum_{\bf h}^{}$$\displaystyle \mathcal {C}$2 = N$\displaystyle \mathcal {C}$2

where N is the number of observed reflections. Now, the target value for $ \chi^{2}_{}$ during the calculation is $ \chi^{2}_{}$ = N. If $ \mathcal {C}$ < 1.0, ie $ \overline{F/\sigma(F)}$ < 1.0, the uniform map will satisfy the $ \chi^{2}_{}$ constraint and GraphEnt will stop immediately. Just because GraphEnt stops, does not necessarily mean that there is no signal in the data : for Gaussian noise, $ \chi^{2}_{}$ $ \approx$ N is the expected value of the distribution $ \chi^{2}_{}$ $ \approx$ N$ \pm$$ \sqrt{(2N)}$. This means that depending on the data, the target of the calculation could as well be significantly lower than the value aimed for by GraphEnt. I think it is worth emphasising this with an example. The following figure compares the conventional (right) and GraphEnt (left) map at the section v = 1/2 of the 20-3Å anomalous Patterson function for horse heart myoglobin crystals (dashed contour at the mean, and then every 0.5 rmsd of the whole map). The data were collected with CuK$\scriptstyle \alpha$ radiation and the anomalous signal comes from the iron atom of heme (this is one of the examples distributed with GraphEnt, file Myoglobin_anom_Patt_no_outliers.in).

\scalebox {0.9}{\includegraphics{/usr/people/glykos/progs/maxent/doc/fig/myo.ps}}

A normal run of GraphEnt with the whole data set would immediately stop with the ``uniformity'' message. Even after rejecting all reflections with F/$ \sigma$(F) < 0.5, GraphEnt would still refuse to co-operate (for 615 reflections with F/$ \sigma$(F) > 0.5, the initial $ \chi^{2}_{}$ --for the uniform map-- was 502.8). The map shown above could only be produced after explicitly setting the TARGet $ \chi^{2}_{}$-value to 100.0 (by editing the MAXENT_AUTO.IN file). As you see, it probably worth the effort16.

Getting around it : Start GraphEnt the usual way. When the program stops with the ``uniformity'' message, edit the MAXENT_AUTO.IN and add a line with a new TARGet value (which should be less than the starting $ \chi^{2}_{}$ value reported by the program if the VERBose flag is set on, see page [*]). Depending on the circumstances, you could also add a line with LIMIt 0.5 to exclude reflections with F/$ \sigma$(F) < 0.5 (this should reduce the amount of computation required for convergence).



Footnotes

... effort16
I should add, however, that I am not convinced that changing the TARGet $ \chi^{2}_{}$ is the correct way around. To continue with the example, even if we take the expected value of $ \chi^{2}_{}$ to be $ \chi^{2}_{}$ $ \approx$ N - 3$ \sqrt{(2N)}$ = 509 (ie 3$ \sigma$ away from the mean), the uniform map is still consistent with the data. The fact that there seems to be some signal in the data when we reduce the TARGet, probably points the way to over-estimated standard deviations.

next up previous contents
Next:  My native Patterson function Up:  Pathology of GraphEnt calculations, Previous:  The GraphEnt map changes   Contents
NMG, March 2000