diff --git a/r/tor_owd.R b/r/tor_owd.R
new file mode 100644
index 0000000..3922804
--- /dev/null
+++ b/r/tor_owd.R
@@ -0,0 +1,20 @@
+library(ggplot2)
+library(sqldf)
+library(dplyr)
+library(cowplot)
+
+tor_owd <- read.csv("tor_owd.csv")
+
+for (run_id in c(0:39)) {
+  print(run_id)
+  sel_owd <- filter(tor_owd, measure_type != "delta", run == paste("./out/nBV7DWLpj8WA6Be3-",run_id,sep=""))
+  ggplot(data = sel_owd, aes(x = packet_id, y=latency_ms, color=measure_type)) +
+    geom_line() +
+    theme_classic() +
+    ggsave(paste("tor_owd_",run_id,".png", sep = ""), dpi=150, dev='png', height=18, width=32, units="cm")
+}
+
+ggplot(data = sqldf("select * from tor_owd where measure_type='delta_abs'"), aes(x=latency_ms)) +
+  stat_ecdf(pad = FALSE) +
+  coord_cartesian(xlim = c(0,500)) +
+  theme_classic()