Skip to main content

Advertisement

Log in

Global Fitting of the Response Surface via Estimating Multiple Contours of a Simulator

  • Original Article
  • Published:
Journal of Statistical Theory and Practice Aims and scope Submit manuscript

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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9

Similar content being viewed by others

References

  1. Ba S (2015) SLHD: maximin-distance (sliced) Latin hypercube designs. R package version 2.1-1

  2. 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

    Article  MathSciNet  Google Scholar 

  3. Bower RG, Benson AJ et al (2006) The broken hierarchy of galaxy formation. Mon Not R Astron Soc 370:645–655

    Article  Google Scholar 

  4. 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

    Article  MathSciNet  Google Scholar 

  5. 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

    Article  MathSciNet  Google Scholar 

  6. Dancik G (2018) mlegp: maximum likelihood estimates of Gaussian processes. R package version 3.1.7

  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

    Article  MathSciNet  Google Scholar 

  8. Dette H, Pepelyshev A (2010) Generalized Latin hypercube design for computer experiments. Technometrics 52:421–429

    Article  MathSciNet  Google Scholar 

  9. Dixon LCW, Szego GP (1978) The global optimization problem: an introduction. Towards Glob Optim 2:1–15

    Google Scholar 

  10. Fang KT, Li R, Sudjianto A (2005) Design and modeling for computer experiments. CRC Press, New York

    Book  Google Scholar 

  11. Gramacy RB, Lee HKH (2008) Bayesian treed Gaussian process models with an application to computer modeling. J Am Stat Assoc 103(483):1119–1130

    Article  MathSciNet  Google Scholar 

  12. Gramacy RB, Lee HKH (2009) Adaptive design and analysis of supercomputer experiments. Technometrics 51(2):130–145

    Article  MathSciNet  Google Scholar 

  13. Gramacy RB, Lee HK (2012) Cases for the nugget in modeling computer experiments. Stat Comput 22(3):713–722

    Article  MathSciNet  Google Scholar 

  14. Gramacy RB, Taddy MA (2016) tgp: Bayesian treed Gaussian process models. R package version 2-4-14

  15. Greenberg D (1979) A numerical model investigation of tidal phenomena in the Bay of Fundy and Gulf of Maine. Mar Geod 2:161–187

    Article  Google Scholar 

  16. Gu M, Palomo J, Berger J (2018) RobustGaSP: robust Gaussian stochastic process emulation. R package version 0.5.6

  17. 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

    Article  Google Scholar 

  18. Johnson M, Moore L, Ylvisaker D (1990) Minimax and maximin distance design. J Stat Plan Inference 26:131–148

    Article  MathSciNet  Google Scholar 

  19. Joseph VR, Dasgupta T, Tuo R, Wu CJ (2015) Sequential exploration of complex surfaces using minimum energy designs. Technometrics 57:64–74

    Article  MathSciNet  Google Scholar 

  20. Joseph VR, Gul E, Ba S (2015) Maximum projection designs for computer experiments. Biometrika 102:371–380

    Article  MathSciNet  Google Scholar 

  21. Joseph VR, Hung Y (2008) Orthogonal-maximin Latin hypercube designs. Stat Sin 18:171–186

    MathSciNet  MATH  Google Scholar 

  22. Lam CQ, Notz WI (2008) Sequential adaptive designs in computer experiments for response surface model fit. Stat Appl 6:207–233

    Google Scholar 

  23. 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

    Article  MathSciNet  Google Scholar 

  24. Loeppky JL, Moore LM, Williams BJ (2010) Batch sequential designs for computer experiments. J Stat Plan Inference 140(6):1452–1464

    Article  MathSciNet  Google Scholar 

  25. Loeppky JL, Sacks J, Welch WJ (2009) Choosing the sample size of a computer experiment: a practical guide. Technometrics 51(4):366–376

    Article  MathSciNet  Google Scholar 

  26. 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

  27. 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

    MathSciNet  MATH  Google Scholar 

  28. Morris MD, Mitchell TJ (1995) Exploratory designs for computer experiments. J Stat Plan Inference 43:381–402

    Article  Google Scholar 

  29. Owen AB (1992) Orthogonal arrays for computer experiments, integration and visualization. Stat Sin 2:439–452

    MathSciNet  MATH  Google Scholar 

  30. 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

    Article  Google Scholar 

  31. Park JS (1991) Tuning complex computer codes to data and optimal designs. Unpublished Ph.D. thesis, University of Illinois, Champaign-Urbana

  32. Roustant O, Ginsbourger D, Deville Y (2018) DiceKriging: kriging methods for computer experiments. R package version 1.5.6

  33. 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)

    Article  MathSciNet  Google Scholar 

  34. Rasmussen CE, Williams CKI (2006) Gaussian processes for machine learning. The MIT Press, Cambridge, MA

    MATH  Google Scholar 

  35. Santner TJ, Williams BJ, Notz WI (2003) The design and analysis of computer experiments. Springer, New York

    Book  Google Scholar 

  36. Sacks J, Schiller SB, Welch WJ (1989) Designs for computer experiments. Technometrics 31:41–47

    Article  MathSciNet  Google Scholar 

  37. 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

    Chapter  Google Scholar 

  38. Sacks J, Welch WJ, Mitchell TJ, Wynn HP (1989) Design and analysis of computer experiments. Stat Sci 4:409–423

    Article  MathSciNet  Google Scholar 

  39. Tang B (1993) Orthogonal array-based Latin hypercubes. J Am Stat Assoc 88:1392–1397

    Article  MathSciNet  Google Scholar 

  40. 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

    Article  Google Scholar 

Download references

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

Authors

Corresponding author

Correspondence to C. Devon Lin.

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Published:

  • DOI: https://doi.org/10.1007/s42519-019-0077-0

Keywords

Navigation