Kriging Models

 

Kriging is a geostatistical method of spatial data interpolation. The mathematical model of kriging is named after D.G. Krige, who first introduced a version of this spatial prediction process. Kriging has been extensively described in the literature since Sacks et al. proposed the application of kriging in computer experiments.

Unlike the real tests, computer analysis codes are deterministic therefore it is not influenced with measurement errors. Hence, the approximate models can be defined as a combination of a regression model plus a departure term:

 

where,  is the approximate model,  is a polynomial type regression model, and  is a Gaussian random process with. If the regression model  globally approximates the design space, the departure term  represents the localized deviations so that the Kriging model interpolates the  sampled points.

The covariance matrix of  is given by

 

 

where,  is the correlation matrix and  is the correlation function between any two of the  sampled points. Hence,  is a  symmetric matrix with ones in the diagonal term. There are many correlation functions . Among them, the Gaussian type is widely used

 

 

where,  are the unknown correlation parameters to fit model. The estimates, of the response  at the untried values of  are given by

 

 

The correlation vector between  and the sampled points  is given by:

 

 

In the estimates, the unknown coefficients of regression model is determined as

 

 

Also, in order to determine the unknown correlation parameters , the estimate of the variance  (not the variance in the observed data)

 

 

is introduced. Hence, the correction parameters  is determined by solving

 

 or

 

While any values for  create an interpolation model, the best kriging model is found by solving the k-dimensional unconstrained optimization problems described in the above.

 

Reference

1.  Matheron G. Principles of geostatistics, Economic Geology 1963; 58:1246-1266.

2.  Sacks J, Welch WJ, Mitchell TJ, Wynn HP, Design and analysis of computer experiments. Statistical Science 1989; 4:409-435.

3.  Sacks J, Susannah SB, Welch WJ. Design for computer experiments. Technometrics 1989; 31:41-47