diff --git a/plots.R b/plots.R
new file mode 100755
index 0000000000000000000000000000000000000000..d838069df785b02610289fa6e33cb3e3232f3055
--- /dev/null
+++ b/plots.R
@@ -0,0 +1,56 @@
+#!R -f
+
+(function(lp) {
+np <- lp[!(lp %in% installed.packages()[,"Package"])]
+if(length(np)) install.packages(np,repos=c("https://cran.rstudio.com/"))
+x <- lapply(lp,function(x){library(x,character.only=TRUE)}) 
+})(c("ggplot2", "ggthemes", "gtable", "reshape2"))
+
+theme <- theme_few(base_size = 24) + 
+theme(axis.title.y=element_text(vjust=0.9), 
+  axis.title.x=element_text(vjust=-0.1),
+  axis.ticks.x=element_blank(),
+  text=element_text(family="serif"),
+  legend.position = "none")
+
+dd <- read.table("results.csv", sep="\t", header=F)
+names(dd) <- c("sys", "dataset", "nrows", "loadtime", "querytime", "readtime")
+ddm <- melt(dd[c(1,2,4,5,6)], id.vars=c("sys","dataset"))
+ddm$sys <- factor(ddm$sys, c("MonetDBLite", "MonetDB", "SQLite"))
+
+print(ddm)
+
+xlf <- function(value) ifelse(value == 0, "DNF", ifelse(value < 10, 
+  paste0(round(value, 1), "s"), paste0(round(value), "s")))
+
+plt <- function(ddp, fname, xmax, title) {
+  pdf(fname, width=10, height=5)
+  p <- ggplot(ddp,aes(x=dataset, y=value, fill=sys)) + 
+    geom_bar(stat="identity", position = "dodge", width=.7)  +
+    geom_text(aes(label = xlf(value), family="serif"), size = 5, vjust=-.4, 
+      position = position_dodge(width=.7)) + 
+    scale_x_discrete(labels=xlabel) +
+    scale_y_continuous(limits=c(0, xmax)) + 
+    scale_fill_manual(values=c( "#1f78b4", "#a6cee3", "#b2df8a")) +
+    xlab("# Tuples") + ylab("Time (s)") + theme + ggtitle(title)
+
+  # haaaack (secondary x axis patched in using gtable)
+  axis <- ggplot(ddp,aes(x=dataset, y=value, fill=sys)) + 
+    geom_text(aes(label=sys, y=0, family="serif"), angle=90, size = 5,  hjust = 0.55, 
+      position = position_dodge(width=.7))
+  annotation <- gtable_filter(ggplotGrob(axis), "panel", trim=TRUE)
+  annotation[["grobs"]][[1]][["children"]][c(1,3)] <- NULL
+  g <- ggplotGrob(p)
+  g <- gtable_add_rows(g, unit(6, "line"), pos=3)
+  g <- gtable_add_grob(g, annotation, t=4, b=4, l=4, r=4)
+  grid.draw(g)
+
+  dev.off()
+}
+
+xlabel <- c("1"="60K", "2"="600K", "3"="6M", "4"="60M", "f"="128M")
+plt(ddm[ddm$variable=='loadtime',], "load.pdf", 10000, "Loading from CSV files")
+plt(ddm[ddm$variable=='querytime',], "query.pdf", 70, "Run HMDA analysis")
+
+xlabel <- c("1"="10K", "2"="100K", "3"="1M", "4"="10M", "f"="14M")
+plt(ddm[ddm$variable=='readtime',], "read.pdf", 1500, "Convert table to data.frame")