rlibkriging (version 0.9-1)

simulate,KM-method: Simulation from a KM Object

Description

The simulate method is used to simulate paths from the kriging model described in object.

Usage

# S4 method for KM
simulate(
  object,
  nsim = 1,
  seed = NULL,
  newdata,
  cond = TRUE,
  nugget.sim = 0,
  checkNames = FALSE,
  ...
)

Value

A numeric matrix with nrow(newdata) rows and

nsim columns containing as its columns the simulated paths at the input points given in newdata.

XXX method simulate KM

Arguments

object

A KM object.

nsim

Integer: number of response vectors to simulate.

seed

Random seed.

newdata

Numeric matrix with it rows giving the points where the simulation is to be performed.

cond

Logical telling wether the simulation is conditional or not. Only TRUE is accepted for now.

nugget.sim

Numeric. A postive nugget effect used to avoid numerical instability.

checkNames

Check consistency between the design data X within object and newdata. The default is FALSE. XXXY Not used!!!

...

Ignored.

Author

Yann Richet

Details

Without a dedicated simulate method for the class "KM", this method would have been inherited from the "km" class. The dedicated method is expected to run faster. A comparison can be made by coercing a KM object to a km object with as.km before calling simulate.

Examples

f <-  function(x) 1 - 1 / 2 * (sin(12 * x) / (1 + x) + 2 * cos(7 * x) * x^5 + 0.7)
plot(f)
set.seed(123)
X <- as.matrix(runif(5))
y <- f(X)
points(X, y, col = 'blue')
k <- KM(design = X, response = y, covtype = "gauss")
x <- seq(from = 0, to = 1, length.out = 101)
s_x <- simulate(k, nsim = 3, newdata = x)
lines(x, s_x[ , 1], col = 'blue')
lines(x, s_x[ , 2], col = 'blue')
lines(x, s_x[ , 3], col = 'blue')