Repository logo
Log in(current)
Repository logoMIT Open ScholarshipDSpace@MIT
  1. Home
  2. MIT Open Access Articles
  3. MIT Open Access Articles
  4. Likelihood-informed dimension reduction for nonlinear inverse problems

Likelihood-informed dimension reduction for nonlinear inverse problems

Thumbnail Image
Download
Name

Marzouk_Likelihood-informed.pdf

Size

1.89 MB

Format

Adobe PDF

Checksum (MD5)

1825e54f043a55241c658ed504a0c0e4

Author(s)
Martin, J.
•
Cui, Tiangang
•
Marzouk, Youssef M.
•
Solonen, Antti
•
Spantini, Alessio
Date Issued
October 2014
Journal
Inverse Problems
Publisher
IOP Publishing
Citation
Cui, T, J Martin, Y M Marzouk, A Solonen, and A Spantini. “Likelihood-Informed Dimension Reduction for Nonlinear Inverse Problems.” Inverse Problems 30, no. 11 (October 29, 2014): 114015.
Version
Author's final manuscript
Abstract
The intrinsic dimensionality of an inverse problem is affected by prior information, the accuracy and number of observations, and the smoothing properties of the forward operator. From a Bayesian perspective, changes from the prior to the posterior may, in many problems, be confined to a relatively low-dimensional subspace of the parameter space. We present a dimension reduction approach that defines and identifies such a subspace, called the 'likelihood-informed subspace' (LIS), by characterizing the relative influences of the prior and the likelihood over the support of the posterior distribution. This identification enables new and more efficient computational methods for Bayesian inference with nonlinear forward models and Gaussian priors. In particular, we approximate the posterior distribution as the product of a lower-dimensional posterior defined on the LIS and the prior distribution marginalized onto the complementary subspace. Markov chain Monte Carlo sampling can then proceed in lower dimensions, with significant gains in computational efficiency. We also introduce a Rao−Blackwellization strategy that de-randomizes Monte Carlo estimates of posterior expectations for additional variance reduction. We demonstrate the efficiency of our methods using two numerical examples: inference of permeability in a groundwater system governed by an elliptic PDE, and an atmospheric remote sensing problem based on Global Ozone Monitoring System (GOMOS) observations.
MIT Department
Massachusetts Institute of Technology. Department of Aeronautics and Astronautics
Terms of Use
Creative Commons Attribution-Noncommercial-Share Alike
http://creativecommons.org/licenses/by-nc-sa/4.0/
Persistent DSpace Link
http://hdl.handle.net/1721.1/96973
DOI of Published Version
https://doi.org/10.1088/0266-5611/30/11/114015
Repository logo
PrivacyPermissionsAccessibilityContact us
Repository logo
Notify us about copyright concerns.