Skip to content

Commit

Permalink
different neutral prop
Browse files Browse the repository at this point in the history
  • Loading branch information
chriscrsmith committed Aug 15, 2023
1 parent b643b08 commit ff1c6a7
Showing 1 changed file with 4 additions and 11 deletions.
15 changes: 4 additions & 11 deletions stdpopsim/catalog/HomSap/dfes.py
Original file line number Diff line number Diff line change
Expand Up @@ -188,17 +188,10 @@ def _KyriazisDFE():
distribution_args=[s],
dominance_coeff=h,
)
# 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 = 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
# mu=1.29e-8 (current mu from catalog), we need a proportion of 0.8139535:
proportion_deleterious = 0.7 # ratio 2.31:1
lethal_prop = proportion_deleterious * 0.003 # 0.03% lethals
gamma_prop = proportion_deleterious - lethal_prop
neutral_prop = 1 - proportion_deleterious
return stdpopsim.DFE(
id=id,
description=description,
Expand Down

0 comments on commit ff1c6a7

Please sign in to comment.