def test_trivial_calc_multilocus_f_statistics(self): ml_stats = VCF.calc_multilocus_f_statistics(self.trivial_Hs_est_dict, self.trivial_Ht_est_dict) self.assertEqual({('pop1', 'pop2'): {'G_prime_st_est': 0.0, 'Gst_est.stdev': 0.0, 'G_double_prime_st_est.stdev': 0.0, 'G_double_prime_st_est': 0.0, 'Gst_est': 0.0, 'D_est.stdev': 0.0, 'G_prime_st_est.stdev': 0.0, 'D_est': 0.0}}, ml_stats)
def test_trivial_2_calc_multilocus_f_statistics(self): ml_stats = VCF.calc_multilocus_f_statistics(self.trivial_2_Hs_est_dict, self.trivial_2_Ht_est_dict) self.assertEqual({('pop1', 'pop2'): {'G_prime_st_est': 0.8333333333333334, 'Gst_est.stdev': 0.09999999999999998, 'G_double_prime_st_est.stdev': 0.060586734693877375, 'G_double_prime_st_est': 0.8888888888888888, 'Gst_est': 0.5, 'D_est.stdev': 0.08928571428571419, 'G_prime_st_est.stdev': 0.07857142857142851, 'D_est': 0.6488650338184054}}, ml_stats)
def test_calc_multilocus_f_statistics(self): ml_stats = VCF.calc_multilocus_f_statistics(self.Hs_est_dict, self.Ht_est_dict) self.assertEqual({('pop1', 'pop2'): {'Gst_est': 0.30312672938572266, 'Gst_est.stdev': 1.3360742705570265, 'G_double_prime_st_est': 0.5003506191636901, 'G_double_prime_st_est.stdev': 2.394045897494315, 'G_prime_st_est': 0.34889353651117816, 'G_prime_st_est.stdev': 1.6091544573608876, 'D_est': 0.0015629851350084287, 'D_est.stdev': 0.25110803099568757}}, ml_stats)