compute_F4DfromF2samples {poolfstat}R Documentation

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

#

[Package poolfstat version 2.2.0 Index]