ENH: Optimize plots for publication

Also refactor ggsave() export to multiple formats.

.gitignore

@@ -6,6 +6,7 @@ z_stack/
 *.csv
 *.pdf
 *.png
+*.svg
 
 # Log files from CellProfiler Analyst
 *.log

plots.R

@@ -11,28 +11,33 @@ binwidth <- 0.025
 ## Bring RNAi to the front of the level list.
 wells_all$type <- relevel(wells_all$type, "RNAi")
 
+theme_set(theme_bw())
+
+plot_save <- function(prefix, extensions = c("png", "pdf", "svg")) {
+    invisible(sapply(paste(prefix, extensions, sep = "."),
+                     ggplot2::ggsave, width = 8, height = 5,
+                     units = "in"))
+}
+
 ggplot(filter(wells_all, type == "RNAi"),
        aes(coloc, fill = significant)) +
     geom_histogram(binwidth = binwidth) +
-    scale_fill_manual(values = c("grey70", "blue")) +
-    ggtitle("Histogram of RNAi wells with significant overlap (p-value < 5%)") +
+    scale_fill_manual(values = c("grey", "blue")) +
+    ggtitle("Histogram of wells with significant overlap (p-value < 5%)") +
     xlab("% colocalization")
-ggsave("rnai-hist-signif.png")
-ggsave("rnai-hist-signif.pdf")
+plot_save("rnai-hist-signif")
 
 ggplot(wells_all, aes(coloc, color = type)) +
     geom_freqpoly(binwidth = binwidth) +
     scale_y_log10() +
     facet_wrap(~ type) +
     ggtitle("Histogram of all well types") +
     xlab("% colocalization")
-ggsave("all-hist.png")
-ggsave("all-hist.pdf")
+plot_save("all-hist")
 
 ggplot(wells_all, aes(coloc, color = type)) +
     geom_density() +
     facet_wrap(~ type) +
     ggtitle("Density plot of wells to show distribution") +
     xlab("% colocalization")
-ggsave("all-density.png")
-ggsave("all-density.pdf")
+plot_save("all-density")