Emulators for Multivariate Deterministic Functions

One of the challenges with emulating the response of a multivariate function to its inputs is the quantity of data that must be assimilated, which is the product of the number of model evaluations and the number of outputs. This paper shows how even large calculations can be made tractable. It is already appreciated that gains can be made when the emulator residual covariance function is treated as separable in the model-inputs and model-outputs. Here an additional simplification on the structure of the regressors in the emulator mean function allows very substantial further gains. The result is that it is now possible to emulate rapidly—on a desktop computer—models with hundreds of evaluations and hundreds of outputs. This is demonstrated through calculating costs in floating-point operations, and in an illustration. Even larger sets of outputs are possible if they have additional structure, e.g. spatial-temporal.

[1]  Holger Wendland,et al.  Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree , 1995, Adv. Comput. Math..

[2]  A. Seheult,et al.  Pressure Matching for Hydrocarbon Reservoirs: A Case Study in the Use of Bayes Linear Strategies for Large Computer Experiments , 1997 .

[3]  Michael Goldstein,et al.  Bayesian Forecasting for Complex Systems Using Computer Simulators , 2001 .

[4]  M. J. Bayarri,et al.  Computer model validation with functional output , 2007, 0711.3271.

[5]  Gene H. Golub,et al.  Matrix computations , 1983 .

[6]  Astrid Maute,et al.  Emulating the Thermosphere-Ionosphere Electrodynamics General Circulation Model (TIE-GCM) , 2007 .

[7]  A. OHagan,et al.  Bayesian analysis of computer code outputs: A tutorial , 2006, Reliab. Eng. Syst. Saf..

[8]  T. J. Mitchell,et al.  Bayesian Prediction of Deterministic Functions, with Applications to the Design and Analysis of Computer Experiments , 1991 .

[9]  D. Higdon,et al.  Computer Model Calibration Using High-Dimensional Output , 2008 .

[10]  D. Nychka,et al.  Covariance Tapering for Interpolation of Large Spatial Datasets , 2006 .

[11]  J. Rougier,et al.  Bayes Linear Calibrated Prediction for Complex Systems , 2006 .

[12]  A. O'Hagan,et al.  Bayesian calibration of computer models , 2001 .

[13]  Thomas J. Santner,et al.  The Design and Analysis of Computer Experiments , 2003, Springer Series in Statistics.

[14]  Michael Goldstein,et al.  Reified Bayesian modelling and inference for physical systems , 2009 .

[15]  M. Genton Separable approximations of space‐time covariance matrices , 2007 .