diff options
| author | Takanori Nakane <nakane.t@gmail.com> | 2014-08-14 14:58:45 +0100 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2014-09-18 16:50:00 +0200 |
| commit | 1000cd8a44b8943e6c248d16a7ef1a5f49dcb147 (patch) | |
| tree | cfd317d44835a0e0e530bf9e100bf26850939a71 /scripts/plot-cc-and-scale.R | |
| parent | 944c01964d87e9f3830a5079a78e0e5d4903a12e (diff) | |
Added scripts/plot-cc-and-scale.R to show how to plot
frame-wise statistics from process_hkl.
Diffstat (limited to 'scripts/plot-cc-and-scale.R')
| -rw-r--r-- | scripts/plot-cc-and-scale.R | 23 |
1 files changed, 23 insertions, 0 deletions
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 |