View source: R/TreeQuartetHypTest.R
quartetTreeTestInd  R Documentation 
Perform a tree hypothesis test for all quartet counts in an input table, as if the counts for different choices of 4 taxa are independent.
quartetTreeTestInd(
rqt,
model = "T3",
lambda = 0,
method = "MLest",
smallcounts = "approximate",
bootstraps = 10^4,
speciestree = NULL
)
rqt 
table of resolved quartet counts, as produced by 
model 

lambda 
power divergence statistic parameter (e.g., 0 for likelihood ratio statistic, 1 for Chisquared statistic) 
method 

smallcounts 

bootstraps 
number of samples for bootstrapping 
speciestree 
species tree, in Newick as text, to determine quartet for T1 test; required for 
This function assumes all quartets are resolved. The test performed and the arguments
are described more fully in quartetTreeTest
.
if model="T3"
, a copy of rqt
with a new column "p_T3"
appended with pvalues for each quartet;
if model="T1"
, a copy of rqt
with 2 columns appended: "p_T1"
with pvalues, and "qindex"
giving index of quartet consistent with specified species tree,
i.e., 1 if 1234 on species tree, 2 if 1324, 3 if 1423
MAR19MSCquartets
quartetTreeTest
, quartetTestPlot
, quartetStarTestInd
, quartetTableResolved
gtrees=read.tree(file=system.file("extdata","dataGeneTreeSample",package="MSCquartets"))
tnames=c("t1","t2","t3","t4","t5","t6")
QT=quartetTable(gtrees,tnames)
RQT=quartetTableResolved(QT)
stree="(((t5,t6),t4),((t1,t2),t3));"
pTable3=quartetTreeTestInd(RQT,"T3")
quartetTablePrint(pTable3[1:6,])
stree="((((t5,t6),t4),t7),((t8,t9),((t1,t2),t3)));"
pTable1=quartetTreeTestInd(RQT,"T1",speciestree=stree)
quartetTablePrint(pTable1[1:6,])
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.