read_f2: Read blocked f2 estimates from disk

In uqrmaie1/admixtools: Inferring demographic history from genetic data

Read blocked f2 estimates from disk

Description

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.

Usage

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

## Not run: 
read_f2(f2_dir, pops = c('pop1', 'pop2', 'pop3'))

## End(Not run)

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