The simulation in phytools assumes a Brownian motion model of trait evolution. Post-hoc tests for pairwise group comparisons were also conducted in phytools based on distributions of pairwise t-statistics derived from the simulations, with a correction for multiple tests.