Skip to contents

Two co-kriging routines live in the package and both have callers that spell the entry point morie_cokriging(): morie_cokrig (coords, z1, z2, s0, cross_vario, model: the ordinary co-kriging system with explicit variograms) and cokrg (x, y, coords, target, ...: the exponential model with fitted sills). The convention is picked from the first argument: a coordinate matrix routes to morie_cokrig, a numeric vector to cokrg. Call either routine directly to avoid the dispatch.

Usage

morie_cokriging(...)

Arguments

...

Arguments for morie_cokrig() or cokrg().

Value

Whatever the selected routine returns.

Examples

set.seed(1)
pts <- rbind(c(0, 0), c(1, 0), c(0, 1), c(1, 1))
morie_cokriging(pts, c(1, 2, 3, 4), rep(9, 4), c(0.5, 0.5),
                cross_vario = function(h) 0 * h)$prediction
#> [1] 2.5
morie_cokriging(rnorm(12), rnorm(12), matrix(runif(24), 12, 2),
                target = c(0.3, 0.7))$estimate
#> [1] 0.4127342