Compute the Weighted Quartet Distance between taxa of YR19;textualMSCquartets from a table specifying a collection of quartets on
n
taxa and the quartets' internal branch lengths.
quartetWeightedDist(dqt)
a pairwise distance matrix on n
taxa
an (n
choose 4) x (n+1
) matrix of the form output by quartetTableDominant
YR19MSCquartets
quartetTableDominant
,
WQDSAdjustLengths
,
WQDS
,
WQDC
,
WQDCrecursive
,
quartetWeightedDist
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)
stree=WQDSAdjustLengths(tree)
write.tree(stree)
Run the code above in your browser using DataLab