phylANOVA(tree, x, y, nsim=1000, posthoc=TRUE, p.adj="holm")
"phylo"
format.p.adjust
.phy.anova
in the "geiger" package as well as pairwise.t.test
.
Harmon, L. J., J. T. Weir, C. D. Brock, R. E. Glor, W. Challenger. (2008) GEIGER: investigating evolutionary radiations. Bioinformatics, 24, 129-131.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.
anova
, pairwise.t.test