GADGET (version 0.2.0)

bo09_toy: Two-Input Toy Model from Bastos and O'Hagan (2009)

Description

A toy model for the demonstration of Gaussian process validation plots and statistics. It is a product of a rational function in x[1] and an exponential function in x[2]. The function is smoother in x[2] than x[1]. The function's domain is the unit square.

Usage

bo09_toy(x)

Arguments

x

A vector of length 2.

Value

A scalar.

References

Bastos, L. S., & O'Hagan, A. (2009). Diagnostics for gaussian process emulators. Technometrics, 51(4), 425<U+2013>438, <doi:10.1198/TECH.2009.08019>.

Examples

# NOT RUN {
print(bo09_toy(c(0.5,0.25)))

n <- 100
a <- seq(0,1,length.out = n)
x <- expand.grid(a,a)
z <- matrix(apply(x,1,bo09_toy),n,n)
graphics::contour(a,a,z, xlab = expression(x[1]), ylab = expression(x[2]))

# }