I share Ben's thought. I gave it a try with simulated trees and an
analysis similar to Stacey's:
R tr - rmtree(1e3, 100)
R system.time(lapply(tr, drop.tip, paste(t, 1:5, sep = )))
utilisateur système écoulé
17.278 0.050 17.359
R system.time(for (i in 1:1000) tr[[i]] -
Hmm. I should try solutions before I post them. You need to make sure
that each tree in the list is of class phylo too. This works:
phylist-read.tree(text=A:0.21,B:0.21):0.28,C:0.49):0.13,D:0.62):0.38,E:1.00);A:0.21,B:0.21):0.28,C:0.49):0.13,D:0.62):0.38,E:1.00);)
newnames-c(G,H,I,J,K)
Try
class(newlist) - multiPhylo
Then use write.tree.
Cheers,
Simon.
On Sun, 2008-08-31 at 20:39 -0400, [EMAIL PROTECTED] wrote:
Hi all,
I have hit an obstacle and I hope someone will know a quick fix. I
want to read
a list of trees, do something to those trees and then write them to a