-## drop.tip.R (2010-11-24)
+## drop.tip.R (2011-11-21)
## Remove Tips in a Phylogenetic Tree
-## Copyright 2003-2010 Emmanuel Paradis
+## Copyright 2003-2011 Emmanuel Paradis
## This file is part of the R-package `ape'.
## See the file ../COPYING for licensing issues.
if (wbl) phy$edge.length <- phy$edge.length[keep]
TIPS <- phy$edge[, 2] <= Ntip
tip <- phy$edge[TIPS, 2]
- phy$tip.label <- phy$tip.label[sort(tip)] # <- added sort to avoid shuffling of tip labels (2010-07-21)
+ ## Fix by Ludovic Mallet and Mahendra Mariadassou (2011-11-21):
+ name <- vector("character", length(tip))
+ name[order(tip)] <- phy$tip.label[tip]
+ phy$tip.label <- name
+ ## End of fix
## keep the ordering so no need to reorder tip.label:
phy$edge[TIPS, 2] <- order(tip)
if (!is.null(phy$node.label))
{
if (!inherits(phy, "phylo"))
stop('object "phy" is not of class "phylo"')
+ if (!length(tip)) return(phy)
Ntip <- length(phy$tip.label)
## find the tips to drop: