Constructing a simulation surrogate with partially observed output

Abstract

Gaussian process surrogates are a popular alternative to directly using computationally expensive simulation models. When the simulation output consists of many responses, dimension-reduction techniques are often employed to construct these surrogates. However, surrogate methods with dimension reduction generally rely on complete output training data. This article proposes a new Gaussian process surrogate method that permits the use of partially observed output while remaining computationally efficient. The new method involves the imputation of missing values and the adjustment of the covariance matrix used for Gaussian process inference. The resulting surrogate represents the available responses, disregards the missing responses, and provides meaningful uncertainty quantification. The proposed approach is shown to offer sharper inference than alternatives in a simulation study and a case study where an energy density functional model that frequently returns incomplete output is calibrated.

arXiv commentsubmitted to Technometrics

References in corpus (19)

• Scikit-learn: Machine Learning in Python
• Nuclear Energy Density Optimization
• Nuclear energy density optimization: Large deformations
• Error Estimates of Theoretical Models: a Guide
• Nuclear isotope shifts within the local energy-density functional approach
• Computer model validation with functional output
• Nuclear energy density optimization: Shell structure
• The Mira-Titan Universe II: Matter Power Spectrum Emulation
• Uncertainty Quantification for Nuclear Density Functional Theory and Information Content of New Measurements
• Towards the global description of nuclear charge radii: exploring the Fayans energy density functional
• Get on the BAND Wagon: A Bayesian Framework for Quantifying Model Uncertainties in Nuclear Dynamics
• Calibrating a large computer experiment simulating radiative shock hydrodynamics
• Energy Density Functional Approach to Superfluid Nuclei
• Uncertainty quantification for computer models with spatial output using calibration-optimal bases
• A Bayesian Approach for Parameter Estimation and Prediction using a Computationally Intensive Model
• Fast dimension-reduced climate model calibration and the effect of data aggregation
• Optimization and Supervised Machine Learning Methods for Fitting Numerical Physics Models without Derivatives
• Fast Nonseparable Gaussian Stochastic Process with Application to Methylation Level Interpolation
• Constructing a simulation surrogate with partially observed output

Cited by in corpus (3)

• On model emulation and closure tests for 3+1D relativistic heavy-ion collisions
• Constructing a simulation surrogate with partially observed output
• Portable, heterogeneous ensemble workflows at scale using libEnsemble