powered by
mast computes the maximum agreement subtree (MAST).
mast
mast(x, y, tree = TRUE, rooted = TRUE)
a tree, i.e. an object of class phylo.
phylo
a logical, if TRUE returns a tree other wise the tip labels of the the maximum agreement subtree.
logical if TRUE treats trees as rooted otherwise unrooted.
mast returns a vector of the tip labels in the MAST.
The code is derived from the code example in Valiente (2009). The version for the unrooted trees is much slower.
G. Valiente (2009). Combinatorial Pattern Matching Algorithms in Computational Biology using Perl and R. Taylor & Francis/CRC Press
SPR.dist
# NOT RUN { tree1 <- rtree(100) tree2 <- rSPR(tree1, 5) tips <- mast(tree1, tree2) # }
Run the code above in your browser using DataLab