Warning: no labels for x. Assuming order of tree$tip.label I am using the phylANOVA() function in the R package phytools, as follows:
