Matthew A Taddy et al 2009 Inverse Problems 25 085001 doi:10.1088/0266-5611/25/8/085001
Matthew A Taddy1, Herbert K H Lee2 and Bruno Sansó2
Show affiliationsComputer 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.
86A22 Inverse problems (See also 35R30)
86A05 Hydrology, hydrography, oceanography (See also 76Bxx, 76E20, 76Q05, 76Rxx, 76U05)
Issue 8 (August 2009)
Received 9 October 2008, in final form 3 June 2009
Published 24 June 2009
Matthew A Taddy et al 2009 Inverse Problems 25 085001
Anton M. Zeitlin JHEP12(2007)098
F. Basile and G. Akemann JHEP12(2007)043
Maite López et al 2008 Metrologia 45 08003
Jérôme Martin et al JHEP03(2005)060
G Vargha et al 2008 Metrologia 45 08004
Alexander Oron et al 2007 J. Phys.: Conf. Ser. 64 012007
C Bambi and F R Urban 2008 Class. Quantum Grav. 25 095006
Valerio Marra and Francesca Rosati JCAP05(2005)011
B Joseph et al 2009 J. Phys.: Condens. Matter 21 432201