Quick search Find article
Quick search
Find article

Fast inference for statistical inverse problems

Matthew A Taddy1, Herbert K H Lee2 and Bruno Sansó2

Show affiliations


Computer models for the simulation of physical and environmental phenomena are often regulated by complicated dependences on unknown variables, and these unobservable inputs must be inferred from a comparison of simulator output against physical data. The standard Bayesian statistical approaches to this inference problem require fitting a complicated statistical model to the existing parameter evaluations, usually through use of a Markov chain Monte Carlo sampling scheme. When there already exists a large bank of simulated values, it may be undesirable to develop a sophisticated statistical surrogate or to sample additional output from the computer simulator. In response to this motivation, we discuss a sampling importance resampling algorithm for Bayesian inference in inverse problems that works in conjunction with kernel density estimation to resample, from the original computer output, an approximate posterior sample for the unobservable variables of interest. Given a sufficiently large bank of computer output, our resampling method is able to provide high-quality results at a much lower cost than the standard Bayesian techniques. We present two applications where unobservable inputs are to be inferred from scarce observations and abundant simulated output. One consists of a climate simulator and the other of a groundwater flow model.


PACS

02.30.Zz Inverse problems

92.60.Ry Climatology

02.50.Ng Distribution theory and Monte Carlo studies

92.40.Kf Groundwater

02.50.Ga Markov processes

02.50.Tt Inference methods

MSC

86A22 Inverse problems (See also 35R30)

86A32 Geostatistics

62G09 Resampling methods

62F15 Bayesian inference

86A05 Hydrology, hydrography, oceanography (See also 76Bxx, 76E20, 76Q05, 76Rxx, 76U05)

60Jxx Markov processes

Subjects

Mathematical physics

Computational physics

Environmental and Earth science

Dates

Issue 8 (August 2009)

Received 9 October 2008, in final form 3 June 2009

Published 24 June 2009



  1. Fast inference for statistical inverse problems

    Matthew A Taddy et al 2009 Inverse Problems 25 085001

  2. Perturbative calculation of two-photon double electron ionization of helium

    I A Ivanov and A S Kheifets 2008 J. Phys. B: At. Mol. Opt. Phys. 41 095002

  3. PSU(2, 2|4) transformations of IIB superstring in AdS5 × S5

    Madoka Nishimura and Yoshiaki Tanii 2009 J. Phys. A: Math. Theor. 42 095401

  4. Author index with titles

    2009 Plasma Sources Sci. Technol. 18 049901

  5. A comparison of finite element and atomistic modelling of fracture

    V R Coffman et al 2008 Modelling Simul. Mater. Sci. Eng. 16 065008

  6. An electrostatic 3-phase linear stepper motor fabricated by vertical trench isolation technology

    Edin Sarajlic et al 2009 J. Micromech. Microeng. 19 074001

  7. Perturbations of the signum-Gordon model

    P Klimas 2008 J. Phys. A: Math. Theor. 41 095403

  8. Weyl functions, the inverse problem and special solutions for the system auxiliary to the nonlinear optics equation

    Alexander Sakhnovich 2008 Inverse Problems 24 025026

  9. Defects in semiconductors and oxides: where are the gaps in first principles theory?

    Marshall Stoneham 2009 Modelling Simul. Mater. Sci. Eng. 17 084009

  10. Virtues and flaws of the Pauli potential

    J Taruna et al 2008 J. Phys. A: Math. Theor. 41 035308

View by subject




Export








Please login to access our web services, or create an account if you don't yet have one.

You must have cookies enabled in your web browser to be able to login.

Username
Password

Forgotten your password? Get a new one here.