def test_hardcoded_fixations(self): game = games.two_times_two(a=3.0, b=0.5, c=0.5, d=2.0) intensity_of_selection = 1.0 population_size = 5 # number_of_samples=100000 number_of_samples = 20000 index_of_the_incumbent = 1 index_of_the_mutant = 0 # numerical_value = numerical.fixation_probability_strategy_a(game, intensity_of_selection, population_size) numerical_value = numerical.fixation_probability( mutant_index=0, resident_index=1, intensity_of_selection=intensity_of_selection, population_size=population_size, game_matrix=game, mapping='EXP') print(numerical_value) mp = MoranProcess(population_size, intensity_of_selection, game_matrix=game, fitness_mapping='exp', mutation_probability=0.1) simulated_value = mp.simulate_fixation_probability( index_of_the_incumbent=index_of_the_incumbent, index_of_the_mutant=index_of_the_mutant, number_of_samples=number_of_samples, seed=None) print(simulated_value) self.assertAlmostEqual(numerical_value, simulated_value, places=2, msg="", delta=None)
def test_fixation_of_neutral_mutant(self): number_of_strategies_value = 2 number_of_samples_=10000 for _ in xrange(0,5): pop_size =np.random.randint(2,11) mp = MoranProcess(population_size=pop_size, intensity_of_selection=0.0, game_matrix=np.random.rand(number_of_strategies_value,number_of_strategies_value), number_of_strategies=number_of_strategies_value, fitness_mapping='exp', mutation_probability=0.1) fix = mp.simulate_fixation_probability(0, 1, number_of_samples=number_of_samples_, seed=None) np.testing.assert_allclose(fix,1.0/pop_size, rtol=0.01, atol=0.01, err_msg="Paila", verbose=True)
def test_fixation_of_neutral_mutant(self): number_of_strategies_value = 2 number_of_samples_ = 10000 for _ in range(0, 5): pop_size = np.random.randint(2, 11) mp = MoranProcess(population_size=pop_size, intensity_of_selection=0.0, game_matrix=np.random.rand(number_of_strategies_value, number_of_strategies_value), number_of_strategies=number_of_strategies_value, fitness_mapping='exp', mutation_probability=0.1) fix = mp.simulate_fixation_probability(0, 1, number_of_samples=number_of_samples_, seed=None) np.testing.assert_allclose(fix, 1.0 / pop_size, rtol=0.01, atol=0.01, err_msg="Paila", verbose=True)
def test_hardcoded_fixations(self): game= games.two_times_two(a=3.0, b=0.5, c=0.5, d=2.0) intensity_of_selection=1.0 population_size=5 #number_of_samples=100000 number_of_samples=20000 index_of_the_incumbent=1 index_of_the_mutant=0 #numerical_value = numerical.fixation_probability_strategy_a(game, intensity_of_selection, population_size) numerical_value = numerical.fixation_probability(mutant_index=0, resident_index=1, intensity_of_selection=intensity_of_selection, population_size=population_size, game_matrix=game, mapping='EXP') print numerical_value mp = MoranProcess(population_size, intensity_of_selection, game_matrix=game, fitness_mapping='exp', mutation_probability=0.1) simulated_value = mp.simulate_fixation_probability(index_of_the_incumbent=index_of_the_incumbent, index_of_the_mutant=index_of_the_mutant, number_of_samples=number_of_samples, seed=None) print simulated_value self.assertAlmostEqual(numerical_value, simulated_value, places=2, msg="", delta=None)