diff --git a/analysis/final/R/plots.R b/analysis/final/R/plots.R
new file mode 100644
index 0000000..704908a
--- /dev/null
+++ b/analysis/final/R/plots.R
@@ -0,0 +1,28 @@
+comparison_plot <- function(data, p25, p75, median, ylab, free_y = FALSE) {
+  ggplot(data, aes(x = network_size, col = label, fill = label, group = label)) +
+    geom_ribbon(
+      aes(ymin = {{ p25 }}, ymax = {{ p75 }}),
+      col = 'lightgray', alpha = 0.5
+    ) +
+    geom_point(aes(y = {{ p25 }}), col = 'darkgray', size = 10, shape = '-') +
+    geom_point(aes(y = {{ p75 }}), col = 'darkgray', size = 10, shape = '-') +
+    geom_line(aes(y = {{ median }})) +
+    geom_point(aes(y = {{ median }})) +
+    ylab(ylab) +
+    xlab('network size') +
+    theme_minimal(base_size = 15) +
+    facet_grid(
+      scales = if (free_y) 'free_y' else 'fixed',
+      file_size ~ seeder_ratio,
+      labeller = labeller(
+        seeder_ratio = as_labeller(function(x) {
+          paste0("seeder ratio: ", scales::percent(as.numeric(x)))
+        }),
+      )
+    ) +
+    scale_color_discrete(name = '') +
+    guides(fill = 'none', alpha = 'none')
+}
+
+Y_BPS <- scale_y_continuous(labels = function(x) paste0(scales::label_bytes()(x), '/s'))
+Y_TIMESPAN <- scale_y_continuous(labels = scales::label_timespan())