plots output of https://gist.github.com/mtholder/81a6378d668dc87204515005b36e8c59

plot-ot-ts.R

d = read.table('times_num_studies.tsv', header=TRUE, sep="\t");
print(d);
dt = as.POSIXct(d$time, origin="1970-01-01");
print(dt);

pdf("ot-time-series.pdf");
plot(dt, d$num_phyle_studies,
     xlab="Date", ylab="num studies",
     ylim=c(0, max(d$num_phyle_studies)),
     type="l",
     col=1);
lines(dt, d$num_synth_trees, col=2);
legend(min(dt),
       max(d$num_phyle_studies),
       c("total # studies in phylesystem", "# trees used in synth"),
       lty=c(1,1),
       col=c(1, 2))
dev.off();