The Ball Pit Algorithm: A Markov Chain Monte Carlo Method Based on Path Integrals

11/05/2021
by   Miguel Fudolig, et al.
0

The Ball Pit Algorithm (BPA) is a novel Markov chain Monte Carlo (MCMC) algorithm for sampling marginal posterior distributions developed from the path integral formulation of the Bayesian analysis for Markov chains. The BPA yielded comparable results to the Hamiltonian Monte Carlo as implemented by the adaptive No U-Turn Sampler (NUTS) in sampling posterior distributions for simulated data from Bernoulli and Poisson likelihoods. One major advantage of the BPA is its significantly lower computational time, which was measured to be at least 95 also applied to a multi-parameter Cauchy model using real data of the height differences of cross- and self-fertilized plants. The posterior medians for the location parameter were consistent with other Bayesian sampling methods. Additionally, the posterior median for the logarithm of the scale parameter obtained from the BPA was close to the estimated posterior median calculated using the Laplace normal approximation. The computational time of the BPA implementation of the Cauchy analysis is 55 Overall, we have found that the BPA is a highly efficient alternative to the Hamiltonian Monte Carlo and other standard MCMC methods.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
11/30/2020

A Markov Chain Monte-Carlo Approach to Dose-Response Optimization Using Probabilistic Programming (RStan)

A hierarchical logistic regression Bayesian model is proposed and implem...
research
08/01/2022

Computing Bayes: From Then 'Til Now'

This paper takes the reader on a journey through the history of Bayesian...
research
11/27/2020

Comparison of Bayesian Nonparametric Density Estimation Methods

In this paper, we propose a nonparametric Bayesian approach for Lindsey ...
research
09/21/2022

Single chain differential evolution Monte-Carlo for self-tuning Bayesian inference

1. Bayesian inference is difficult because it often requires time consum...
research
06/05/2023

Gibbs Sampling the Posterior of Neural Networks

In this paper, we study sampling from a posterior derived from a neural ...
research
06/13/2019

Direct Sampling of Bayesian Thin-Plate Splines for Spatial Smoothing

Radial basis functions are a common mathematical tool used to construct ...
research
08/30/2018

An Introduction to Inductive Statistical Inference -- from Parameter Estimation to Decision-Making

These lecture notes aim at a post-Bachelor audience with a backgound at ...

Please sign up or login with your details

Forgot password? Click here to reset