Table 1 Results from asymptoticMK for simulation runs conducted with SLiM 2
ModelαtrueαasymptoticαoriginalΔasymptoticΔoriginalρexp
Baseline0.329 ± 0.0150.307 ± 0.0580.164 ± 0.0350.0450.1650.75
L = 1080.327 ± 0.0080.301 ± 0.0130.174 ± 0.0120.0250.1521.00
L = 1060.321 ± 0.0670.246 ± 0.1340.142 ± 0.1410.1200.1910.15
μ = 10−80.306 ± 0.0050.287 ± 0.0160.173 ± 0.0090.0190.1321.00
μ = 10−100.317 ± 0.0570.288 ± 0.1690.145 ± 0.0740.1340.1730.05
rb = 0.00100.493 ± 0.0180.481 ± 0.0450.378 ± 0.0250.0410.1140.70
rb = 0.00010.091 ± 0.0140.115 ± 0.080−0.103 ± 0.0530.0710.1940.55
sb = 0.200.477 ± 0.0160.451 ± 0.0320.366 ± 0.0250.0290.1110.70
sb = 0.020.096 ± 0.0110.090 ± 0.068−0.119 ± 0.0470.0570.2150.50
sd = −0.2000.424 ± 0.0240.422 ± 0.0420.289 ± 0.0360.0320.1350.60
sd = −0.0020.233 ± 0.0110.234 ± 0.0570.104 ± 0.0390.0450.1290.50
T = 2 × 1060.324 ± 0.0060.302 ± 0.0140.173 ± 0.0120.0220.1511.00
T = 2 × 1040.345 ± 0.0630.369 ± 0.1830.225 ± 0.1130.1260.1200.05
  • The first row shows the averaged results from 20 replicate runs of the baseline SLiM model supplied on GitHub (see text). These runs used parameter values of mutation rate μ = 10−9 per base position per generation, chromosome length L = 107 base positions, beneficial mutation rate rb = 0.0005, beneficial mutation selection coefficient sb = 0.1, deleterious mutation selection coefficient sd = −0.02, and time after burn-in T = 2 × 105 generations. Each subsequent row shows the results from 20 replicate runs using the nonbaseline parameter value shown. αtrue specifies the true value of α averaged across the 20 replicates in each row; αasymptotic and αoriginal specify the asymptoticMK estimate and the estimate from the original test, respectively. SDs across the 20 replicates of each row are shown as ± values. Δasymptotic = |αasymptoticαtrue| and Δoriginal = |αoriginalαtrue| specify the mean absolute errors between true α values and the estimates from asymptoticMK and the original test, respectively, in each run, averaged over the 20 replicates. ρexp specifies the fraction of runs in which the exponential fit was chosen.