Table 2

Maximum-likelihood estimates (MLEs) and the 90% highest posterior density (HPD) intervals of model parameters from IMa analyses
θ1 θ2 θA m1 m2 t N1 N2 NA m1 m2 2N1m1 2N2m2 T (year)
(1) A. linearis vs. (2) A. apiculata
  MLE 0.095 0.263 0.718 2.665 0.388 0.065 8643 23831 65189 0.007 0.001 0.127 0.051 24782
 HPD90Lo 0.040 0.132 0.004 0.655 0.005 0.018 4010 12483 9709 0.000 0.000 0.013 0.000 9074
 HPD90Hi 0.207 0.454 7.171 6.538 1.802 8.628 21473 45014 704721 0.023 0.008 0.677 0.409 3131867
(1) A. oblonga vs. (2) A. macroclinidioides (north)
 MLE 0.053 0.492 0.772 19.827 10.165 0.025 4579 44819 70084 0.057 0.023 0.526 2.499 9074
 HPD90Lo 0.022 0.223 0.006 4.777 2.925 0.005 1902 20152 501 0.009 0.005 0.053 0.326 1815
 HPD90Hi 0.170 0.912 9.641 48.573 23.715 9.985 15007 82918 873072 0.168 0.063 4.118 10.810 3624008
(1) North vs. (2) South
 MLE 0.436 0.585 0.158 0.005 0.208 1.605 60536 81181 21950 0.000 0.000 0.001 0.061 890679
 HPD90Lo 0.264 0.354 0.013 0.005 0.035 0.545 36578 49127 1742 0.000 0.000 0.001 0.006 302442
 HPD90Hi 0.668 0.927 12.359 0.215 0.488 26.835 92622 128566 1714557 0.001 0.001 0.072 0.226 14891815

N1, N2, NA: effective population sizes of population 1, 2 and ancestral. θ=4Nu.

m1: the rate of gene flow from population 2 to population 1 per gene per generation. m1=m1/u.

m2: the rate of gene flow from population 1 to population 2 per gene per generation. m2=m2/u.

2N1m1: the effective rate of gene flow from population 2 to population 1 per generation.

2N2m2: the effective rate of gene flow from population 1 to population 2 per generation.

Each value is the average of three independent runs of the simulations. Demographic quantities, N1, N2, NA, m1, m2, and T, were converted based on a mutation rate (u) of 7.49 × 10–9 substitutions/site/year.

Mitsui and Setoguchi

Mitsui and Setoguchi BMC Evolutionary Biology 2012 12:254   doi:10.1186/1471-2148-12-254

Open Data