Table 1  Comparisons between theoretical predictions (from coalescent simulations) and forward-time Wright–Fisher simulations for Tajima’s D (Tajima 1989) in a sample of size 10, D10
Ub/sNsD10 theoryD10 simulations
0.20005000−3.3199−3.3378
0.100010000−3.3489−3.3569
0.020050000−3.3533−3.3322
0.0100100000−3.3571−3.4188
0.0020500000−3.3665−3.3024
0.00101000000−3.3717−3.3670
  • Here Ub = 10−4 and N = 107 while s is varied. Theoretical predictions are obtained by sampling 107 backward coalescent simulations. Forward-time simulation results are an average over 56 forward simulation runs, with 106 samples of π and S10 used to compute D10.