Abstract:In Bayesian molecular-clock dating of species divergences, rate models are used to construct the prior on the molecular evolutionary rates for branches in the phylogeny, with independent and autocorrelated rate models being commonly used. The two class of models, however, can result in markedly different divergence time estimates for the same dataset, and thus Bayesian model selection appears necessary to select for the best rate model and obtain reliable inferences of divergence times. However, the properties of Bayesian rate model selection are not well understood, in particular when the number of sequence partitions analysed increases and when fossil calibrations are misspecified. Furthermore, Bayesian rate model selection is computationally expensive as it requires calculation of marginal likelihoods by MCMC sampling, and therefore methods that can speed up the model selection procedure without compromising its accuracy are desirable. In this study, we use a combination of computer simulations and real data analysis to investigate the statistical behavior of Bayesian rate model selection and we also explore approximations of the likelihood to improve computational efficiency in large phylogenomic datasets. Our simulations demonstrate that the posterior probability for the correct rate model converges to one as more molecular sequence partitions are analyzed and when no fossil calibrations are used, as expected due to asymptotic Bayesian model selection theory. Furthermore, we also show the model selection procedure is robust to slight misspecification of fossil calibrations, and reliable inference of the correct rate model is possible in this case. However, we show that when fossil calibrations are seriously misspecified, calculated model probabilities are completely wrong and may converge to one for the wrong rate model. Finally, we demonstrate that approximating the phylogenetic likelihood under an arcsine branch-length transform can dramatically reduce the computational cost of rate model selection without compromising accuracy. We test the approximate procedure on two large phylogenies of primates (372 species) and flowering plants (644 species), replicating results obtained on smaller datasets using exact likelihood. Our findings and methodology can assist users in selecting the optimal rate model for estimating times and rates along the Tree of Life.

The utility of a Bayesian analysis of complex models and the study of archeological ${}^{14}$C data

A Bayesian Approach to Inferring Depositional Ages Applied to a Late Tonian Reference Section in Svalbard

Bayesian radiocarbon modelling for beginners

Modelling diffusion of innovation curves using radiocarbon data

Discussion of: Bayesian views of an archaeological find

Central European Early Bronze Age chronology revisited: A Bayesian examination of large-scale radiocarbon dating

Bayesian Selection of Relaxed-clock Models: Distinguishing Between Independent and Autocorrelated Rates

Mixture Models for Dating with Confidence

Principles, Error Sources and Application Suggestions of Prevailing Molecular Dating Methods

Matching Bayesian and frequentist coverage probabilities when using an approximate data covariance matrix

A Bayesian hierarchical age model for optical dating of single grains of quartz

Probabilistic modelling improves relative dating from gene phylogenies

Quantifying Artifacts over Time: Interval Estimation of a Poisson Distribution using the Jeffreys Prior

Complex-valued Bayesian parameter estimation via Markov chain Monte Carlo

Bayesian classification for dating archaeological sites via projectile points

Technical note: A Bayesian mixing model to unravel isotopic data and quantify trace gas production and consumption pathways for time series data – Time-resolved FRactionation And Mixing Evaluation (TimeFRAME)

A Bayesian hierarchical age model for single-grain optical dating of feldspars

A macroevolutionary analysis of European Late Upper Palaeolithic stone tool shape using a Bayesian phylodynamic framework

Radiocarbon chronology of Iron Age Jerusalem reveals calibration offsets and architectural developments

Bayesian inference for a wavefront model of the Neolithisation of Europe

Relaxed Phylogenetics and Dating with Confidence