compute_F3fromF2samples {poolfstat}

compute_F3fromF2samples

Description

Compute all F3 from F2 values obtained from each block-jackknife bloc

Usage

.compute_F3fromF2samples(blockF2, blockHet, 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.)

blockHet

Numeric Matrix with npop rows and nblocks columns containing all within pop heterozygosity estimates 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 F3 and F3star estimates and their s.e. based on block-jackknife estimates of all F2 (and heterozygosities)

Value

Return a matrix with nF3=npops*(npops-1)*(npops-2)/2 rows and four columns corresponding to the mean and the s.e. of F3 and the mean and s.e. of F3star

Examples

#