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

thunder_ms <- read.csv("thunder_configure_2_ms.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() +
  ylim(0,1000) +
  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")