plot_histograms.R

if (!"ggplot2" %in% rownames(installed.packages())) {
  install.packages("ggplot2")
}

library("ggplot2")

args <- commandArgs(trailingOnly = TRUE)
before_file <- args[1]
after_file <- args[2]

if (!file.exists(before_file)) stop(paste(before_file, "does not exist!"))
if (!file.exists(after_file)) stop(paste(after_file, "does not exist!"))

df1 <- read.table(before_file,  sep="\t", header = TRUE)
df2 <- read.table(after_file,  sep="\t", header = TRUE)

df1$test <- "before"
df2$test <- "after"

df <- rbind(df1, df2)

svg("histogram_interleaved.svg", width = 7, height = 4)

ggplot(df, aes(ttime, fill = test)) +
  scale_x_continuous(breaks = seq(min(df$ttime), max(df$ttime))) +
  geom_histogram(binwidth = 1, position="dodge")

svg("histogram_facet.svg", width = 7, height = 4)

ggplot(df, aes(ttime, fill = test)) +
  geom_histogram(binwidth = 1) +
  facet_wrap(~test, ncol = 1)

dev.off()