ρ | ||||||
---|---|---|---|---|---|---|
3 | 1 | 0.5 | 0.1 | 0.05 | 0.01 | |
Actual recombination events ≥1 | ||||||
9990 | 9474 | 8035 | 2925 | 1680 | 46 | |
Detectable recombination events ≥1 | ||||||
θ = 7 | 6763 | 3401 | 1864 | 442 | 220 | 48 |
θ = 6 | 6270 | 3057 | 1712 | 381 | 186 | 38 |
θ = 5 | 5725 | 2667 | 1527 | 327 | 162 | 33 |
θ = 4 | 4868 | 2202 | 1171 | 249 | 119 | 28 |
θ = 3 | 3965 | 1733 | 866 | 182 | 85 | 22 |
θ = 2 | 2687 | 1073 | 548 | 108 | 62 | 10 |
θ = 1 | 1120 | 430 | 211 | 35 | 15 | 7 |
θ = 0.5 | 410 | 156 | 64 | 10 | 4 | 1 |
θ = 0.1 | 17 | 7 | 3 | 1 | 0 | 1 |
For each value of ρ, 10,000 realizations of the ancestral recombination graph for n = 20 were generated, and it was noted whether recombination had taken place in each. To each graph, infinite-sites mutations were then added, once for each value of θ, and it was noted whether four gametes could be detected or not. For example, the probability that a genealogy of n = 20 fragments contains recombination events is 8035/10,000 = 0.80 if ρ= 0.5 for the fragment. However, the probability that such a sample contains a detectable recombination event is only 1864/10,000 = 0.19 if θ = 7 and 211/10,000 = 0.02 if θ= 1.