Skip to main content

Table 4 IMa2 estimates and 95% highest probability density intervals (HPD) of demographic parameters

From: Inferring the evolutionary histories of divergences in Hylobates and Nomascus gibbons through multilocus sequence data

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
  1. 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).