if (!inherits(phy, "phylo")) stop('object "phy" is not of class "phylo"')
if (!is.null(names(S))) {
if (all(names(S) %in% phy$tip.label)) S <- S[phy$tip.label]
- else warning('the names of argument "S" and the names of the tip labels
-did not match: the former were ignored in the analysis.')
+ else warning('the names of argument "S" and the tip labels
+did not match: the former were ignored.')
}
N <- length(S)
x <- branching.times(phy)