Hi Liam,
It is fixed and on the SVN:
https://svn.mpl.ird.fr/ape/dev/ape/R/drop.tip.R
Do not hesitate to test it because this function has been written some
years ago and has been patched quite a few times.
Thanks for the accurate report.
Emmanuel
On Sat, 20 Nov 2010 15:57:05 -0500
"Liam J
Hi list.
The following seems to reveal a bug in drop.tip{ape}:
require(geiger)
tr<-read.tree(text="(((A:1.0,B:1.0):1.0,C:2.0):1.0,(D:1.0,E:1.0):2.0);")
plot(drop.tip(tr,c("A","B"),trim.internal=FALSE))
# or, equivalently
write.tree(drop.tip(tr,c("A","B"),trim.internal=FALSE))
# which yields: [1]
