TABLE 1

Fluctuate simulation results

Θ = 0.01Θ = 0.1
Locibp = 500bp = 1000bp = ∞bp = 500bp = 1000bp = ∞
A. Estimate of Θ
g = −1010.0140.0130.010NDNDND
20.0120.0110.010NDNDND
g = 010.0130.0120.0110.1120.1070.113
20.0120.0110.0100.1040.1060.104
g = 1001NDND0.0110.0970.0920.110
2NDND0.0110.1030.0970.106
B. Standard deviation of Θ
g = −1010.0090.0060.002NDNDND
20.0040.0030.002NDNDND
g = 010.0090.0040.0020.0320.0270.028
20.0040.0030.0020.0230.0170.018
g = 1001NDND0.0030.0430.0380.031
2NDND0.0020.0330.0290.021
C. Estimate of g
g = −101257.3165.850.0NDNDND
2130.271.327.8NDNDND
g = 01360.2145.945.114.66.212.7
2128.482.146.85.16.25.3
g = 1001NDND227.473.669.7119.1
2NDND187.153.452.7110.2
D. Standard deviation of Θ
g = −101567.1286.6116.9NDNDND
2463.2149.267.8NDNDND
g = 011215.1248.395.319.315.115.7
2298.5144.688.18.110.28.5
g = 1001NDND214.873.669.745.9
2NDND144.753.452.727.8
E. Number of samples (out of 100) in which the true values were rejected at the 95% level
g = −1011561NDNDND
27122NDNDND
g = 0116102243
21375913
g = 1001NDND6658
2NDND7743
  • Estimates of Θ and g based on 100 simulated data sets each, with 25 sequences of the given number of base pairs. Columns headed bp = ∞ were created by assuming that the genealogy could be reconstructed without error. Table 1E shows the number of times that the true values of Θ and g could be rejected at the nominal 95% level, out of 100 data sets. ND, not determined.