BAMBI: blind accelerated multimodal Bayesian inference

Abstract

In this paper we present an algorithm for rapid Bayesian analysis that combines the benefits of nested sampling and artificial neural networks. The blind accelerated multimodal Bayesian inference (BAMBI) algorithm implements the MultiNest package for nested sampling as well as the training of an artificial neural network (NN) to learn the likelihood function. In the case of computationally expensive likelihoods, this allows the substitution of a much more rapid approximation in order to increase significantly the speed of the analysis. We begin by demonstrating, with a few toy examples, the ability of a NN to learn complicated likelihood surfaces. BAMBI's ability to decrease running time for Bayesian inference is then demonstrated in the context of estimating cosmological parameters from Wilkinson Microwave Anisotropy Probe and other observations. We show that valuable speed increases are achieved in addition to obtaining NNs trained on the likelihood functions for the different model and data combinations. These NNs can then be used for an even faster follow-up analysis using the same likelihood and different priors. This is a fully general algorithm that can be applied, without any pre-processing, to other problems with computationally expensive likelihood functions.

12 pages, 8 tables, 17 figures; accepted by MNRAS; v2 to reflect minor changes in published version

References in corpus (8)

• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Power Spectra and WMAP-Derived Parameters
• Bayesian Selection of sign(mu) within mSUGRA in Global Fits Including WMAP5 Results
• Use of the MultiNest algorithm for gravitational wave data analysis
• Fast cosmological parameter estimation using neural networks
• {\sc CosmoNet}: fast cosmological parameter estimation in non-flat models using neural networks
• Classifying LISA gravitational wave burst signals using Bayesian evidence
• Caching and Interpolated Likelihoods: Accelerating Cosmological Monte Carlo Markov Chains
• Pico: Parameters for the Impatient Cosmologist

Cited by in corpus (30)

• Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library
• Deep Neural Networks to Enable Real-time Multimessenger Astrophysics
• The First Two Years of Electromagnetic Follow-Up with Advanced LIGO and Virgo
• Bayesian analysis of cosmic-ray propagation: evidence against homogeneous diffusion
• Parameter estimation for binary neutron-star coalescences with realistic noise during the Advanced LIGO era
• Nested Sampling with Normalising Flows for Gravitational-Wave Inference
• Exploring the Bayesian parameter estimation of binary black holes with LISA
• Learning Bayesian posteriors with neural networks for gravitational-wave inference
• Nested Sampling Methods
• The Missing Link: Bayesian Detection and Measurement of Intermediate-Mass Black-Hole Binaries
• Bilby-MCMC: An MCMC sampler for gravitational-wave inference
• Gravitational wave parameter estimation with compressed likelihood evaluations
• Accelerating parameter inference with graphics processing units
• Applications of physics informed neural operators
• Importance nested sampling with normalising flows
• Cosmological parameter estimation via iterative emulation of likelihoods
• ECoPANN: A Framework for Estimating Cosmological Parameters using Artificial Neural Networks
• Accelerated Bayesian inference using deep learning
• Neural Network Reconstruction of $H^{\prime }(z)$ and its application in Teleparallel Gravity
• Exploring Photometric Redshifts as an Optimization Problem: An Ensemble MCMC and Simulated Annealing-Driven Template-Fitting Approach
• Neural Networks Optimized by Genetic Algorithms in Cosmology
• Density estimation with Gaussian processes for gravitational-wave posteriors
• Machine Learning Model of the Swift/BAT Trigger Algorithm for Long GRB Population Studies
• Detection of Dipole Modulation in CMB Temperature Anisotropy Maps from WMAP and Planck using Artificial Intelligence
• Deep Learning and genetic algorithms for cosmological Bayesian inference speed-up
• Estimation of Full Sky Power Spectrum between Intermediate to Large Angular Scales from Partial Sky CMB Anisotropies using Artificial Neural Network
• Predicting the motion of a high-Q pendulum subject to seismic perturbations using machine learning
• ParamANN: A Neural Network to Estimate Cosmological Parameters for $\Lambda$CDM Universe Using Hubble Measurements
• Optical Follow-up of Planck Cluster Candidates with Small Instruments
• Statistically-informed deep learning for gravitational wave parameter estimation