positive selection for HomSap

popsim-consortium · Jan 19, 2025 · f3ac566 · f3ac566
1 parent 46ae171
commit f3ac566
Showing 1 changed file with 8 additions and 4 deletions.
diff --git a/stdpopsim/catalog/HomSap/dfes.py b/stdpopsim/catalog/HomSap/dfes.py
@@ -234,6 +234,10 @@ def _ZhenDFE():
         distribution_type="g",  # gamma distribution
         distribution_args=[gamma_mean, gamma_shape],
     )
+    # p. 2 in supplement says that the total sequence length of synonymous sites LS
+    # related to the total sequence length LNS by LNS = 2.31 * LS
+    prop_synonymous = 1 / (1 + 2.31)
+
     sel_coeff = 10 ** (-3.949)
     prop_beneficial = 1.55e-2
     positive = stdpopsim.MutationType(
@@ -248,10 +252,10 @@ def _ZhenDFE():
         long_description=long_description,
         mutation_types=[neutral, negative, positive],
         proportions=[
-            0.3,
-            0.7 * (1 - prop_beneficial),
-            0.7 * prop_beneficial,
-        ],  # [0.3 and 0.7 are from Gamma_H17
+            prop_synonymous,
+            (1 - prop_synonymous) * (1 - prop_beneficial),
+            (1 - prop_synonymous) * prop_beneficial,
+        ],
         citations=citations,
     )