Adjust edge lengths on tree built from Weighted Quartet distance to estimate metric tree
Modify edge lengths of a tree built from a distance table produced by quartetWeightedDist,
to remove scaling factors related to the size of the split associated to the edge.
WQDSAdjustLengths(tree)
tree |
an unrooted metric tree, of type phylo |
As explained by Yourdkhani and Rhodes (2020), a metric tree produced from the weighted quartet distance has edge lengths inflated by a factor dependent on the associated split size. Removing these factors yields a consistent estimate of the metric species tree displaying the weighted quartets, if such a tree exists.
This function should not be used on trees output from WQDS, WQDC,
or WQDCrecursive, as
their edges are already adjusted. It can be used on trees built from the distance
computed by quartetWeightedDist.
an unrooted metric tree, of type phylo
Yourdkhani S, Rhodes JA (2020). “Inferring metric trees from weighted quartets via an intertaxon distance.” Bul. Math. Biol., 82(97). doi: 10.1007/s11538-020-00773-4, https://doi.org/10.1007/s11538-020-00773-4.
gtrees=read.tree(file=system.file("extdata","dataGeneTreeSample",package="MSCquartets"))
tnames=taxonNames(gtrees)
QT=quartetTable(gtrees,tnames[1:6])
RQT=quartetTableResolved(QT)
DQT=quartetTableDominant(RQT,bigweights="finite")
D=quartetWeightedDist(DQT)
tree=NJ(D)
write.tree(tree)
plot(tree)
stree=WQDSAdjustLengths(tree)
write.tree(stree)
plot(stree)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.