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", measure_type != 'd_client_server', measure_type != 'd_server_client', measure_type != 'd_server_client_abs', measure_type != 'd_client_server_abs', run == paste("./out/nBV7DWLpj8WA6Be3-",run_id,sep=""))
  ggplot(data = sel_owd, aes(x = packet_id, y=latency_ms, color=measure_type)) +
    coord_cartesian(ylim=c(0,600)) +
    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()

ggplot(data = sqldf("select * from tor_owd where measure_type='d_client_server_abs' or measure_type='d_server_client_abs'"), aes(x=latency_ms)) +
  stat_ecdf(pad = FALSE) +
  coord_cartesian(xlim = c(0,500)) +
  theme_classic()

ggplot(data = sqldf("select * from tor_owd where packet_id < 500 and run = './out/nBV7DWLpj8WA6Be3-16' and (measure_type='d_client_server' or measure_type='d_server_client')"), aes(y=latency_ms,x=packet_id,color=measure_type)) +
  geom_line() +
  coord_cartesian(ylim=c(-250,250)) +
  geom_point() +
  theme_classic()