Comparison of the exact CLR(t) (obtained by solving the full recursions numerically) with the analytic approximation (33) when the selected locus is between the two neutral loci
| rLS = rRS = 0.01s | rLS = rRS = 0.05s | rLS = rRS = 0.1s | ||||
|---|---|---|---|---|---|---|
| t | Exact | Equation 33 | Exact | Equation 33 | Exact | Equation 33 |
| 0 | 2.42 | 2.42 | 2.42 | 2.42 | 2.42 | 2.42 |
| 500 | 2.21 | 2.21 | 1.48 | 1.48 | 0.90 | 0.90 |
| 1000 | 2.18 | 2.18 | 0.98 | 0.98 | 0.36 | 0.36 |
| 1500 | 3.67 | 3.70 | 1.11 | 1.12 | 0.25 | 0.25 |
| 2000 | 6.91 | 6.89 | 1.39 | 1.39 | 0.19 | 0.19 |
| 2500 | 1.59 | 1.53 | 0.21 | 0.21 | 0.02 | 0.02 |
| 3000 | 0.14 | 0.13 | 0.01 | 0.01 | 0.00 | 0.00 |
| 3500 | 0.01 | 0.01 | 0.00 | 0.00 | 0.00 | 0.00 |
| 3982 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
The values shown are in units of 10−2. We used s = 0.01, pS(0) = 0.00005, pL(0) = 0.38, pR(0) = 0.41, and CLR(0) = 0.0242. In the exact numerical computation, pS(t) = 1 – pS(0) at t = 3982.