I'm not sure this is a particularly rigorous way to get "noisier bits" less noisy. Surely the way to do it is to refine one model using the other model as restraints. Both refmac and buster do it, (maybe phenix.refine too). Cheers phx On 10/03/2011 17:15, Stephen Graham wrote:
Hi there,
I'd like to average 2Fo-Fc maps from 2 (approximately) isomorphous data sets (one with ligand, one without) to see whether I get better density for some of the noisier bits of the structure. I don't want to do any kind of solvent flattening, histogram matching, cyclical averaging, etc; I just want to sum the equivalent pixels and divide by 2.
Is it possible to do this easily using any of the phenix tools? It seems like all the required plumbing is there with get_cc_mtz_mtz or superpose_maps, but I can't find a phenix utility for averaging two superposed maps.
Thanks,
Stephen