User Manual Privacy Policy Disclaimer Contact us
  Advanced SearchBrowse





OBAMA: OBAMA for Bayesian amino-acid model averaging


Bouckaert,  Remco
Linguistic and Cultural Evolution, Max Planck Institute for the Science of Human History, Max Planck Society;

Tatarinova,  Tatiana
Linguistic and Cultural Evolution, Max Planck Institute for the Science of Human History, Max Planck Society;

External Ressource

(Supplementary material)

Fulltext (public)

(Publisher version), 827KB

Supplementary Material (public)
There is no public supplementary material available

Bouckaert, R. (2020). OBAMA: OBAMA for Bayesian amino-acid model averaging. PeerJ, 8: e9460. doi:10.7717/peerj.9460.

Cite as: http://hdl.handle.net/21.11116/0000-0006-D212-8
Background: Bayesian analyses offer many benefits for phylogenetic, and have been popular for analysis of amino acid alignments. It is necessary to specify a substitution and site model for such analyses, and often an ad hoc, or likelihood based method is employed for choosing these models that are typically of no interest to the analysis overall. - Methods: We present a method called OBAMA that averages over substitution models and site models, thus letting the data inform model choices and taking model uncertainty into account. It uses trans-dimensional Markov Chain Monte Carlo (MCMC) proposals to switch between various empirical substitution models for amino acids such as Dayhoff, WAG, and JTT. Furthermore, it switches base frequencies from these substitution models or use base frequencies estimated based on the alignment. Finally, it switches between using gamma rate heterogeneity or not, and between using a proportion of invariable sites or not. - Results: We show that the model performs well in a simulation study. By using appropriate priors, we demonstrate both proportion of invariable sites and the shape parameter for gamma rate heterogeneity can be estimated. The OBAMA method allows taking in account model uncertainty, thus reducing bias in phylogenetic estimates. The method is implemented in the OBAMA package in BEAST 2, which is open source licensed under LGPL and allows joint tree inference under a wide range of models.