Species 1 × species 2
|
m
1
|
m
2
|
t (years)
|
Ne1
|
Ne2
|
NeA
|
2N1M1
|
2N2M2
|
---|
Hylobates
comparsions
|
H.lar × H. agilis
|
Peak value
|
0.0035
|
0.1025*
|
2487549
|
63334
|
42819
|
39479
|
0.0079
|
0.1009*
|
Lower 95% HPD
|
0.0005
|
0.0125
|
1686037
|
49260
|
30415
|
8468
|
0.0004
|
0.0153
|
Upper 95% HPD
|
0.1055
|
0.2605
|
3758517
|
80748
|
59278
|
82894
|
0.1435
|
0.2417
|
H. lar × H. moloch
|
Peak value
|
0.0528*
|
0.0003a
|
2315183
|
61664
|
32084
|
43296
|
0.0744*
|
0.0003a
|
Lower 95% HPD
|
0.0068
|
0.0003a
|
1557155
|
47590
|
21588
|
16817
|
0.0101
|
0.0003a
|
Upper 95% HPD
|
0.1522
|
0.1177
|
3319799
|
79555
|
47113
|
84326
|
0.2040
|
0.0836
|
H. lar × H. muelleri
|
Peak value
|
0.0898*
|
0.0003a
|
2205588
|
55342
|
161665
|
34623
|
0.1175*
|
0.0021a
|
Lower 95% HPD
|
0.0013
|
0.0003a
|
1685015
|
40076
|
97326
|
16630
|
0.0090
|
0.0021a
|
Upper 95% HPD
|
0.2557
|
0.0873
|
2881421
|
73881
|
307791
|
61885
|
0.2914
|
0.3667
|
H. lar × H. pileatus
|
Peak value
|
0.0003a
|
0.0003a
|
3000898
|
63214
|
11075
|
67645
|
0.0002a
|
0.0002a
|
Lower 95% HPD
|
0.0003a
|
0.0003a
|
1770005
|
49583
|
5623
|
13461
|
0.0002a
|
0.0002a
|
Upper 95% HPD
|
0.0411
|
0.2091
|
5624895
|
79913
|
18913
|
141253
|
0.0577
|
0.0508
|
H. moloch × H. agilis
|
Peak value
|
0.0005a
|
0.0545*
|
2138864
|
30608
|
50112
|
22304
|
0.0004a
|
0.0650*
|
Lower 95% HPD
|
0.0005a
|
0.0015
|
1318541
|
23180
|
36015
|
96.55
|
0.0004a
|
0.0018
|
Upper 95% HPD
|
0.1455
|
0.2395
|
2968456
|
45863
|
69229
|
53008
|
0.0986
|
0.2593
|
H. moloch × H. muelleri
|
Peak value
|
0.0516
|
0.0396
|
1520623
|
29171
|
215645
|
34623
|
0.0423
|
0.2063
|
Lower 95% HPD
|
0.0004
|
0.0004
|
1100511
|
17175
|
122408
|
17175
|
0.0006
|
0.0048
|
Upper 95% HPD
|
0.3740
|
0.2916
|
2169057
|
46073
|
423929
|
60795
|
0.2239
|
1.6740
|
H. muelleri × H. agilis
|
Peak value
|
0.0006a
|
0.2586*
|
1867672
|
229719
|
37929
|
23616
|
0.0051a
|
0.2606*
|
Lower 95% HPD
|
0.0006a
|
0.0330
|
1429294
|
137163
|
22185
|
8349
|
0.0051a
|
0.0584
|
Upper 95% HPD
|
0.1710
|
0.8238
|
2506973
|
416738
|
58444
|
47948
|
1.0010
|
0.5714
|
Nomascus
comparisons
|
N. gabriellae × N. siki
|
Peak value
|
0.0730
|
0.1350
|
1768642
|
26922
|
21696
|
37599
|
0.1002
|
0.1790
|
Lower 95% HPD
|
0.0010
|
0.0010
| |
13063
|
10110
| |
0.0018
|
0.0020
|
Upper 95% HPD
|
0.9590
|
1.1910
| |
49640
|
42597
| |
0.5158
|
0.5411
|
N. leucogenys × N. gabriellae
|
Peak value
|
0.0053
|
0.1477
|
1676632
|
97877
|
27950
|
22554
|
0.0153
|
0.1512
|
Lower 95% HPD
|
0.0008
|
0.0008
|
1093218
|
63345
|
13489
|
2698
|
0.0022
|
0.0017
|
Upper 95% HPD
|
0.2542
|
0.8363
|
2614455
|
155503
|
53417
|
48669
|
0.5773
|
0.4736
|
N. leucogenys × N. siki
|
Peak value
|
0.0028a
|
1.1960*
|
1414914
|
102194
|
10468
|
15216
|
0.0061a
|
0.3641*
|
Lower 95% HPD
|
0.0028a
|
0.2447
|
864215
|
64856
|
3561
|
108
|
0.0061a
|
0.1268
|
Upper 95% HPD
|
0.5803
|
3.9410
|
2309117
|
166510
|
25575
|
37878
|
1.4320
|
0.8551
|
- Missing values are where the HPD of the parameters could not be reliably estimated by IMa2; m1 is migration rate into species 1 from species 2; m2 is migration rate into species 2 from species 1; t is the time since the species 1 and 2 split in years; Ne1, Ne2 and NeA are effective population sizes for species 1, species 2 and their ancestral population, respectively; 2N1M1 is population migration rate into species 1 from species 2 per generation; 2N2M2 is population migration rate into species 2 from species 1 per generation; aThe value corresponds to the first bin of the parameter space and hence could be interpreted as zero. *Asterisks indicates the estimates of migration rate that are significantly different from zero by the LLR tests [49, 63] at the P < 0.05 level and bold text indicates the significance after Bonferroni correction at the P < 0.0035 level (for the Hylobates comparisons) or at the P < 0.008 level (for the Nomascus comparisons).