Abstract:It is often desirable to build a statistical emulator of a complex computer simulator in order to perform analysis which would otherwise be computationally infeasible. We propose methodology to model multivariate output from a computer simulator taking into account output structure in the responses. The utility of this approach is demonstrated by applying it to a chemical and biological hazard prediction model. Predicting the hazard area which results from an accidental or deliberate chemical or biological release is imperative in civil and military planning and also in emergency response. The hazard area resulting from such a release is highly structured in space, and we therefore propose the use of a thin-plate spline to capture the spatial structure and fit a Gaussian process emulator to the coefficients of the resultant basis functions. We compare and contrast four different techniques for emulating multivariate output: dimension reduction using (i) a fully Bayesian approach with a principal component basis, (ii) a fully Bayesian approach with a thin-plate spline basis, assuming that the basis coefficients are independent, and (iii) a plug-in Bayesian approach with a thin-plate spline basis and a separable covariance structure; and (iv) a functional data modeling approach using a tensor-product (separable) Gaussian process. We develop methodology for the two thin-plate spline emulators and demonstrate that these emulators significantly outperform the principal component emulator. Further, the separable thin-plate spline emulator, which accounts for the dependence between basis coefficients, provides substantially more realistic quantification of uncertainty, and is also computationally more tractable, allowing fast emulation. For high resolution output data, it also offers substantial predictive and computational advantages over the tensor-product Gaussian process emulator.

Emulation of Multivariate Simulators Using Thin-Plate Splines with Application to Atmospheric Dispersion

Residual Gaussian Process: A Tractable Nonparametric Bayesian Emulator for Multi-Fidelity Simulations

Neural Network Emulation of Spatio-temporal Data Using Linear and Nonlinear Dimensionality Reduction

Gaussian Process Emulation for High-Dimensional Coupled Systems

Evolving Bayesian Emulators for Structured Chaotic Time Series, with Application to Large Climate Models

Manifold Learning for the Emulation of Spatial Fields from Computational Models.

Spline-Based Bayesian Emulators for Large Scale Spatial Inverse Problems

Deep Gaussian Process Emulation and Uncertainty Quantification for Large Computer Experiments

Gaussian process emulation for discontinuous response surfaces with applications for cardiac electrophysiology models

Bayesian Emulation of Geotechnical Deterioration Curves Using Quadratic and B-Spline Hierarchical Models

Computer emulation with non-stationary Gaussian processes

Bayesian Emulation for Computer Models with Multiple Partial Discontinuities

Why not a thin plate spline for spatial models? A comparative study using Bayesian inference

Combining Data and Physical Models for Probabilistic Analysis: A Bayesian Augmented Space Learning Perspective

Computer Model Emulation with High-Dimensional Functional Output in Large-Scale Observing System Uncertainty Experiments

A Fast and Calibrated Computer Model Emulator: An Empirical Bayes Approach

Mesh-clustered Gaussian Process Emulator for Partial Differential Equation Boundary Value Problems

Emulation‐based inference for spatial infectious disease transmission models incorporating event time uncertainty

Emulation and History Matching Using the hmer Package

Reduced Dimensional Gaussian Process Emulators Of Parametrized Partial Differential Equations Based On Isomap

Embed and Emulate: Contrastive representations for simulation-based inference