Merge pull request #82 from andrewkern/sheehan_song_fix

fixed a bug with the times in sheehan and song model
popsim-consortium · May 28, 2019 · 59a37c7 · 59a37c7
2 parents 2542d14 + 1337fc7
commit 59a37c7
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/stdpopsim/drosophila_melanogaster.py b/stdpopsim/drosophila_melanogaster.py
@@ -97,7 +97,7 @@ def __init__(self):
         # Times are provided in 4N_ref generations, so we convert into generations.
         # generation_time = 10 / year
         t_1 = t_1_coal * 4 * N_ref
-        t_2 = t_2_coal * 4 * N_ref
+        t_2 = (t_1_coal + t_2_coal) * 4 * N_ref
         # Single population in this model
         self.population_configurations = [
             msprime.PopulationConfiguration(initial_size=N_R),