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