Using plyr to run anova on each phyla in each horizon