library(ggplot2)
library(sqldf)
library(plyr)
library(cowplot)

thunder_ms <- read.csv("thunder_configure_10.csv")

thunder_ms <- sqldf("select run,jmax,links,latency, CAST(latency as real) / 1000. as lat_ms from thunder_ms")
thunder_ms$links <- as.factor(thunder_ms$links)
thunder_ms$jmax <- as.factor(thunder_ms$jmax)

v1 <- ggplot(data = thunder_ms, aes(x = jmax, y=lat_ms, fill=links)) +
  geom_boxplot() +
  scale_fill_grey() +
  #scale_y_log10() +
  ylab("latency (ms)") +
  xlab("max allowed jitter") +
  theme_classic()

v1 + ggsave("thunder_configure_ms.png", dpi=300, dev='png', height=5, width=15, units="cm")

thunder_rcv <- sqldf("select run,jmax,links,(6000-COUNT(latency)*1.0)/6000 as dlv from thunder_ms group by jmax,links,run")
thunder_rcv$jmax <- as.factor(thunder_rcv$jmax)
thunder_rcv$links <- as.factor(thunder_rcv$links)
v2 <- ggplot(data = thunder_rcv, aes(x = jmax, y=dlv, fill=links)) +
  geom_boxplot() +
  scale_y_continuous(labels = scales::percent) +
  scale_fill_grey() +
  ylab("dropped packets") +
  xlab("max allowed jitter") +
  theme_classic()

thunder_bw <- read.csv("thunder_configure_10_bw.csv")
thunder_bw <- sqldf("select run,jmax,links,udp_sent,udp_rcv,cells_sent,cells_rcv,1.0*cells_sent/udp_sent as sent_ratio,1.0*cells_rcv/udp_rcv as rcv_ratio from thunder_bw where udp_sent > 4000")
thunder_bw$jmax <- as.factor(thunder_bw$jmax)
thunder_bw$links <- as.factor(thunder_bw$links)

v3 <- ggplot(data = thunder_bw, aes(x = jmax, y=sent_ratio, fill=links)) +
  geom_boxplot() +
  #scale_y_log10() +
  #scale_y_log10(labels = scales::percent) +
  scale_fill_grey() +
  ylab("bandwidth ratio") +
  xlab("max allowed jitter") +
  theme_classic()


t1 <- plot_grid(v1, v2, v3, labels = c('A', 'B', 'C'), ncol=1)
t1 + ggsave("thunder_configure.png", dpi=300, dev='png', height=15, width=15, units="cm")