Research reports

Sparsity in Bayesian Inversion of Parametric Operator Equations

by Cl. Schillings and Ch. Schwab

(Report number 2013-17)

Abstract
We establish posterior sparsity in Bayesian inversion for systems with distributed parameter uncertainty subject to noisy data. We generalize the particular case of scalar diffusion problems with random coefficients in [29] to broad classes of operator equations. For countably parametric, deterministic representations of uncertainty in the forward problem which belongs to a certain sparsity class, we quantify analytic regularity of the (countably parametric) Bayesian posterior density and prove that the parametric, deterministic density of the Bayesian posterior belongs to the same sparsity class. Generalizing [32, 29], the considered forward problems are parametric, deterministic operator equations, and computational Bayesian inversion is to evaluate expectations of quantities of interest (QoIs) under the Bayesian posterior, conditional on given data. The sparsity results imply, on the one hand, sparsity of Legendre (generalized) polynomial chaos expansions of the Bayesian posterior and, on the other hand, convergence rates for data-adaptive Smolyak integration algorithms for computational Bayesian estimation which are independent of dimension of the parameter space. The convergence rates are, in particular, superior to Markov Chain Monte-Carlo sampling of the posterior, in terms of the number N of instances of the parametric forward problem to be solved.

Keywords: Bayesian Inverse Problems, Parametric Operator Equations, Smolyak Quadrature, Sparsity, Uniform Prior Measures

BibTeX
@Techreport{SS13_513,
  author = {Cl. Schillings and Ch. Schwab},
  title = {Sparsity in Bayesian Inversion of Parametric Operator Equations},
  institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
  number = {2013-17},
  address = {Switzerland},
  url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2013/2013-17.pdf },
  year = {2013}
}

Disclaimer
© Copyright for documents on this server remains with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use. The administrators respectfully request that authors inform them when any paper is published to avoid copyright infringement. Note that unauthorised copying of copyright material is illegal and may lead to prosecution. Neither the administrators nor the Seminar for Applied Mathematics (SAM) accept any liability in this respect. The most recent version of a SAM report may differ in formatting and style from published journal version. Do reference the published version if possible (see SAM Publications).

JavaScript has been disabled in your browser