Abstract:A new multi-level Markov Chain Monte Carlo algorithm for Approximate Bayesian Computation, ABC-SubSim, has recently appeared that exploits the Subset Simulation method for efficient rare-event simulation. ABC-SubSim adaptively creates a nested decreasing sequence of data-approximating regions in the output space that correspond to increasingly closer approximations of the observed output vector in this output space. At each level, multiple samples of the model parameter vector are generated by a component-wise Metropolis algorithm so that the predicted output corresponding to each parameter value falls in the current data-approximating region. Theoretically, if continued to the limit, the sequence of data-approximating regions would converge on to the observed output vector and the approximate posterior distributions, which are conditional on the data-approximation region, would become exact, but this is not practically feasible. In this paper we study the performance of the ABC-SubSim algorithm for Bayesian updating of the parameters of dynamical systems using a general hierarchical state-space model. We note that the ABC methodology gives an approximate posterior distribution that actually corresponds to an exact posterior where a uniformly distributed combined measurement and modeling error is added. We also note that ABC algorithms have a problem with learning the uncertain error variances in a stochastic state-space model and so we treat them as nuisance parameters and analytically integrate them out of the posterior distribution. In addition, the statistical efficiency of the original ABC-SubSim algorithm is improved by developing a novel strategy to regulate the proposal variance for the component-wise Metropolis algorithm at each level. We demonstrate that Self-regulated ABC-SubSim is well suited for Bayesian system identification by first applying it successfully to model updating of a two degree-of-freedom linear structure for three cases: globally, locally and un-identifiable model classes, and then to model updating of a two degree-of-freedom nonlinear structure with Duffing nonlinearities in its interstory force-deflection relationship.

Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems

Approximate Bayesian Computation by Subset Simulation for Parameter Inference of Dynamical Models

Inference for SDE models via Approximate Bayesian Computation

Approximate Bayesian Computation for Smoothing

ABC-CDE: Towards Approximate Bayesian Computation with Complex High-Dimensional Data and Limited Simulations

Lack of confidence in ABC model choice

Efficient exact inference for dynamical systems with noisy measurements using sequential approximate Bayesian computation

Approximate Bayesian Computation by Subset Simulation Using Hierarchical State-Space Models

Approximate Bayesian Forecasting

An automatic adaptive method to combine summary statistics in approximate Bayesian computation

Approximate Bayesian Computation via Classification

Stratified sampling and bootstrapping for approximate Bayesian computation

Parameter Estimation for Hidden Markov Models with Intractable Likelihoods

Accelerating inference for diffusions observed with measurement error and large sample sizes using Approximate Bayesian Computation

A comparison of emulation methods for Approximate Bayesian Computation

Calibrating the Discrete Boundary Conditions of a Dynamic Simulation: A Combinatorial Approximate Bayesian Computation Sequential Monte Carlo (ABC-SMC) Approach

Generalized Posteriors in Approximate Bayesian Computation

Monte Carlo algorithms for model assessment via conflicting summaries

Parallelization and High-Performance Computing Enables Automated Statistical Inference of Multi-scale Models.

An improved Approximate Bayesian Computation approach for high-dimensional posterior exploration of hydrological models

An easy-to-use empirical likelihood ABC method