| phylANOVA | R Documentation |
Simulation based phylogenetic ANOVA following Garland et al. (1993), with post-hoc tests.
phylANOVA(tree, x, y, nsim=1000, posthoc=TRUE, p.adj="holm")
tree |
a phylogenetic tree in |
x |
a vector containing the groups. |
y |
a vector containing the response variable (continuously valued). |
nsim |
an integer specifying the number of simulations (including the observed data). |
posthoc |
a logical value indicating whether or not to conduct posthoc tests to compare the mean among groups. |
p.adj |
method to adjust P-values for the posthoc tests to account for multiple testing. Options same as |
This function performs the simulation-based phylogenetic ANOVA of Garland et al. (1993) and (optionally) conducts all posthoc comparisons of means among groups (also obtaining the P-values by phylogenetic simulation).
This function uses a little bit of code from both phy.anova in the geiger package and pairwise.t.test.
An object of class "phylANOVA" containing the following elements:
F |
F from observed data. |
Pf |
P-value for F from simulation. |
T |
matrix of t-values. |
Pt |
matrix of multiple test corrected P-values from posthoc t-tests. |
Liam Revell liam.revell@umb.edu
Garland, T., Jr., A. W. Dickerman, C. M. Janis, & J. A. Jones (1993) Phylogenetic analysis of covariance by computer simulation. Systematic Biology, 42, 265-292.
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. (2024) phytools 2.0: an updated R ecosystem for phylogenetic comparative methods (and other things). PeerJ, 12, e16505.
anova, pairwise.t.test
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.