MCMC for Bayesian uncertainty quantification from time-series data

05/28/2020
by   Philip Maybank, et al.
0

Many problems in science and engineering require uncertainty quantification that accounts for observed data. For example, in computational neuroscience, Neural Population Models (NPMs) are mechanistic models that describe brain physiology in a range of different states. Within computational neuroscience there is growing interest in the inverse problem of inferring NPM parameters from recordings such as the EEG (Electroencephalogram). Uncertainty quantification is essential in this application area in order to infer the mechanistic effect of interventions such as anaesthesia. This paper presents C++ software for Bayesian uncertainty quantification in the parameters of NPMs from approximately stationary data using Markov Chain Monte Carlo (MCMC). Modern MCMC methods require first order (and in some cases higher order) derivatives of the posterior density. The software presented offers two distinct methods of evaluating derivatives: finite differences and exact derivatives obtained through Algorithmic Differentiation (AD). For AD, two different implementations are used: the open source Stan Math Library and the commercially licenced dco/c++ tool distributed by NAG (Numerical Algorithms Group). The use of derivative information in MCMC sampling is demonstrated through a simple example, the noise-driven harmonic oscillator. And different methods for computing derivatives are compared. The software is written in a modular object-oriented way such that it can be extended to derivative based MCMC for other scientific domains.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
05/02/2018

BayesLands: A Bayesian inference approach for parameter uncertainty quantification in Badlands

Bayesian inference provides a principled approach towards uncertainty qu...
research
04/17/2018

On Barotropic Mechanisms of Uncertainty Propagation in Estimation of Drake Passage Transport

Uncertainty in estimation of Drake Passage transport is analyzed in a He...
research
06/18/2020

Bayesian Changepoint Analysis

In my PhD thesis, we elaborate upon Bayesian changepoint analysis, where...
research
07/29/2021

Nonparametric posterior learning for emission tomography with multimodal data

In this work we continue studies of the uncertainty quantification probl...
research
03/19/2020

Systematic statistical analysis of microbial data from dilution series

In microbial studies, samples are often treated under different experime...
research
11/13/2017

Uncertainty quantification for radio interferometric imaging: II. MAP estimation

Uncertainty quantification is a critical missing component in radio inte...

Please sign up or login with your details

Forgot password? Click here to reset