using homsap mu

popsim-consortium · Aug 14, 2023 · a60b85b · a60b85b
1 parent 2296e68
commit a60b85b
Showing 1 changed file with 2 additions and 3 deletions.
diff --git a/stdpopsim/catalog/HomSap/dfes.py b/stdpopsim/catalog/HomSap/dfes.py
@@ -191,11 +191,10 @@ def _KyriazisDFE():
     # To get 0.63 deleterious mutations per 30Mb diploid coding sequence using 
     # mu from catalog, we need a deleterious proportion of:
     proportion_deleterious = 0.63 / (3e7 * _species.genome.mean_mutation_rate * 2)
-
 
     # Next, 0.03% of deleterious are lethals:
-    gamma_prop = 0.8139535 * 0.997
-    lethal_prop = 0.8139535 * 0.003
+    gamma_prop = proportion_deleterious * 0.997
+    lethal_prop = proportion_deleterious * 0.003
     neutral_prop = 1 - gamma_prop - lethal_prop
 
     # To get 0.63 deleterious mutations per 30Mb diploid coding sequence with