library(ggplot2) uptime = data.frame( date=seq(from=as.Date("2013-11-04"), to=as.Date("2014-11-05"), by="2 days"), frac=c(999,957,0,20,999,999,999,999,999,999,999,999,999,999,999,999,999,999,999,999,999,999,520,0,228,874,353,999,999,999,457,0,0,0,0,0,853,999,999,999,749,749,0,457,999,999,999,999,999,999,999,NA,999,999,999,478,0,0,0,185,999,999,999,999,999,999,999,999,999,999,999,999,999,999,999,999,999,999,478,0,0,NA,NA,NA,NA,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,166,208,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,104,395,999,999,999,790,541,0,0,686,999,999,999,999,999,999,999,999,999,999,999,166,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,20,0,0,0,0,66)/999 ) c <- read.csv("clients.csv", colClasses=c(date="Date")) p <- ggplot(c[c$date >= as.Date("2014-04-01") & c$date < as.Date("2014-11-01") & c$node=="bridge" & (c$transport=="obfs3" | c$transport==""), ]) p <- p + geom_rect(data=uptime, aes(xmin=uptime$date, xmax=uptime$date+2, ymin=-Inf, ymax=Inf), fill="lightblue", alpha=uptime$frac) p <- p + geom_line(aes(date, clients, color=transport)) p <- p + xlim(c(as.Date("2014-04-01"), as.Date("2014-11-01"))) p <- p + theme_bw() ggsave("obfs3-or-pacificsunset.png", p, width=4.8, height=3.2, dpi=120)