calculate_fst: Calculate FST

In GenomeAdmixR: Simulate Admixture of Genomes

Calculate FST

Description

The FST value between two populations is calculated, given a number of markers. Markers are superimposed upon the (known) ancestry along the chromosome for all sampled individuals. Markers can be chosen to be regularly spaced, or randomly distributed.

Usage

calculate_fst(
  pop1,
  pop2,
  sampled_individuals = 10,
  number_of_markers = 100,
  random_markers = FALSE
)

Arguments

pop1	Population object
pop2	Population object
sampled_individuals	Number of individuals to base the FST upon. Individuals are randomly drawn from each population, a lower number speeds up calculations.
number_of_markers	Number of markers along the chromosome used to calculate FST metrics.
random_markers	If TRUE, markers are randomly spaced along the chromosome, if FALSE, markers are equidistantly spaced along the chromosome.

Details

Uses the function wc from the package hierfstat to calculate the FST. The function wc computes the Weir and Cockerham F statistic.

Value

FST value

Examples

two_populations <- simulate_admixture(
                        module = ancestry_module(),
                       migration = migration_settings(migration_rate = 0.01,
                                           population_size = c(100, 100)))

FST <- calculate_fst(pop1 = two_populations$population_1,
                     pop2 = two_populations$population_2,
                     sampled_individuals = 10,
                     number_of_markers = 100,
                     random_markers = TRUE)

GenomeAdmixR documentation built on March 18, 2022, 5:40 p.m.