Skip to main content

Table 1 PAML results: log likelihoods, parameter estimates, Bayes Empirical Bayes identified sites, and likelihood ratio tests.

From: Long-wavelength sensitive visual pigments of the guppy (Poecilia reticulata): six opsins expressed in a single individual

Selected Branch

Model

n.p.

ln L

Parameter Estimates

Bayes Empirical Bayes Identified Sites *

Likelihood Ratio Test P-value

N/A

M0 (one-ratio)

33

-2127.37

ω = 0.08737

none

N/A

 

Site-specific models

 

M1a (neutral)

34

-2097.31

p0 = 0.91742, ω0 = 0.05054

N/A

N/A

 

M2a (selection)

36

-2097.31

p0 = 0.91768, ω0 = 0.05052, p2 = 0.04090, ω2 = 0.99980

275

P = 1 (vs. m1a)

 

M3 (discrete, K = 3)

37

-2088.38

p0 = 0.53422, p1 = 0.38154, ω0 = 0.00000, ω1 = 0.12511, ω2 = 0.63618

none

P < 0.000001 (vs. m0)

 

M7 (β)

34

-2089.48

p = 0.21896, q = 1.82556

N/A

N/A

 

M8 (β & ω)

36

-2088.91

p = 0.24125, q = 2.23372, p0 = 0.99249, ω = 1.66242

275

P = 0.568523 (vs. m7)

 

Branch-specific models

X 'green'

two-ratio

34

-2123.10

ω0 = 0.0806, ω1 = 0.4259

N/A

P = 0.003463 (vs. m0)

Y 'transducin'

two-ratio

34

-2125.68

ω0 = 0.0832, ω1 = 0.4523

N/A

P = 0.066388 (vs. m0)

 

Branch-site (BS) models

X 'green'

BS model A (null)

35

-2094.55

p0 = 0.73148, ω0 = 0.04521, p1 = 0.06817

N/A

N/A

 

BS model A (alternative)

36

-2093.93

p0 = 0.86007, ω0 = 0.04576, p1 = 0.06817, ω2 = 16.53309

271, 275, 277, 279, 285, 286, 302

P = 0.266644 (vs. BS model A [null])

Y 'transducin'

BS model A (null)

35

-2093.10

p0 = 0.31761, ω0 = 0.04432, p1 = 0.03013

N/A

N/A

 

BS model A (alternative)

36

-2093.10

p0 = 0.32723, ω0 = 0.04442, p1 = 0.03012, ω2 = 1.00000

248, 249, 256, 261, 263

P = 1 (vs. BS model A [null])

  1. n.p. = number of parameters (e.g., proportions, ω values, branch lengths).
  2. ln L = log likelihood value for model given the tree, dataset, and parameter estimates.
  3. * = sites in bold received posterior probability support > 0.9.
  4. Test statistic compared to ac2 distribution with d.o.f. = n.p.alternative - n.p.null
  5. See Fig. 4