constrainedKriging-package: constrainedKriging a package for nonlinear spatial predictions with local change of support

Description

The constrainedKriging package provides functions for tow-dimensional spatial interpolation by constrained, covariance-matching constrained and universal (external drift) kriging for points or block of any shape for data with a nonstationary mean function and an isotropic weakly stationary variogram. The linear spatial interpolation methods, constrained and covariance-matching constrained kriging, provide approximately unbiased prediction for nonlinear target values under change of support.

In principle, the package provides two user functions, preCKrige and CKrige to calculate spatial prediction in two steps:

1. Call of preCKrige to calculate the variance-covariance matrices for defined sets of points or polygons (blocks). 2. Call of CKrige by using the output of preCKrige to calculate the spatial interpolation by one of the three kriging methods.

Arguments

eqn

$\boldsymbol{\Sigma}$

deqn

$$MSPE[ \widehat{\mathbf{Y}}_{CMCK} ] = MSPE[ \widehat{\mathbf{Y}}_{UK} ] + (\mathbf{P}_{1}-\mathbf{Q}_{1})(\mathbf{P}_{1}-\mathbf{Q}_{1}).$$

itemize

MSPE of the universal kriging predictor for 1 block:

item

MSPE of the constrained kriging predictor for 1 block:
MSPE of the covariance-matching constrained kriging predictor for m blocks:

Details

The constrained kriging predictor proposed by Cressie (1993) and the covariance-matching constrained kriging predictor proposed by Aldworth and Cressie (2003) are linear in the data like the universal kriging predictor. However, the constrained kriging predictor satisfies in addition to the unbiasedness constraint of universal kriging a second constraint that matches the variances of the predictions to the variances of the prediction target (either a point value or a block mean).The covariance-matching constrained kriging predictor matches for a set of blocks both the variances and covariances of predictions and prediction targets ( points or block means) and is the extended version of the constrained kriging predictor.Like constrained kriging, covariance-matching constrained kriging is less biased than universal kriging for nonlinear predictions and exactly unbiased if the spatial variable is Gaussian.

The formulas of the three kriging predictors are given below:

universal kriging predictor for 1 block:$$\widehat{Y}_{\mathrm{UK}}(B) = \mathbf{x}(B)^{\prime}\widehat{\beta}_{\mathrm{GLS}} + \mathbf{ c }^{\prime} \boldsymbol{\Sigma} ( \mathbf{ Z } - \mathbf{X}\widehat{\beta}_{\mathrm{GLS}} ) ,$$

constrained kriging predictor for 1 block: $$\widehat{Y}_{\mathrm{CK}}(B) = \mathbf{x}(B)^{\prime}\widehat{\beta}_{\mathrm{GLS}} + K \mathbf{ c }^{\prime} \boldsymbol{\Sigma} ( \mathbf{ Z } - \mathbf{X}\widehat{\beta}_{\mathrm{GLS}}) ,$$ covariance-matching constrained kriging predictor for m blocks: $$\widehat{ \mathbf{ Y } }_{\mathrm{CMCK}} = \mathbf{X}_{m}\widehat{\beta}_{\mathrm{GLS}} + \mathbf{ K }^{\prime} \mathbf{C}^{\prime}\boldsymbol{\Sigma} (\mathbf{ Z } - \mathbf{X}\widehat{\beta}_{\mathrm{GLS}}),$$

References

Aldworth, J. and Cressie, N. (2003). Prediction of nonlinear spatial functionals. Journal of Statistical Planning and Inference, 112, 3--41

Cressie, N. (1993). Aggregation in geostatistical problems. In A. Soares, editor, Geostatistics Troia 92, 1, pages 25--36, Dordrecht. Kluwer Academic Publishers.

Hofer, C. and Papritz, A. (2010). Predicting threshold exceedance by local block means in soil pollution surveys. Mathematical Geosciences. 42, 631-656, doi: 10.1007/s11004-010-9287-4

Hofer, C. and Papritz, A. (in preparation). constrainedKriging: an R-package for customary, constrained and covariance-matching constrained point or block kriging. Computers & Geosciences.