This functions calculates pairwise FST values according to Hudson's estimator
following Bathia. FST values are calculated for each pairwise combination of
the populations present in the data (defined by the nPops
parameter).
pairFST(nPops, Pop_pi, coverage)
a upper triangular matrix with the pairwise FST values between each population.
is an integer indicating how many populations are present in the data.
is a matrix of allele frequencies. Each row of that matrix should correspond to a different population and each column to a different site.
is a matrix containing depths of coverage. Each row of that matrix should correspond to a different population and each column to a different site.
This functions performs pairwise comparisons of multiple populations and thus returns multiple FST values, one for each comparison. However, this function computes FST for a single locus.