RobustGaSP: Robust Gaussian Stochastic Process Emulation in R

Gaussian stochastic process (GaSP) emulation is a powerful tool for approximating compu tationally intensive computer models. However, estimation of parameters in the GaSP emulator is a challenging task. No closed-form estimator is available and many numerical problems arise with standard estimates, e.g., the maximum likelihood estimator. In this package, we implement a marginal posterior mode estimator, for special priors and parameterizations. This estimation method that meets the robust parameter estimation criteria was discussed in Gu et al. (2018); mathematical reasons are provided therein to explain why robust parameter estimation can greatly improve predictive performance of the emulator. In addition, inert inputs (inputs that almost have no effect on the variability of a function) can be identified from the marginal posterior mode estimation at no extra computational cost. The package also implements the parallel partial Gaussian stochastic process (PP GaSP) emulator (Gu and Berger (2016)) for the scenario where the computer model has multiple outputs on, for example, spatial-temporal coordinates. The package can be operated in a default mode, but also allows numerous user specifications, such as the capability of specifying trend functions and noise terms. Examples are studied herein to highlight the performance of the package in terms of out-of-sample prediction.

Mengyang Gu , Jesus Palomo , James O. Berger

Supplementary materials

Supplementary materials are available in addition to this article. It can be downloaded at

CRAN packages used

DiceKriging, GPfit, mleGP, spatial, fields, RobustGaSP, lhs, sensitivity, nloptr

CRAN Task Views implied by cited packages

ExperimentalDesign, Spatial, Distributions, Environmetrics, Optimization, SocialSciences


Text and figures are licensed under Creative Commons Attribution CC BY 4.0. The figures that have been reused from other sources don't fall under this license and can be recognized by a note in their caption: "Figure from ...".


For attribution, please cite this work as

Gu, et al., "The R Journal: RobustGaSP: Robust Gaussian Stochastic Process Emulation in R", The R Journal, 2019

BibTeX citation

  author = {Gu, Mengyang and Palomo, Jesus and Berger, James O.},
  title = {The R Journal: RobustGaSP: Robust Gaussian Stochastic Process Emulation in R},
  journal = {The R Journal},
  year = {2019},
  note = {},
  doi = {10.32614/RJ-2019-011},
  volume = {11},
  issue = {1},
  issn = {2073-4859},
  pages = {112-136}