code for manipulating ballgown output in R - View it on GitHub