TABLE 2

Simulation results: simulated data conform to the analysis model

Mean parameter estimates (SD)
Simulated βSimulated N2E(s)N1 assumedβN2E(s)
a. Constant population
0.25200.19 (0.054)5.0 (1.45)
1000.21 (0.057)4.6 (1.14)
1200.84 (0.18)5.6 (1.03)
1000.88 (0.63)5.7 (0.85)
5205.5 (2.19)5.2 (0.72)
1006.0 (2.33)4.6 (0.20)
b. Twofold population contraction
0.210200.11 (0.022)138 (186)
1000.23 (0.035)8.4 (2.34)
1200.93 (0.47)18.5 (5.28)
1000.93 (0.080)11.1 (0.98)
5203.3 (0.94)15.7 (2.93)
1005.6 (2.01)12.4 (2.26)
c. Twofold population expansion
0.210200.22 (0.066)8.9 (3.57)
1000.22 (0.076)11.0 (7.15)
1201.0 (0.18)9.7 (1.95)
1001.0 (0.22)9.3 (1.52)
5205.3 (1.19)9.6 (0.98)
1005.6 (1.43)9.0 (0.63)
d. Fourfold population expansion
0.220200.18 (0.018)17.2 (1.4)
1000.18 (0.016)23.4 (4.8)
1201.1 (0.12)19.3 (2.0)
1001.2 (0.13)17.0 (1.4)
5207.2 (2.35)21.7 (3.5)
1009.4 (4.10)20.5 (3.0)
  • Replicate data sets were generated by Monte Carlo simulation for three values of β and E(s) = 0.1 under four demographic models: (a) N1 = 50, N2 = 50; (b) N1 = 200, N2 = 100; (c) N1 = 50, N2 = 100; and (d) N1 = 50, N2 = 200. Other parameter values were t2 = 200 and f0 = 0.95. The data sets consisted of 40 alleles sampled for 10,000 neutral sites and 2000 selected sites. Each data set was analyzed separately assuming two different N1's (20 and 100), and estimates of N2, t, β, E(s), and f0 were obtained. There were five replicates for β = 1 and β = 5 and 10 replicates for β = 0.2.