Re: [R-sig-phylo] Problem date consensus phylogeny

2019-01-23 Thread Emmanuel Paradis
Hi Alina, Did you try multi2di() to remove polytomies? Best, Emmanuel Le 23/01/2019 à 17:41, Alina van dijk a écrit : Hi everyone, My name is Alina, I'm attending master degree in ecology and I have a little problem with my consensus phylogeny. I used 1000 phylogeny of birdtree to construct

Re: [R-sig-phylo] Problem date consensus phylogeny

2019-01-23 Thread Joseph Brown
I typically use a non-R tool like treeannotator (from the BEAST package) to do this type of thing, especially when involving large numbers of big trees. HTH. Joseph. On Wed, Jan 23, 2019 at 4:42 PM Alina van dijk wrote: > Hi everyone, > > My name is Alina, I'm attending master degree in

[R-sig-phylo] Problem date consensus phylogeny

2019-01-23 Thread Alina van dijk
Hi everyone, My name is Alina, I'm attending master degree in ecology and I have a little problem with my consensus phylogeny. I used 1000 phylogeny of birdtree to construct the consensus tree but now I need to insert the dates in this consensus tree. To accomplish this idea, I tried different

Re: [R-sig-phylo] exponential time with write.tree

2019-01-23 Thread Jérémie Bardin
Thanks a lot Emmanuel. It works perfectly. Amitiés, Jérémie. Jérémie Bardin, Dr. CR2P - Centre de Recherche en Paléontologie - ParisSorbonne Université - MNHN - CNRS Site Jussieu, Tour 46-56, 5°et. 4 place Jussieu, 75252 Paris Cedex 05 tel. +331.44.27.51.77. jeremie.bar...@upmc.fr /

Re: [R-sig-phylo] exponential time with write.tree

2019-01-23 Thread Emmanuel Paradis
Hi Jérémie, You're correct: the current version of write.tree() scales with n^2. This seems to be related to calling the generic '[[' operator. You can modify the code with the usual fix(write.tree), then insert these two lines: phy <- .uncompressTipLabel(phy) class(phy) <- NULL

[R-sig-phylo] exponential time with write.tree

2019-01-23 Thread Jérémie Bardin
Hi all, I tried to export 2 millions trees to a .tre file with write.tree. This never ends, thus i tried to get why and realized that the time to write the .tre depends exponentially on the number of trees. trees<-allTrees(9,rooted=T)to<-seq(1,200,1)durations<-NAfor(i in