Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2013 Feb;30(2):239-43.
doi: 10.1093/molbev/mss243. Epub 2012 Oct 22.

Accurate model selection of relaxed molecular clocks in bayesian phylogenetics

Guy Baele  1 , Wai Lok Sibon Li, Alexei J Drummond, Marc A Suchard, Philippe Lemey
Affiliations

Accurate model selection of relaxed molecular clocks in bayesian phylogenetics

Guy Baele et al. Mol Biol Evol. 2013 Feb.

Abstract

Recent implementations of path sampling (PS) and stepping-stone sampling (SS) have been shown to outperform the harmonic mean estimator (HME) and a posterior simulation-based analog of Akaike's information criterion through Markov chain Monte Carlo (AICM), in bayesian model selection of demographic and molecular clock models. Almost simultaneously, a bayesian model averaging approach was developed that avoids conditioning on a single model but averages over a set of relaxed clock models. This approach returns estimates of the posterior probability of each clock model through which one can estimate the Bayes factor in favor of the maximum a posteriori (MAP) clock model; however, this Bayes factor estimate may suffer when the posterior probability of the MAP model approaches 1. Here, we compare these two recent developments with the HME, stabilized/smoothed HME (sHME), and AICM, using both synthetic and empirical data. Our comparison shows reassuringly that MAP identification and its Bayes factor provide similar performance to PS and SS and that these approaches considerably outperform HME, sHME, and AICM in selecting the correct underlying clock model. We also illustrate the importance of using proper priors on a large set of empirical data sets.

PubMed Disclaimer

Figures

F<sc>ig</sc>. 1.
Fig. 1.
Results for the analysis of a large set of mammalian genes: (A) model selection results for the analysis of 872 genes using improper priors using the different estimators; (B) model selection results for the analysis of 961 genes using proper priors; (C) comparison of the log Bayes factor estimates for a common set of 872 genes for the MAP and sHME (gray) and for the MAP and SS (black), assuming improper priors; and (D) comparison of the log Bayes factor estimates for a common set of 961 genes for the MAP and sHME (gray) and for the MAP and SS (black), assuming proper priors.
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Baele G, Lemey P, Bedford T, Rambaut A, Suchard MA, Alekseyenko AV. Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty. Mol Biol Evol. 2012;29:2157–2167. - PMC - PubMed
    1. Carlin BP, Chib S. Bayesian model choice via Markov chain Monte Carlo. J R Stat Soc B. 1995;57:473–484.
    1. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–1973. - PMC - PubMed
    1. Friel N, Petitt AN. Marginal likelihood estimation via power posteriors. J R Stat Soc B. 2008;70:589–607.
    1. Gelman A, Meng XL. Simulating normalizing constants: from importance sampling to bridge sampling to path sampling. Stat Sci. 1998;13:163–185.

Publication types

-