data(woodmouse)
otree <- root(fastme.ols(dist.dna(woodmouse)), "No305", resolve.root=TRUE)
breps <- 500
trees <- boot.phylo(otree, woodmouse, B=breps, function(x)
root(fastme.ols(dist.dna(x)), "No305", resolve.root=TRUE),trees=TRUE)
combined.trees <- c(list(otree), trees$trees)
bin.multiPhylo(combined.trees)
Run the code above in your browser using DataLab