pvalHist: Plot histogram of log p-values in table

In MSCquartets: Analyzing Gene Tree Quartets under the Multi-Species Coalescent

Plot histogram of log p-values in table

Description

Graphical exploration of extreme p-values from quartet hypothesis tests, to aid in choosing critical values for hypothesis tests. Log base 10 of p-values exceeding some minimum are plotted, to explore gaps in the tail of the distribution.

Usage

pvalHist(pTable, model, pmin = 0)

Arguments

pTable	a quartet table with p-values, such as from NANUQ, quartetTreeTestInd, or quartetStarTestInd
model	one of "T1", "T3", or "star", where pTable contains a column p_model of p-values
pmin	include only p-values above pmin in the histogram

Details

Since logarithms are plotted, p-values close to 0 will appear as negative numbers of large magnitude, putting the tail of the distribution to the left in the histogram.

When exploring possible critical values for the hypothesis tests in the NANUQ algorithm, use model= "T3" to choose values for alpha which distinguish treelikeness from hybridization, and model= "star" to choose values for beta which distinguish polytomies from resolved trees. In general, alpha should be chosen to be small and beta to be large so that most quartets are interpreted as resolved trees.

See Also

NANUQ, NANUQdist

Examples

data(pTableYeastRokas)
pvalHist(pTableYeastRokas,"T3")

MSCquartets documentation built on Oct. 31, 2024, 1:08 a.m.