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
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
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:
>
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/
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.