James, I must admit not (yet) having read the "data processing paper PMID: 20389587" nor the 88 author paper. Nevertheless, I would like to comment on your remark "Personally, I was quite impressed by how good the R factors were, all things considered."
If I am not mistaken, perfectly twinned intensity data, such as you seem to have been dealing with, will generate data sets with compressed dynamic ranges, which in turn means that crystallographic R factors computed based on such data sets will be lower by roughly sqrt(2) or 1.4 compared to a non-twinned data set. So a somewhat iffy R(cryst) of 30% would look quite nice (21.4%) when computed on a hemihedrally twinned data set. But maybe this is all discussed in your paper... Cheers, Hudel