[buster-discuss] Error for rms gradient in Buster or Gelly after sanity check
Thomas Womack
twomack at globalphasing.com
Wed Mar 3 16:55:52 CET 2010
Dear buster users,
The message
' *** ERROR rms gradient is NaN - not a number'
is almost always caused by having two atoms in exactly the same place, often by inserting a water into a peak in coot twice.
The next release will regard this as a sanity error, with message
1 dist= 0.00 idealD= 2.74 sigma= 0.20 abs(Delta/sigma)= 13.70
Between atoms A|121:O (WAT) A|150:O (WAT)
distance is exactly zero. So atoms are exactly coincident.
Advice: check this carefully
Probably the easiest way to figure out which atoms are on top of each other is to use the ccp4 remark-500 tool: run 'r500 HKLIN my.pdb', which will produce a file 'rem500' containing a message of the form
REMARK 500 THE FOLLOWING ATOMS ARE IN CLOSE CONTACT.
REMARK 500
REMARK 500 ATM1 RES C SSEQI ATM2 RES C SSEQI DISTANCE
REMARK 500 O HOH A 121 - O HOH A 150 0.00
Hope this helps,
Tom Womack (Global Phasing)
More information about the buster-discuss
mailing list