Orthogonal contrasts were made between each treatment and the 2008 baseline. In addition, all 2013 treatments
were compared to each other. Significance of individual contrasts were evaluated using false discovery rate corrections from the lsmeans function in the ‘lsmeans’ package (Lenth, 2013).