Perform hypothesis tests for a species quartet star tree vs. any alternative for all quartet counts in an input table, as if the quartets are independent.
quartetStarTestInd(rqt)
the same table as the input rqt
with column "p_star"
appended, containing p-values for
judging fit to MSC on a star tree
Table of resolved quartet counts, as produced by quartetTableResolved
, or quartetTreeTestInd
This function assumes all quartets are resolved.
The test performed is described in quartetStarTest
.
quartetStarTest
, quartetTreeTest
, quartetTreeTestInd
,
quartetTableResolved
, quartetTestPlot
gtrees=read.tree(file=system.file("extdata","dataGeneTreeSample",package="MSCquartets"))
tnames=taxonNames(gtrees)
QT=quartetTable(gtrees,tnames[1:6])
RQT=quartetTableResolved(QT)
pTable=quartetStarTestInd(RQT)
quartetTablePrint(pTable[1:6,])
Run the code above in your browser using DataLab