Higher order Quasi-Monte Carlo integration for Bayesian Estimation
METADATA ONLY
Loading...
Author / Producer
Date
2016
Publication Type
Report
ETH Bibliography
yes
Citations
Altmetric
METADATA ONLY
Data
Rights / License
Abstract
We analyze combined Quasi-Monte Carlo quadrature and Finite Element approximations in Bayesian estimation of solutions to countably-parametric operator equations with holomorphic dependence on the parameters as considered in [Cl. Schillings and Ch. Schwab: Sparsity in Bayesian Inversion of Parametric Operator Equations. Inverse Problems, 30, (2014)]. Such problems arise in numerical uncertainty quantification and in Bayesian inversion of operator equations with distributed uncertain inputs, such as uncertain coefficients, uncertain domains or uncertain source terms and boundary data. We show that the parametric Bayesian posterior densities belong to a class of weighted Bochner spaces of functions of countably many variables, with a particular structure of the QMC quadrature weights: up to a (problem-dependent, and possibly large) finite dimension S product weights can be used, and beyond this dimension, weighted spaces with so-called SPOD weights, recently introduced in [J. Dick, F.Y. Kuo, Q.T. Le Gia, D. Nuyens and Ch. Schwab, Christoph Higher order QMC Petrov-Galerkin discretization for affine parametric operator equations with random field inputs. SIAM J. Numer. Anal. 52 (2014), 2676--2702.], are used to describe the solution regularity. We establish error bounds for higher order Quasi-Monte Carlo quadrature for the Bayesian estimation based on [J. Dick, Q.T. LeGia and Ch. Schwab, Higher order Quasi-Monte Carlo integration for holomorphic, parametric operator equations, Report 2014-23, SAM, ETH Zürich]. It implies, in particular, regularity of the parametric solution and of the countably-parametric Bayesian posterior density in SPOD weighted spaces. This, in turn, implies that the Quasi-Monte Carlo quadrature methods in [J. Dick, F.Y. Kuo, Q.T. Le Gia, D. Nuyens, Ch. Schwab, Higher order QMC Galerkin discretization for parametric operator equations, SINUM (2014)] are applicable to these problem classes, with dimension-independent convergence rates O(N−1/p) of N-point HoQMC approximated Bayesian estimates, where 0
Permanent link
Publication status
published
Editor
Book title
Journal / series
Volume
2016-13
Pages / Article No.
Publisher
Seminar for Applied Mathematics, ETH Zurich
Event
Edition / version
Methods
Software
Geographic location
Date collected
Date created
Subject
Quasi Mont Carlo; Lattice rules; Digital nets; Parametric operator equations; Infinite-dimensional quadrature; Bayesian inverse problems; Uncertainty quantification; CBC construction; SPOD weights
Organisational unit
03435 - Schwab, Christoph / Schwab, Christoph