Abstract
Computer simulators are widely used to understand complex physical systems in many areas such as aerospace, renewable energy, climate modelling, and manufacturing. One fundamental aspect of the study of computer simulators is known as experimental design, that is, how to select the input settings where the computer simulator is run and the corresponding response is collected. Extra care should be taken in the selection process because computer simulators can be computationally expensive to run. The selection should acknowledge and achieve the goal of the analysis. This article focuses on the goal of producing more accurate prediction which is important for risk assessment and decision making. We propose two new methods of design approaches that sequentially select input settings to achieve this goal. The approaches make novel applications of simultaneous and sequential contour estimations. Numerical examples are employed to demonstrate the effectiveness of the proposed approaches.
Similar content being viewed by others
References
Ba S (2015) SLHD: maximin-distance (sliced) Latin hypercube designs. R package version 2.1-1
Bayarri MJ, Berger JO, Calder ES, Dalbey K, Lunagomez S, Patra AK, Pitman EB, Spiller ET, Wolpert RL (2009) Using statistical and computer models to quantify volcanic hazards. Technometrics 51:402–413
Bower RG, Benson AJ et al (2006) The broken hierarchy of galaxy formation. Mon Not R Astron Soc 370:645–655
Chipman H, Ranjan P, Wang W (2012) Sequential design for computer experiments with a flexible Bayesian additive model. Can J Stat 40(4):663–678
Currin C, Mitchell T, Morris M, Ylvisaker D (1991) Bayesian prediction of deterministic functions, with applications to the design and analysis of computer experiments. J Am Stat Assoc 86(416):953–963
Dancik G (2018) mlegp: maximum likelihood estimates of Gaussian processes. R package version 3.1.7
Deng X, Lin CD, Liu KW, Rowe RK (2017) Additive Gaussian process for computer models with qualitative and quantitative factors. Technometrics 59:283–292
Dette H, Pepelyshev A (2010) Generalized Latin hypercube design for computer experiments. Technometrics 52:421–429
Dixon LCW, Szego GP (1978) The global optimization problem: an introduction. Towards Glob Optim 2:1–15
Fang KT, Li R, Sudjianto A (2005) Design and modeling for computer experiments. CRC Press, New York
Gramacy RB, Lee HKH (2008) Bayesian treed Gaussian process models with an application to computer modeling. J Am Stat Assoc 103(483):1119–1130
Gramacy RB, Lee HKH (2009) Adaptive design and analysis of supercomputer experiments. Technometrics 51(2):130–145
Gramacy RB, Lee HK (2012) Cases for the nugget in modeling computer experiments. Stat Comput 22(3):713–722
Gramacy RB, Taddy MA (2016) tgp: Bayesian treed Gaussian process models. R package version 2-4-14
Greenberg D (1979) A numerical model investigation of tidal phenomena in the Bay of Fundy and Gulf of Maine. Mar Geod 2:161–187
Gu M, Palomo J, Berger J (2018) RobustGaSP: robust Gaussian stochastic process emulation. R package version 0.5.6
Iman RL, Conover WJ (1982) A distribution-free approach to inducing rank correlation among input variables. Commun Stat Part B Simul Comput 11:311–334
Johnson M, Moore L, Ylvisaker D (1990) Minimax and maximin distance design. J Stat Plan Inference 26:131–148
Joseph VR, Dasgupta T, Tuo R, Wu CJ (2015) Sequential exploration of complex surfaces using minimum energy designs. Technometrics 57:64–74
Joseph VR, Gul E, Ba S (2015) Maximum projection designs for computer experiments. Biometrika 102:371–380
Joseph VR, Hung Y (2008) Orthogonal-maximin Latin hypercube designs. Stat Sin 18:171–186
Lam CQ, Notz WI (2008) Sequential adaptive designs in computer experiments for response surface model fit. Stat Appl 6:207–233
Linkletter C, Bingham D, Hengartner N, Higdon D, Ye KQ (2006) Variable selection for Gaussian process models in computer experiments. Technometrics 48(4):478–490
Loeppky JL, Moore LM, Williams BJ (2010) Batch sequential designs for computer experiments. J Stat Plan Inference 140(6):1452–1464
Loeppky JL, Sacks J, Welch WJ (2009) Choosing the sample size of a computer experiment: a practical guide. Technometrics 51(4):366–376
MacDonald B, Ranjan P, Chipman H (2015) GPfit: an R package for Gaussian process model fitting using a new optimization algorithm. R package version 1.0-1
McKay MD, Beckman RJ, Conover WJ (1979) A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics 21:239–45
Morris MD, Mitchell TJ (1995) Exploratory designs for computer experiments. J Stat Plan Inference 43:381–402
Owen AB (1992) Orthogonal arrays for computer experiments, integration and visualization. Stat Sin 2:439–452
Palomo J, Paulo R, Garcia-Donato G (2015) SAVE: an r package for the statistical analysis of computer models. J Stat Softw 64(13):1–23
Park JS (1991) Tuning complex computer codes to data and optimal designs. Unpublished Ph.D. thesis, University of Illinois, Champaign-Urbana
Roustant O, Ginsbourger D, Deville Y (2018) DiceKriging: kriging methods for computer experiments. R package version 1.5.6
Ranjan P, Bingham D, Michailidis G (2008) Sequential experiment design for contour estimation from complex computer codes. Technometrics 50(4):527–541 (Errata (2011), Technometrics, 53:1, 109–110)
Rasmussen CE, Williams CKI (2006) Gaussian processes for machine learning. The MIT Press, Cambridge, MA
Santner TJ, Williams BJ, Notz WI (2003) The design and analysis of computer experiments. Springer, New York
Sacks J, Schiller SB, Welch WJ (1989) Designs for computer experiments. Technometrics 31:41–47
Sacks J, Schiller S (1988) Spatial designs. In: Gupta, Berger (eds) Statistical decision theory and related topics IV, vol 2. Springer, New York, pp 385–399
Sacks J, Welch WJ, Mitchell TJ, Wynn HP (1989) Design and analysis of computer experiments. Stat Sci 4:409–423
Tang B (1993) Orthogonal array-based Latin hypercubes. J Am Stat Assoc 88:1392–1397
Zhang R, Lin CD, Ranjan P (2018) Local approximate Gaussian process model for large-scale dynamic computer experiments. J Comput Graph Stat 27:798–807
Acknowledgements
The authors would like to thank the chief editor, the handling editor, and the reviewers for their helpful comments. Yang’s research was funded by China Scholarship Council. Lin’s research was funded in part by the Discovery Grant from the Natural Sciences and Engineering Research Council of Canada. Ranjan’s research was partially supported by the Extra Mural Research Fund (EMR/2016/003332/MS) from the Science and Engineering Research Board (SERB), Government of India.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Part of special issue guest edited by Pritam Ranjan and Min Yang—Algorithms, Analysis and Advanced Methodologies in the Design of Experiments.
Rights and permissions
About this article
Cite this article
Yang, F., Lin, C.D. & Ranjan, P. Global Fitting of the Response Surface via Estimating Multiple Contours of a Simulator. J Stat Theory Pract 14, 9 (2020). https://doi.org/10.1007/s42519-019-0077-0
Published:
DOI: https://doi.org/10.1007/s42519-019-0077-0