qpwave_pairs: Compute all pairwise qpwave p-values

In uqrmaie1/admixtools: Inferring demographic history from genetic data

Compute all pairwise qpwave p-values

Description

For all pairs of left populations qpwave rank 0 Chi-squared statistics and p-values will be computed. This tests for each pair of left populations whether it forms a clade with respect to the right populations.

Usage

qpwave_pairs(f2_data, left, right)

Arguments

f2_data	A 3d array of blocked f2 statistics, output of f2_from_precomp. Alternatively, a directory with f2 statistics.
left	Left populations
right	Right populations

Examples

## Not run: 
left = pops[5:7]
right = pops[1:4]
f2_blocks = f2_from_precomp('/my/f2/dir/', pops = left, pops2 = right, afprod = TRUE)
qpwave_pairs(f2_blocks, left, right)
# If f2-stats are too big to load them into memory,
# the following will read the data for each pair from disk:
qpwave_pairs('/my/f2/dir/', left, right)

## End(Not run)

uqrmaie1/admixtools documentation built on Nov. 3, 2024, 12:56 a.m.