compute_F4DfromF2samples: compute_F4DfromF2samples

In poolfstat: Computing f-Statistics and Building Admixture Graphs Based on Allele Count or Pool-Seq Read Count Data

compute_F4DfromF2samples

Description

Compute all F4 and Dstat from F2 values obtained from each block-jackknife bloc

Usage

.compute_F4DfromF2samples(blockF2, blockDenom, npops, verbose)

Arguments

blockF2	Numeric Matrix with nF2=(npop*(npop-1))/2 rows and nblocks columns matrix containing pairwise-pop F2 estimates for each block-jackknife sample (l.o.o.)
blockDenom	Numeric Matrix with nF4=(npops*(npops-1)/2)*((npops-2)*(npops-3)/2)/2 rows and nblocks containing the estimates of the denominator of Dstat (see compute_blockDdenom) for each block-jackknife sample (l.o.o.)
npops	Integer giving the number of populations
verbose	Logical (if TRUE progression bar is printed on the terminal)

Details

Compute F4 and D estimates and their s.e. based on block-jackknife estimates of all F2 (and heterozygosities)

Value

Return a matrix with nF4=(npops*(npops-1)/2)*((npops-2)*(npops-3)/2)/2 rows and four columns corresponding to the mean and the s.e. of F4 and the mean and s.e. of Dstat

Examples

#

poolfstat documentation built on Sept. 8, 2023, 5:49 p.m.