Read blocked f2 estimates from disk

This function reads blocked f2 estimates (or allele frequency products) which were writtend to disk by write_f2 and returns them as a 3d array.

read_f2(
  f2_dir,
  pops = NULL,
  pops2 = NULL,
  type = "f2",
  counts = FALSE,
  remove_na = TRUE,
  verbose = FALSE
)

Arguments

f2_dir

Directory from which to read files

pops

Populations for which f2 statistics should be read. Defaults to all populations, which may require a lot of memory.

pops2

Specify this if you only want to read a subset of all population pairs. The resulting array will differ on 1st and 2nd dimension and will not work with all functions.

counts

Return allele counts instead of f2 estimates

remove_na

Remove blocks with missing values

verbose

Print progress updates

afprod

Return allele frequency products instead of f2 estimates

Value

A 3d array of block jackknife estimates

See also

write_f2

Examples

if (FALSE) {
read_f2(f2_dir, pops = c('pop1', 'pop2', 'pop3'))
}