Re: [R-sig-phylo] Randomly Resolved Polytomies

2017-06-21 Thread Liam J. Revell
Hi Laura. I think multi2di also performs random resolutions, so if you wanted to generate, say, 100 random resolutions could you not do: trees<-replicate(100,multi2di(tree),simplify=FALSE) class(trees)<-"multiPhylo" All the best, Liam Liam J. Revell, Associate Professor of Biology

Re: [R-sig-phylo] Randomly Resolved Polytomies

2017-06-21 Thread Laura Jackson
Hi Liam, I was trying to get this to work on my tree, but it seems that it doesn't work for larger phylogenies (~1000 tips)? I get the following error message,..."Error in allTrees(n, TRUE, dd) : That would generate 7905853580625 trees, and take up more than 7905853581 MB of memory! " I have

Re: [R-sig-phylo] Randomly Resolved Polytomies

2017-06-21 Thread Laura Jackson
Liam, Thanks so much for your reply and for the blog post, I was able to get it working. -Laura On Wed, Jun 21, 2017 at 10:35 AM, Liam J. Revell wrote: > Dear Laura. > > I also just posted a possible solution to my blog here: >

Re: [R-sig-phylo] Randomly Resolved Polytomies

2017-06-21 Thread Liam J. Revell
Dear Laura. I also just posted a possible solution to my blog here: http://blog.phytools.org/2017/06/generating-set-of-random-resolutions-of.html. All the best, Liam Liam J. Revell, Associate Professor of Biology University of Massachusetts Boston web: http://faculty.umb.edu/liam.revell/

[R-sig-phylo] Randomly Resolved Polytomies

2017-06-21 Thread Laura Jackson
Hi all, I am using the ape package to randomly resolve polytomies using 'multi2di' and wondering if there is a way to use this function to get a single output tree file that contains multiple different randomly resolved trees using some number of resamplings? Thanks, -Laura -- *Laura M.