diff options
-rw-r--r-- | scripts/plot-cc-and-scale.R | 18 |
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)) |