Augmented universal kriging system

The system with matrix \((C, X; X', 0)\), unknowns \((\lambda, \mu)\) and right-hand side \((c_0, x_0)\), with the data covariance matrix \(C\), covariances \(c_0\) with the target and trend design \(X\) (default ordinary kriging); returns the matrix, right-hand side, weights and Lagrange multipliers (Cressie 1993, section 3.4).

Usage

KrigingSystem(coords, target, model, X = NULL, x0 = NULL)

Arguments

coords

Matrix of data locations.

target

Target location.

model

Covariance model (KrigingCovariance).

X

Trend design at the data.

x0

Trend design at the target.

Value

List with matrix, rhs, weights, lagrange.

Examples

KrigingSystem(rbind(c(0, 0), c(1, 0)), c(0.5, 0), list(model = "Exp", psill = 1, range = 1))$weights
#> [1] 0.5 0.5