update.Kriging: Update a 'Kriging' Object with New Points

View source: R/KrigingClass.R

update.KrigingR Documentation

Update a Kriging Object with New Points

Description

Update a Kriging model object with new points

Usage

## S3 method for class 'Kriging'
update(object, newy, newX, ...)

Arguments

object

S3 Kriging object.

newy

Numeric vector of new responses (output).

newX

Numeric matrix of new input points.

...

Ignored.

Value

No return value. Kriging object argument is modified.

Caution

The method does not return the updated object, but instead changes the content of object. This behaviour is quite unusual in R and differs from the behaviour of the methods update.km in DiceKriging and update,KM-method.

Author(s)

Yann Richet yann.richet@irsn.fr

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(10))
y <- f(X)
points(X, y, col = "blue")

k <- Kriging(y, X, "matern3_2")

x <- seq(from = 0, to = 1, length.out = 101)
p <- predict(k, x)
lines(x, p$mean, col = "blue")
polygon(c(x, rev(x)), c(p$mean - 2 * p$stdev, rev(p$mean + 2 * p$stdev)),
 border = NA, col = rgb(0, 0, 1, 0.2))

newX <- as.matrix(runif(3))
newy <- f(newX)
points(newX, newy, col = "red")

## change the content of the object 'k'
update(k, newy, newX)

x <- seq(from = 0, to = 1, length.out = 101)
p2 <- predict(k, x)
lines(x, p2$mean, col = "red")
polygon(c(x, rev(x)), c(p2$mean - 2 * p2$stdev, rev(p2$mean + 2 * p2$stdev)),
 border = NA, col = rgb(1, 0, 0, 0.2))

rlibkriging documentation built on Sept. 22, 2022, 5:07 p.m.