1 args <- commandArgs(TRUE)
8 d1 <- read.csv(file1, header=F)
9 d2 <- read.csv(file2, header=F)
11 sink(file="summary1.txt", split=T)
15 con <- file("summary1.txt")
17 summary_text <- readLines(con)
20 sink(file="summary2.txt", split=T)
24 con <- file("summary2.txt")
26 summary_text <- readLines(con)
30 png(sprintf("%s.png", file1))
33 par(mar = c(4.5, 4.5, 4.5, 5.5))
37 truehist(d1$V1, xlim=c(0,max), col="#66ff6640", border="#66ff66", axes = FALSE, xlab = "", ylab = "", prob=FALSE)
38 axis(side = 1) # x axis
39 axis(side = 2, col.axis = "#66ff66", col = "#66ff66") # left y axis
40 mtext("floodgate-900-0", side = 2, line=3)
44 truehist(d2$V1, xlim=c(0,max), col="#6699ff40", border="#6699ff", axes = FALSE, xlab = "", ylab = "", prob=FALSE)
45 axis(side = 4, col.axis = "#6699ff", col = "#6699ff") # right y axis
46 mtext("floodgate-600-0", side = 4, line=3)
48 mtext(title, side=1, line=3)