Table 2 MLEs of various DFEs
Data setDFEParameter MLEsLog-likelihoodAIC|ΔAIC|a
LuCampGammaα = 0.215, β = 562.1−3334.76673.413
Neu+gammapneu = 0.164, α = 0.338, β = 367.7−3327.26660.40
Neu+exp+letpneu = 0.304, pexp = 0.613, λ = 66.56−3337.86681.621.2
Discretebp1 = 0.309, p2 = 0.024, p3 = 0.247, p4 = 0.372−3334.16676.215.8
1kG EURGammaα = 0.186, β = 875.0−1450.52905.00
Neu+gammapneu = 0.031, α = 0.199, β = 820.6−1450.82907.62.6
Neu+exp+letpneu = 0.312, pexp = 0.509, λ = 41.48−1472.02950.045
Discretebp1 = 0.286, p2 = 0.099, p3 = 0.222, p4 = 0.305−1453.42914.89.8
ESP EURGammaα = 0.169, β = 1327.4−3012.66029.22.6
Neu+gammapneu = 0.092, α = 0.207, β = 1082.3−3010.36026.60
Neu+exp+letpneu = 0.341, pexp = 0.504, λ = 63.90−3071.66149.2122.6
Discretebp1 = 0.334, p2 = 0.041, p3 = 0.201, p4 = 0.306−3029.56067.040.4
  • These results are reported assuming LNS/LS = 2.31 and μ = 1.5×10−8. See Table S4 in File S1 for additional information. The shape and scale parameters of the gamma distribution are denoted with α and β, respectively, and the rate parameter of the exponential distribution is denoted with λ. Neu, neutral; exp, exponential; let, lethal; 1kG, 1000 Genomes.

  • a Change in AIC relative to the model with the lowest AIC.

  • b In terms of |s|, these parameters correspond to the ranges of |s| corresponding to: [0, 10−5), [10−5, 10−4), [10−4, 10−3), and [10−3, 10−2), respectively. The mass in the [10−2, 1] range is the complement of the total mass of the four aforementioned categories.