data(bird.orders)
x <- rep(0, 5)
names(x) <- c("A", "B", "C", "D", "E")
anc1 <- evolve.phylo(bird.orders, x, 1)
anc2 <- evolve.phylo(bird.orders, x, 1)
cor(anc1$tip.character, anc2$tip.character)
Run the code above in your browser using DataCamp Workspace