library(ape)
data(bird.orders)
tree <- bird.orders
mat.list <- RandomMatrix(5, length(tree$tip.label))
names(mat.list) <- tree$tip.label
sample.sizes <- runif(length(tree$tip.label), 15, 20)
phylo.state <- PhyloW(tree, mat.list, sample.sizes)
phylo.comparisons <- PhyloCompare(tree, phylo.state)
# plotting results on a phylogeny:
if (FALSE) {
library(phytools)
plotBranchbyTrait(tree, phylo.comparisons[[2]])
}
Run the code above in your browser using DataLab