A nonparametric ensemble transform method for Bayesian inference

From MaRDI portal
Publication:2855653

DOI10.1137/130907367zbMATH Open1362.65007arXiv1210.0375OpenAlexW2066601679MaRDI QIDQ2855653FDOQ2855653


Authors: Sebastian Reich Edit this on Wikidata


Publication date: 28 October 2013

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

Abstract: Many applications, such as intermittent data assimilation, lead to a recursive application of Bayesian inference within a Monte Carlo context. Popular data assimilation algorithms include sequential Monte Carlo methods and ensemble Kalman filters (EnKFs). These methods differ in the way Bayesian inference is implemented. Sequential Monte Carlo methods rely on importance sampling combined with a resampling step while EnKFs utilize a linear transformation of Monte Carlo samples based on the classic Kalman filter. While EnKFs have proven to be quite robust even for small ensemble sizes, they are not consistent since their derivation relies on a linear regression ansatz. In this paper, we propose another transform method, which does not rely on any a prior assumptions on the underlying prior and posterior distributions. The new method is based on solving an optimal transportation problem for discrete random variables.


Full work available at URL: https://arxiv.org/abs/1210.0375




Recommendations





Cited In (30)





This page was built for publication: A nonparametric ensemble transform method for Bayesian inference

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2855653)