Compute the F_ST statistic for pool-sequenced data of Kofler et al as used in PoPoolation2, for two ranges of SampleCountss.
The approach is called the "classical" or "conventional" estimator in PoPoolation2 [1], and follows Hartl and Clark [2].
[1] PoPoolation2: identifying differentiation between populations using sequencing of pooled DNA samples (Pool-Seq).
Kofler R, Pandey RV, Schlotterer C.
Bioinformatics, 2011, 27(24), 3435–3436. https://doi.org/10.1093/bioinformatics/btr589
[2] Principles of Population Genetics.
Hartl DL, Clark AG.
Sinauer, 2007.
Definition at line 70 of file fst_pool_kofler.hpp.
Compute the SNP-based Theta Pi values used above.
See there for details. The tuple returns Theta Pi for an individual position, which is simply the heterozygosity() at this position, for both samples p1
and p2
, as well as their combined (average frequency) heterozygosity, in that order.
Definition at line 174 of file fst_pool_kofler.hpp.