English

A Bayesian Approach for Parameter Estimation and Prediction using a Computationally Intensive Model

Data Analysis, Statistics and Probability 2015-02-06 v2 Nuclear Theory

Abstract

Bayesian methods have been very successful in quantifying uncertainty in physics-based problems in parameter estimation and prediction. In these cases, physical measurements y are modeled as the best fit of a physics-based model η(θ)\eta(\theta) where θ\theta denotes the uncertain, best input setting. Hence the statistical model is of the form y=η(θ)+ϵy = \eta(\theta) + \epsilon, where ϵ\epsilon accounts for measurement, and possibly other error sources. When non-linearity is present in η()\eta(\cdot), the resulting posterior distribution for the unknown parameters in the Bayesian formulation is typically complex and non-standard, requiring computationally demanding computational approaches such as Markov chain Monte Carlo (MCMC) to produce multivariate draws from the posterior. While quite generally applicable, MCMC requires thousands, or even millions of evaluations of the physics model η()\eta(\cdot). This is problematic if the model takes hours or days to evaluate. To overcome this computational bottleneck, we present an approach adapted from Bayesian model calibration. This approach combines output from an ensemble of computational model runs with physical measurements, within a statistical formulation, to carry out inference. A key component of this approach is a statistical response surface, or emulator, estimated from the ensemble of model runs. We demonstrate this approach with a case study in estimating parameters for a density functional theory (DFT) model, using experimental mass/binding energy measurements from a collection of atomic nuclei. We also demonstrate how this approach produces uncertainties in predictions for recent mass measurements obtained at Argonne National Laboratory (ANL).

Keywords

Cite

@article{arxiv.1407.3017,
  title  = {A Bayesian Approach for Parameter Estimation and Prediction using a Computationally Intensive Model},
  author = {Dave Higdon and Jordan D. McDonnell and Nicolas Schunck and Jason Sarich and Stefan M. Wild},
  journal= {arXiv preprint arXiv:1407.3017},
  year   = {2015}
}

Comments

19 pages, 10 figures, 1 table; Invited paper for the Journal of Physics G: Nuclear and Particle Physics focus section entitled "Enhancing the interaction between nuclear experiment and theory through information and statistics". In press

R2 v1 2026-06-22T05:01:30.813Z