Dear Kay, > Concerning usage of programs, everybody has his/her preferences, but what > could be simpler than a 2-liner XDSCONV.INP like > INPUT_FILE=XDS_ASCII.HKL > OUTPUT_FILE=temp.hkl CCP4 ! or CCP4_F or CCP4_I or SHELX or CNS > and then running XDSCONV by running "xdsconv"? At least there's not much room > for mistakes.
Exactly! >> I beg to disagree, and this was the point of my question. The output of >> XDSCONV literally says that 190093 reflections are read, and [of those, my >> interpretation] 44047 are accepted. I may be a pedant but I can't read that >> output in any other way. To me it looks like it is reading only the >> reflections that already fall into the asymmetric unit and is ignoring all >> the others. So if XDSCONV is really doing what it is supposed to, I would >> suggest rephrasing that output line. > > good point about the rephrasing. I'll see to making the wording consistent > between XDSCONV and XDS. Thanks for that. > But irrespective of the wording, it does take all observations into account > when calculating the intensity (and amplitude) of the unique reflections (as > it should - ignoring reflections would not make sense, and would produce > significantly worse data). Great, this was just the reassurance I was looking for. > However, it does so not by calculating the geometric mean (which you seem to > assume), but by calculating the weighted mean. Weighting is done with the > variances, and here it also does not differ from (c)truncate or other > programs. Sorry, that was a typo born of thinking I could quote the manual from memory. >> But yes, pedantry aside, I will start using pointless in my csh pipelines. > > please report back whether that changes (or even improves) your results! It > is always very good when people compare programs in a meaningful way, but > from my own experience I can say that meaningful comparisons are sometimes > not entirely straightforward to get right. (for the German-speaking: "wer > misst, misst Mist!") Probably I will still be in too much of a hurry each time to make a meaningful and thorough comparison, but if time allows I will compare XDSCONV with pointless/scala. /Derek