aboutsummaryrefslogtreecommitdiff
path: root/scripts/plot-cc-and-scale.R
diff options
context:
space:
mode:
authorTakanori Nakane <nakane.t@gmail.com>2014-09-19 11:31:29 +0100
committerThomas White <taw@physics.org>2014-09-19 13:07:15 +0200
commit2c959daa7a46b99a10dd5a1998b62ccb8def97de (patch)
tree5037024a60320f18825b0834858e87b4fdd88db1 /scripts/plot-cc-and-scale.R
parent0a533227b1006b4019ea9fccd0957e8be7138f8d (diff)
Fix variable name in plot-cc-and-scale.R
Plot histogram of CC. Conflicts: scripts/plot-cc-and-scale.R
Diffstat (limited to 'scripts/plot-cc-and-scale.R')
-rw-r--r--scripts/plot-cc-and-scale.R18
1 files changed, 13 insertions, 5 deletions
diff --git a/scripts/plot-cc-and-scale.R b/scripts/plot-cc-and-scale.R
index 25df0a63..b67ab058 100644
--- a/scripts/plot-cc-and-scale.R
+++ b/scripts/plot-cc-and-scale.R
@@ -2,7 +2,7 @@
#
# by Takanori Nakane (nakane.t@gmail.com)
-scale_cc <- read.table("stats.dat", sep=" ", col.names=c("ImageName", "ScaleFactor", "CC"))
+scale_cc <- read.table("stats.dat", col.names=c("ImageName", "ScaleFactor", "CC"))
# scatter plot
# FIXME: sometimes white lines appear on the plot.
@@ -14,10 +14,18 @@ smoothScatter(scale_cc$ScaleFactor, scale_cc$CC,
# histogram
-par(mfrow=c(2,1))
+par(mfrow=c(2,2))
plot(scale_cc$ScaleFactor, type='h', ylab="Scale Factor", xlab="#Crystal")
-abline(h=mean(sf[,2])+(-5:5)*sd(sf[,2]), col=2)
+abline(h=mean(scale_cc$ScaleFactor) + (-5:5) * sd(scale_cc$ScaleFactor), 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)
+hist(scale_cc$ScaleFactor, breaks=100, main="Distribution of Scale Factor",
+ xlab="Scale Factor")
+abline(v=mean(scale_cc$ScaleFactor) + (-5:5) * sd(scale_cc$ScaleFactor), col=2)
+
+plot(scale_cc$CC, type='h', ylab="CC", xlab="#Crystal")
+abline(h=mean(scale_cc$CC) + (-5:5) * sd(scale_cc$CC), col=2)
+
+hist(scale_cc$CC, breaks=100,
+ main="Distribution of Correlation Coefficient", xlab="CC")
+abline(v=mean(scale_cc$CC) + (-5:5) * sd(scale_cc$CC), col=2)
par(mfrow=c(1,1))