def test_calc_fstats_normal_allele_counts(self):
f_statistics = VCF.calc_fstats(self.normal_allele_counts)
self.assertEqual(
{'G_prime_st_est': 0.9140159767610748,
'D_est': 0.7581699346405228,
'G_double_prime_st_est': 0.9477124183006534,
'Gst_est': 0.6444444444444445,
'Hs_est': 0.1729729729729729,
'Ht_est': 0.48648648648648646},
f_statistics[('pachi', 'outgroups')])
def test_calc_fstats_trivial_allele_counts(self):
f_statistics = VCF.calc_fstats(self.trivial_allele_counts)
self.assertEqual(
{('pop2', 'pop1'):
{'G_prime_st_est': 0.6803049722304385,
'D_est': 0.517460317460318,
'G_double_prime_st_est': 0.7609710550887027,
'Gst_est': 0.33747412008281613,
'Hs_est': 0.3368421052631577,
'Ht_est': 0.508421052631579}},
f_statistics)
