def test_estimate_nm(self):
"""Test Nm estimation."""
ctrl = GenePopController()
mean_sample_size, mean_priv_alleles, mig10, mig25, mig50, mig_corrected =\
ctrl.estimate_nm("PopGen" + os.sep + "big.gen")
assert (mean_sample_size, mean_priv_alleles, mig10, mig25, mig50, mig_corrected) == \
(28.0, 0.016129, 52.5578, 15.3006, 8.94583, 13.6612)
def test_estimate_nm(self):
"""Test Nm estimation."""
ctrl = GenePopController()
mean_sample_size, mean_priv_alleles, mig10, mig25, mig50, mig_corrected =\
ctrl.estimate_nm("PopGen" + os.sep + "big.gen")
assert (mean_sample_size, mean_priv_alleles, mig10, mig25, mig50, mig_corrected) == \
(28.0, 0.016129, 52.5578, 15.3006, 8.94583, 13.6612)
def test_estimate_nm(self):
"""Test Nm estimation."""
ctrl = GenePopController()
path = os.path.join("PopGen", "big.gen")
mean_sample_size, mean_priv_alleles, mig10, mig25, mig50, mig_corrected = ctrl.estimate_nm(path)
self.assertAlmostEqual(mean_sample_size, 28.0)
self.assertAlmostEqual(mean_priv_alleles, 0.016129)
self.assertAlmostEqual(mig10, 52.5578)
self.assertAlmostEqual(mig25, 15.3006)
self.assertAlmostEqual(mig50, 8.94583)
self.assertAlmostEqual(mig_corrected, 13.6612)
