The phytools functions resolveNode and resolveAllNodes could be used.
They resolve either a single node in all possible ways; or all nodes in
all possible ways, returning a "multiPhylo" object. Note that in the
latter case the number of trees can be quite large. From the phytools
doc: "For resolveNode applied to a multifurcation with n descendants,
the number of resolved trees will be equal to the number of possible
rooted trees of n taxa. (For instance, three for a trifurcation, 15 for
a quadrifurcation, and so on.) For resolveAllNodes the number of fully
resolved trees will be equal to the product of numbers for resolveNode
applied to each multifurcation separately. (For instance, 45 for a tree
containing one trifurcation and one quadrifurcation.)" To get a random
set of resolved trees, rather than all possible trees (because this
number can be high), you could start with one multifurcating node in the
tree, resolve it in all ways using resolveNode, pick one of those
resolutions, then move to the next node & so on. I will try to post a
formal solution & then share it to the list.
All the best, Liam
Liam J. Revell, Associate Professor of Biology
University of Massachusetts Boston
On 6/21/2017 10:38 AM, Laura Jackson wrote:
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?
R-sig-phylo mailing list - Rfirstname.lastname@example.org
Searchable archive at http://email@example.com/