## ace.R (2009-06-19)
## Ancestral Character Estimation
## Copyright 2009 Johan Nylander and Emmanuel Paradis
## This file is part of the R-package `ape'.
## See the file ../COPYING for licensing issues.
is.monophyletic <-
function(phy, tips, reroot = !is.rooted(phy), plot = FALSE, ...)
{
if (!inherits(phy, "phylo"))
stop("object 'phy' is not of class 'phylo'")
if (length(tips) == 1) return(TRUE)
n <- length(phy$tip.label)
if (length(tips) == n) return(TRUE)
ROOT <- n + 1
if (is.numeric(tips)) {
if (any(tips > n))
stop("incorrect tip#: should not be greater than the number of tips")
tips <- sort(tips)
}
if (is.character(tips))
tips <- which(phy$tip.label %in% tips)
if (reroot) {
outgrp <- phy$tip.label[-tips][1]
phy <- root(phy, outgroup = outgrp, resolve.root = TRUE)
rerooted <- TRUE
} else rerooted <- FALSE
phy <- reorder(phy)
seq.nod <- .Call("seq_root2tip", phy$edge, n, phy$Nnode, PACKAGE = "ape")
sn <- seq.nod[tips]
newroot <- ROOT
i <- 2
repeat {
x <- unique(unlist(lapply(sn, "[", i)))
if (length(x) != 1) break
newroot <- x
i <- i + 1
}
desc <- which(unlist(lapply(seq.nod, function(x) any(x %in% newroot))))
if (plot) {
zoom(phy, tips, subtree = FALSE, ...)
if (rerooted)
mtext("Input tree arbitrarily rerooted", side = 1, cex = 0.9)
}
## assuming that both vectors are sorted:
identical(tips, desc)
} # end of is.monophyletic
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.