View source: R/StarQuartetHypTest.R
quartetStarTestInd | R Documentation |
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)
rqt |
Table of resolved quartet counts, as produced by |
This function assumes all quartets are resolved.
The test performed is described in quartetStarTest
.
the same table as the input rqt
with column "p_star"
appended, containing p-values for
judging fit to MSC on a star tree
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,])
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.