TABLE 6

Estimating ρ in the presence of sequencing errors

XNTρ0 = 20ρ0 = 40ρ0 = 60
16×321.2 (6.1)41.6 (11.1)62.4 (15.3)
521.2 (6.2)41.8 (11.0)62.9 (15.8)
721.1 (6.2)41.9 (11.0)62.9 (15.8)
321.4 (6.6)41.2 (10.2)61.5 (14.1)
519.9 (6.3)39.3 (11.2)58.9 (14.9)
720.3 (6.7)39.6 (11.4)57.7 (14.2)
320.7 (6.9)41.9 (12.0)60.6 (15.7)
519.9 (6.6)39.5 (11.3)58.2 (14.6)
718.8 (5.9)35.0 (11.2)51.4 (15.2)
319.0 (6.1)37.0 (10.6)53.6 (15.4)
512.5 (8.4)22.8 (19.9)27.6 (23.2)
7NANANA
  • Here, θ = 100 and NT = number of reads required to use homozygous genotype calls in the inference. NA, not applicable, due to there being too few two-locus pairs for which calls are made. We report mean (and RMSE) of the estimates over 1000 simulations.