Isolation model fitting
Species 1 | Species 2 | Test | ![]() | ![]() | ![]() | T | Pχ2 | PWWH |
---|---|---|---|---|---|---|---|---|
sim | mau | 1 | 143.7 | 79.5 | 139.9 | 0.42 | 0.138 | 0.827 |
75.7-280.5 | 46.0-129.7 | 85.5-203.0 | 0.22-0.64 | |||||
sim | mau | 2 | 152.3 | 83.3 | 127.1 | 0.44 | 0.118 | 0.781 |
87.2-293.3 | 50.7-134.0 | 75.7-186.1 | 0.23-0.67 | |||||
sim | sec | 3 | 32.3 | 7.75 | 221.1 | 0.51 | 0.216 | 0.465 |
3.76-187.1 | 0.95-25.8 | 0.00-283.4 | 0.30-0.97 | |||||
sim | sec | 4 | 55.7 | 7.98 | 201.7 | 0.56 | 0.434 | 0.810 |
17.1-303.7 | 2.33-17.0 | 0.0-259.3 | 0.35-1.52 |
For each contrast, the data were fit to the isolation model as described (Wang et al. 1997). The estimated values for the primary parameters are shown, along with the 95% confidence intervals determined by simulation (see text). Four contrasts are shown: (1) 12 loci (all excluding ase and ci); (2) same as (1) but observed shared polymorphisms were decreased by that amount expected by chance (see Table 6); (3) 11 loci (all excluding ase, ci, and janus); and (4) same as (3) except that one D. sechellia sequence for In(2L)t data was excluded (see text). The P values, for both the χ2 and the Wang, Wakely, and Hey (WWH) test statistics, are the proportion of simulated values greater than or equal to the observed. The test is one-tailed because the focus is on detecting a departure from the model in the direction expected if historical gene flow had occurred.