From 6264ff6e0ac13fcc22afcb4ce9f98a920444a92c Mon Sep 17 00:00:00 2001 From: Thomas White Date: Mon, 12 Jul 2010 10:06:47 +0200 Subject: compare_hkl: Add Pearson correlation --- src/compare_hkl.c | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) (limited to 'src/compare_hkl.c') diff --git a/src/compare_hkl.c b/src/compare_hkl.c index 3e313733..0dffd3a0 100644 --- a/src/compare_hkl.c +++ b/src/compare_hkl.c @@ -50,7 +50,7 @@ int main(int argc, char *argv[]) char *afile = NULL; char *bfile = NULL; signed int h, k, l; - double scale, R2, Rmerge; + double scale, R2, Rmerge, pearson; unsigned int *c1; unsigned int *c2; int i; @@ -145,6 +145,8 @@ int main(int argc, char *argv[]) STATUS("R2 = %5.4f %% (scale=%5.2e)\n", R2*100.0, scale); Rmerge = stat_rmerge(ref1, c1, ref2, c2, &scale); STATUS("Rmerge = %5.4f %% (scale=%5.2e)\n", Rmerge*100.0, scale); + pearson = stat_pearson(ref1, c1, ref2, c2); + STATUS("Pearson r = %5.4f\n", pearson); if ( outfile != NULL ) { write_reflections(outfile, NULL, out, NULL, 1, cell, 1); -- cgit v1.2.3