A Laplace method for under-determined Bayesian optimal experimental designs

Quan Long, Marco Scavino, Raul Tempone, Suojin Wang

Research output: Contribution to journalArticlepeer-review

20 Scopus citations

Abstract

In Long et al. (2013), a new method based on the Laplace approximation was developed to accelerate the estimation of the post-experimental expected information gains (Kullback–Leibler divergence) in model parameters and predictive quantities of interest in the Bayesian framework. A closed-form asymptotic approximation of the inner integral and the order of the corresponding dominant error term were obtained in the cases where the parameters are determined by the experiment. In this work, we extend that method to the general case where the model parameters cannot be determined completely by the data from the proposed experiments. We carry out the Laplace approximations in the directions orthogonal to the null space of the Jacobian matrix of the data model with respect to the parameters, so that the information gain can be reduced to an integration against the marginal density of the transformed parameters that are not determined by the experiments. Furthermore, the expected information gain can be approximated by an integration over the prior, where the integrand is a function of the posterior covariance matrix projected over the aforementioned orthogonal directions. To deal with the issue of dimensionality in a complex problem, we use either Monte Carlo sampling or sparse quadratures for the integration over the prior probability density function, depending on the regularity of the integrand function. We demonstrate the accuracy, efficiency and robustness of the proposed method via several nonlinear under-determined test cases. They include the designs of the scalar parameter in a one dimensional cubic polynomial function with two unidentifiable parameters forming a linear manifold, and the boundary source locations for impedance tomography in a square domain, where the unknown parameter is the conductivity, which is represented as a random field.
Original languageEnglish (US)
Pages (from-to)849-876
Number of pages28
JournalComputer Methods in Applied Mechanics and Engineering
Volume285
DOIs
StatePublished - Dec 19 2014

ASJC Scopus subject areas

  • General Physics and Astronomy
  • Mechanics of Materials
  • Mechanical Engineering
  • Computational Mechanics
  • Computer Science Applications

Fingerprint

Dive into the research topics of 'A Laplace method for under-determined Bayesian optimal experimental designs'. Together they form a unique fingerprint.

Cite this