From 1000cd8a44b8943e6c248d16a7ef1a5f49dcb147 Mon Sep 17 00:00:00 2001
From: Takanori Nakane <nakane.t@gmail.com>
Date: Thu, 14 Aug 2014 14:58:45 +0100
Subject: Added scripts/plot-cc-and-scale.R to show how to plot frame-wise
 statistics from process_hkl.

---
 scripts/plot-cc-and-scale.R | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)
 create mode 100644 scripts/plot-cc-and-scale.R

(limited to 'scripts')

diff --git a/scripts/plot-cc-and-scale.R b/scripts/plot-cc-and-scale.R
new file mode 100644
index 00000000..75850e85
--- /dev/null
+++ b/scripts/plot-cc-and-scale.R
@@ -0,0 +1,23 @@
+# Example R script to plot scale factor and CC
+#
+#      by Takanori Nakane (nakane.t@gmail.com)
+
+scale_cc <- read.table("scale.csv", sep=",", col.names=c("ImageName", "ScaleFactor", "CC"))
+
+# scatter plot
+# FIXME: sometimes white lines appear on the plot.
+#        seems to be a bug in smoothScatter
+
+smoothScatter(scale_cc$ScaleFactor, scale_cc$CC,
+              nbin=200, xlim=c(0, 2), ylim=c(0, 1),
+              xlab="Scale Factor", ylab="CC", main="Scale factor & CC")
+
+# histogram
+
+par(mfrow=c(2,1))
+plot(scale_cc$ScaleFactor, type='h', ylab="Scale Factor", xlab="#Crystal")
+abline(h=mean(sf[,2])+(-5:5)*sd(sf[,2]), col=2)
+
+hist(scale_cc$ScaleFactor, breaks=100, main="Distribution of Scale Factor", xlab="Scale Factor")
+abline(v=mean(sf[,2])+(-5:5)*sd(sf[,2]), col=2)
+par(mfrow=c(1,1))
\ No newline at end of file
-- 
cgit v1.2.3