Model selection and parameter inference in phylogenetics using Nested Sampling

Abstract

Bayesian inference methods rely on numerical algorithms for both model selection and parameter inference. In general, these algorithms require a high computational effort to yield reliable estimates. One of the major challenges in phylogenetics is the estimation of the marginal likelihood. This quantity is commonly used for comparing different evolutionary models, but its calculation, even for simple models, incurs high computational cost. Another interesting challenge relates to the estimation of the posterior distribution. Often, long Markov chains are required to get sufficient samples to carry out parameter inference, especially for tree distributions. In general, these problems are addressed separately by using different procedures. Nested sampling (NS) is a Bayesian computation algorithm which provides the means to estimate marginal likelihoods together with their uncertainties, and to sample from the posterior distribution at no extra cost. The methods currently used in phylogenetics for marginal likelihood estimation lack in practicality due to their dependence on many tuning parameters and the inability of most implementations to provide a direct way to calculate the uncertainties associated with the estimates. To address these issues, we introduce NS to phylogenetics. Its performance is assessed under different scenarios and compared to established methods. We conclude that NS is a competitive and attractive algorithm for phylogenetic inference. An implementation is available as a package for BEAST 2 under the LGPL licence, accessible at https://github.com/BEAST2-Dev/nested-sampling.

23 pages, 12 figures, 3 tables

References in corpus (9)

• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• Importance Nested Sampling and the MultiNest Algorithm
• PolyChord: next-generation nested sampling
• A Nested Sampling Algorithm for Cosmological Model Selection
• Model selection and parameter inference in phylogenetics using Nested Sampling
• Properties of Nested Sampling
• Fast Bayesian Inference for Exoplanet Discovery in Radial Velocity Data
• Foundations of Inference
• DNest4: Diffusive Nested Sampling in C++ and Python

Cited by in corpus (7)

• Model selection and parameter inference in phylogenetics using Nested Sampling
• Nested sampling for physical scientists
• Nested sampling cross-checks using order statistics
• The stepping-stone sampling algorithm for calculating the evidence of gravitational wave models
• Computational Techniques for Parameter Estimation of Gravitational Wave Signals
• Marginal likelihoods in phylogenetics: a review of methods and applications
• Bayesian support for Evolution: detecting phylogenetic signal in a subset of the primate family