Stochastic Local Interaction Model: Geostatistics without Kriging
  Classical geostatistical methods face serious computational challenges if
they are confronted with large sets of spatially distributed data. We present a
simplified stochastic local interaction (SLI) model for computationally
efficient spatial prediction that can handle large data. The SLI method
constructs a spatial interaction matrix (precision matrix) that accounts for
the data values, their locations, and the sampling density variations without
user input. We show that this precision matrix is strictly positive definite.
The SLI approach does not require matrix inversion for parameter estimation,
spatial prediction, and uncertainty estimation, leading to computational
procedures that are significantly less intensive computationally than kriging.
The precision matrix involves compact kernel functions (spherical, quadratic,
etc.) which enable the application of sparse matrix methods, thus improving
computational time and memory requirements. We investigate the proposed SLI
method with a data set that includes approximately 11500 drill-hole data of
coal thickness from Campbell County (Wyoming, USA). We also compare SLI with
ordinary kriging (OK) in terms of estimation performance, using cross
validation analysis, and computational time. According to the validation
measures used, SLI performs slightly better in estimating seam thickness than
OK in terms of cross-validation measures. In terms of computation time, SLI
prediction is 3 to 25 times (depending on the size of the kriging neighborhood)
faster than OK for the same grid size.
