This function computes FST according to Hudson's estimator following Bathia. Note that the frequencies for the two populations should be entered as separate inputs.
getFst(freq1, freq2, ss1, ss2)a numeric value which is the FST between the two populations.
is a numeric vector with the allele frequencies for population
Each entry of the vector should be a different site.
is a numeric vector with the allele frequencies for population 2. Each entry of the vector should be a different site.
vector with the sample size for population 1. Each entry of the vector should contain the number of reads for a different site.
vector with the sample size for population 2. Each entry of the vector should contain the number of reads for a different site.
Note that this functions computes a single FST value between two populations and does not perform pairwise comparisons of multiple populations.