CGP-package {CGP}R Documentation

The composite Gaussian process model package


Build nonstationary surrogate models for emulating computationally-expensive computer simulations (computer models).


Package: CGP
Type: Package
Version: 2.1-1
Date: 2018-06-11
License: LGPL-2.1

This package contains functions for fitting the composite Gaussian process (CGP) model, which consists of two Gaussian processes (GPs). The first GP captures the smooth global trend and the second one models local details. The model also incorporates a flexible variance model, which makes it more capable of approximating surfaces with varying volatility. It can be used as an emulator (surrogate model) for approximating computationally expensive functions that are not second-order stationary. When the underlying surface is stationary, the fitted CGP model should degenerate to a standard (stationary) GP model (\hat{\lambda}\approx 0).

The package implements maximum likelihood method to estimate model parameters and also provides predictions (with 5% and 95% prediction quantiles) for unobserved input locations. Leave-one-out cross validation diagnostic plot is also supported.

Gaussian correlation functions

g(\mathbf{h})=\exp(-\sum_{j=1}^p \theta_j h_j^2), \qquad l(\mathbf{h})=\exp(-\sum_{j=1}^p \alpha_j h_j^2)

(with unknown parameters \mathbf{\theta} and \mathbf{\alpha}) are used to describe the correlations in the global and local processes respectively.

For a complete list of functions, please use help(package="CGP"). Important functions are CGP, print.CGP, summary.CGP, predict.CGP and plotCGP.


Shan Ba and V. Roshan Joseph

Maintainer: Shan Ba <>


Ba, S. and V. Roshan Joseph (2012). Composite Gaussian Process Models for Emulating Expensive Functions. Annals of Applied Statistics, 6, 1838-1860.

[Package CGP version 2.1-1 Index]