Dear Daniele, you could 'misuse' RSTATS:
cad hklin1 data1.mtz hklin2 data2.mtz hklout cad.mtz <<e LABI FILE 1 E1=DANO E2=SIGDANO LABO FILE 1 E1=DANO1 E2=SIGDANO2 LABI FILE 2 E1=DANO E2=SIGDANO LABO FILE 2 E1=DANO2 E2=SIGDANO2 e rstats hklin cad.mtz <<e LABI FP=DANO1 SIGFP=SIGDANO1 FC=DANO2 SIGFC=SIGDANO2 NOAB CYCL 0 OUTP NOHKL END e _should_ work and give you those tables ... Or use SFTOOLS on that cad.mtz with the "CORREL" command. Cheers Clemens On Wed, Oct 01, 2008 at 09:30:32AM +0200, Daniele de Sanctis wrote: > Dear all, > > I'm analyzing some SAD data sets and I would like to compare them in > terms of anomalous correlation. Is there a way in scala (or in some > other program) to calculate a correlation plot like the > file_correlplot.xmgr using two different data sets (or better two > different runs) instead of having scala using the randomly generated > two-half data sets ? > > Cheers > > Daniele > > > > -- > Daniele de Sanctis, PhD > > Macromolecular Crystallography Group > ESRF, Grenoble, France > Tel 33 (0)4 76 88 2869 -- *************************************************************** * Clemens Vonrhein, Ph.D. vonrhein AT GlobalPhasing DOT com * * Global Phasing Ltd. * Sheraton House, Castle Park * Cambridge CB3 0AX, UK *-------------------------------------------------------------- * BUSTER Development Group (http://www.globalphasing.com) ***************************************************************