estimateTheta: Estimates Fst and quantiles of Fst

In Ahhgust/MMDIT: Mitochondrial mixture project

Description

This function computes Buckleton's estimator of Fst (theta) and it computes an upper-bound. It's a friendly wrapper for fst_buckleton

Usage

estimateTheta(
  alleles,
  populations,
  quantile = 0.95,
  nJack = 0,
  approximate = FALSE
)

Arguments

alleles	strings ; alleles (haplotypes); 1 row per chromosome sampled
populations	strings ; same length as alleles (parallel array); population label associated with said haplotype
quantile	; [0,1] ; what quantile in the Fst distribution would you like?
nJack	; non-negative integer ; number of leave-one-out jackknifes taken
approximate	; boolean ; whether or not allele frequencies (approximate=TRUE) or site hetero/homozygosities (FALSE) should be used

Ahhgust/MMDIT documentation built on Jan. 27, 2021, 11:48 a.m.