[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

*To*: "Patrick J. Loll" <loll@pharm.med.upenn.edu>*Subject*: Re: [ccp4bb]: Phasing Phoiled by Phistats in Philadelphia. Pfui.*From*: Bart Hazes <bhazes@ualberta.ca>*Date*: Thu, 8 Feb 2001 15:34:40 -0700 (MST)*cc*: ccp4bb@dl.ac.uk*In-Reply-To*: <l03130306b6a8c13cda2b@[130.91.170.247]>*Sender*: owner-ccp4bb@dl.ac.uk

*** For details on how to be removed from this list visit the *** *** CCP4 home page http://www.dl.ac.uk/CCP/CCP4/main.html *** SFTOOLS has an option specially for this situation. It is called PHASHFT and you only have to tell it which two columns to compare. It can print the phase difference in either degrees or as a figure of merit and it will give the stats as function of resolution. E.g. (assuming the phase columns have labels phi1 and phi2): sftools yourfile.mtz phashft col phi1 phi2 stop Phoila... Bart On Thu, 8 Feb 2001, Patrick J. Loll wrote: > *** For details on how to be removed from this list visit the *** > *** CCP4 home page http://www.dl.ac.uk/CCP/CCP4/main.html *** > > I'm confused about output I obtain from phistats. I'm attempting to > calculate the phase errors associated with various phasing procedures. > I've written a little program to calculate the quantity: > > <|deltaPhi|> = the mean value over all reflxns of the absolute value of the > difference between Phi1 and Phi2. I calculate this for all > reflxns and for acentrics only > > In all cases, the DELF number obtained from phistats (described in the > output as "average difference in degrees" for acentric reflxns) is > substantially lower than the number I calculate independently (by, say, 20 > degrees). > > I've tried all the controls I can think of. The most stringent control was > to write h, k, l, F, sigF, Phi1, & Phi2 into a single ASCII file, one line > per reflxn. I made certain Phi1 and Phi1 both lay between 0 and 359 deg > (ie, not -180 to 180). I can then use F2mtz to create an mtz file from > this & feed it into phistats. This should cover most of the obvious > problems (eg, reflxns from different sources reduced to different asym > units). The problem still persists. > > I'm assuming DELF is equivalent to quantity I describe above. Am I wrong? > I think it should mean the average absolute difference, rather than the > "average difference," since the phase errors are uncorrelated, and over a > large dataset will average to zero. > > Of course one might suspect the little program I've written; but I've > checked it every way I can imagine; and I think that even my meager coding > skills are up to the task of calculating a difference between two columns > of numbers and then averaging it... > > > > Pat, Puzzled in Pennsylvania > > ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ > > Patrick J. Loll, Ph. D. (215) 898-1294 > Department of Pharmacology FAX: (215) 573-2236 > University of Pennsylvania School of Medicine > 3620 Hamilton Walk > Philadelphia, PA 19104-6084 USA > > loll@pharm.med.upenn.edu > http://pharm.med.upenn.edu/Loll/loll_lab.html > > =============================================================================== Dept. of Medical Microbiology & Immunology University of Alberta 1-15 Medical Sciences Building Edmonton, Alberta, T6G 2H7, Canada phone: 1-780-492-0042 fax: 1-780-492-7521 ===============================================================================

**References**:**[ccp4bb]: Phasing Phoiled by Phistats in Philadelphia. Pfui.***From:*"Patrick J. Loll" <loll@pharm.med.upenn.edu>

- Prev by Date:
**[ccp4bb]: Phasing Phoiled by Phistats in Philadelphia. Pfui.** - Next by Date:
**Re: [ccp4bb]: Phasing Phoiled by Phistats in Philadelphia. Pfui.** - Prev by thread:
**[ccp4bb]: Phasing Phoiled by Phistats in Philadelphia. Pfui.** - Next by thread:
**PHabulously-PHriendly PHrequently-PHrostbitten PHortran PHringe-PHossil PHinds PHormidable PHrigging PHrightening PHrivolous PHundamental PHunctional PHormula PHalsehood (Re: [ccp4bb]: Phasing Phoiled by Phistats in Philadelphia. Pfui.)** - Index(es):