Converts table of counts of resolved quartets on
n taxa to show only dominant one, with
maximum likelihood estimate of internal edge weight under the MSC.
quartetTableDominant(rqt, bigweights = "infinite")
a table, as produced by
bigweights="finite", when for a set of 4 taxa the quartet counts are (m,0,0) then
the edge weight is computed as if the relative frequency of the dominant topology were m/(m+1).
an (n choose 4)x(n+1) array with dominant quartet topology encoded by 1,1,-1,-1 in
taxon columns, with signs indicating cherries; the (n+1)th column
"weight" contains the maximum likelihood estimates,
under MSC on a 4-taxon tree, of the quartets' central edge lengths, in coalescent units
gtrees=read.tree(file=system.file("extdata","dataGeneTreeSample",package="MSCquartets")) tnames=taxonNames(gtrees) QT=quartetTable(gtrees,tnames[1:6]) RQT=quartetTableResolved(QT) RQT[1:6,] DQT=quartetTableDominant(RQT) DQT[1:6,]
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.