powered by
Perform a hypothesis test for all quartet counts in an input table, as if the counts for different choices of 4 taxa are independent.
quartetCutTestInd( rqt, lambda = 0, method = "MLest", smallcounts = "approximate", bootstraps = 10^4 )
a copy of rqt with two columns appended: "p_cut" with p-values and "cutindex"
rqt
"p_cut"
"cutindex"
giving index 1,2, or 3 of ML estimate of quartet tree of blobs (1 if 12|34, 2 if 13|24, 3 if 14|23) under Cut model.
table of resolved quartet counts, as produced by quartetTableResolved, or quartetStarTestInd
quartetTableResolved
quartetStarTestInd
power divergence statistic parameter (e.g., 0 for likelihood ratio statistic, 1 for Chi-squared statistic)
"MLest", "conservative", or "bootstrap"; see quartetCutTest for explanation
"MLest"
"conservative"
"bootstrap"
quartetCutTest
"bootstrap" or "approximate", method of obtaining p-value when some counts are small, so the chosen method is inappropriate
"approximate"
method
number of samples for bootstrapping
This function assumes all quartets are resolved. The test performed and the arguments are described more fully in quartetCutTest.
ABMR24MSCquartets
quartetCutTest, quartetTestPlot, quartetStarTestInd, quartetTableResolved
quartetTestPlot
data(pTableYeastRokas) pT=pTableYeastRokas[1:10,1:11] pTable=quartetCutTestInd(pT)
Run the code above in your browser using DataLab