Questions tagged [kriging]

Kriging is spatial prediction based on a stochastic model of a spatial random field. Such models and methods can also be used in non-spatial context, and then often known as gaussian processes.

67 questions
Filter by
Sorted by
Tagged with
32 views

Can I replace the distribution in Gaussian Process Regression with a different regression?

Sorry for the confusing title. Let me try to clarify: I have a time series of wind speeds with some missing points here and there. I want to interpolate these points and have tried mainly polynomial ...
13 views

Iterative solution to estimate nugget effect

I have to implement the equation to fit an analytical curve to an empirical covariance function. Until now I used the equation $$c = ae^{-b \Delta t}$$ with $a$ and $b$ are the parameters I want ...
100 views

Is the spherical covariance function not positive definite for d > 3?

I read in a textbook (Japanese one) that the spherical covariance function is only valid for dimensions $d = 1,$ $2,$ and $3.$ I have the following questions: Does that mean the spherical covariance ...
15 views

Cokriging variances differ using cross validation

I'm investigating cokriging using various metals in the Meuse dataset but the variances output by R when I predict values at gridded points differ substantially from the variances produced by cross ...
13 views

Why do you need a variogram for Kriging? goldingn/gpe package?

I am using golingn/gpe (github) package, and it does not provide a variogram and instead look at co-variances. Is it possible to do kriging without providing variograms?
17 views

Kriging with Weighted Data

I have a point-level dataset of apartment building level rent. I have the average rent per square foot per building and the number of units in that building. I would like to krige a surface of points ...
57 views

How to handle multiple points at the same location in spatial interpolation?

I am new to the topic of spatial interpolation and would appreciate your opinion on a general question which has arisen. Suppose I have a data set containing rental rates for different apartments in ...
48 views

Poor performance from gstat::krige with a noise predictor

I'm new to kriging, and I'm considering replacing a use of inverse-distance weighting (IDW) in a spatial modeling project (implemented with gstat::idw in R) with a ...
34 views

How do nugget interactions work in Gaussian Processes/Kriging?

How do nuggets and nugget interactions fit into the variogram framework? I am especially interested in the case where there is more than one distance term being used (e.g. space and time): where you ...
53 views

Gaussian process with interval observations

The stochastic process $(X_t)_{t \in T}$ is a Gaussian process if the marginal distribution of $X_{t_1}, \ldots, X_{t_n}$ is a multivariate Gaussian distribution for all $t_1, \ldots, t_n \in T$. Let ...
41 views

Cokriging, zero distance semivariance [gstat]

Trying cokriging with simulated data, I faced a problem that did not seem one in the demo(cokriging) with the meuse dataset: I can't use ...
163 views

Correct way to compute the error variance of ordinary kriging

I'm learning ordinary kriging and I found some discrepancies in the method of computing error variance among different materials I read, in general there are 2 different formula: Method 1: ...
138 views

spacetime R: How to handle missing data in a space-time-full data structure for spatio-temporal kriging purposes?

I am using R and the spacetime package. I am having problems using STFDF. I want to use the STF data-structure since I have spacetime data with recurrent observations for fixed spatial coordinates. ...
27 views

Noisy conditional simulation

A conditional random field $Z_C(x)$ is a random field whose realisations $z_C(x)$ always take the same values $z_C(x_a)$ at locations $x_a$. Realisations of $Z_C(x)$ can be produced as follows (...
2k views

Ordinary kriging example step by step?

I have followed tutorials online for spatial kriging with both geoR and gstat (and also ...
4k views

What is the nugget effect?

I don't understand exactly what is meant by the term "nugget effect" in geostatistics. When looking at empirical variograms plotting the variogram $\gamma(h)$ vs. the lag $h$, the nugget is defined as ...
65 views

Are Kriging's residuals (i.e. $Z-\hat{Z}$) spatially independent?

Suppose data ${Z(s_i ):i=1, ..., n}$ are observed at spatial locations ${s_i :i=1, ..., n}$. To carry out the spatial prediction (predict un unknown $Z(s_0)$ at a known location $s_0$) we can use a ...
377 views

Is a function describable by a Gaussian process smooth?

I understand that a stochastic process or function is considered a Gaussian process if sampling from it at any point some set of times yields a set of observations that match a Gaussian random ...
115 views

How do you interpret this variogram?

Description: 8000 spatial data points spanned over an entire state 200 bins are used My question: Is the variogram telling something about the nature of the data? Why is it fluctuating? Should I do ...
34 views

What is the statistically correct way of performing the following interpolation?

I have a table with thousands of entries similar to the following: It is desired to determine Property 2 at the unknown locations. Property 1 and Property 2 are spatially correlated. By that I mean ...
385 views

Spherical vs. Exponential Kriging Covariance Functions

A statistical epidemiologist colleague of mine told me that in comparing spherical vs exponential kriging covariance functions, only the latter (i.e., exponential) function is generally a valid model. ...
530 views

How do I estimate the prediction interval of back transformed log-normal data from Gaussian process?

I have some data that are clearly positively skewed and follow a log-normal distribution, lets assume the initial data is $Z = exp(Y)$, where $Y \sim N(\mu,\sigma^2)$. A Gaussian process assumes ...
307 views

What is the difference between accounting for anisotropy and trend removal when performing Kriging?

Without being geostatistician, I read a bit about anisotropy detection, mostly from ArcGIS documentation and the R gstat package tutorial. But still, it is hard to have a confident understanding of ...
80 views

What is 5th order kriging?

For our foray into geostatistics we got data that consists of measurements taken from the soil. The dataset has like concentrations of various different minerals. We were divided into a number of ...
213 views

Kriging variance results

I'm quite a newb at statistics and interpolation, and I cannot understand how to interpret the error estimation computed by Kriging. For example, I performed kriging on temperature values (Celsius ...
204 views

Determining covariance of irregularly spaced spatial data

I'm comparing concentration $C$ of a contaminant in the same spatial region at two time point 2000 and 2010 with sample size of $N_{2000}$ = 51 and $N_{2010}$ = 26 (not all the samples are from the ...
74 views

Explain Like I'm Five version for Variograms in R's gstat package

Long story short, I asked a question on StackOverflow about Variograms in the gstat package in R. The person who answered gave me some tips on creating the variogram using the package. My dataset is ...
2k views

what is the difference between Bayesian optimization and kriging?

Both methods use Gaussian process, and kriging uses the Best Linear Unbiased Predictor (BLUP) to predict the mean (this is not seen in Bayesian optimization?). At the bottom line, they also have ...
714 views

Is kriging suitable for high dimensional regression problems?

I would like to point out that I am new to this field, so if I am not clear please forgive me (and correct me). I set up a DoE (Design of Experiment) with 11 inputs and 121 runs. I used a STOA (...
106 views

Performing inference across multiple target variables with Gaussian Processes regression

Suppose I have a design matrix $\mathbf{X}$ with targets $\mathbf{y}_1$, when plotted looks like so: the data is very sparse, and the magnitude of the targets is shown by color bar, where each ...
408 views

How to implement a Gaussian Markov Random field (GMRF) model in R?

I would like to model a (conditional) GMRF using a linear mixed effects model without having grid Data but only a neighbourhood matrix $W$. My model is given by $$Y=X\beta+ \epsilon$$ and the error ...
72 views

Estimation of time for a specific value of a variable

I have a data set: ...
44 views

Doubts on the methodology used for build and evaluate a Kriging model

I have run 3000 finite element (deterministic) simulations of a physical model, which has 10 design variables and one output variable. I have generated the sampling plan using the LH method. My ...
58 views

Kriging with mean data instead of point data

Kriging is a technique to predict a realisation of a Gaussian process. The values of the realisation are known at a finite subset of points and we would like to optimally predict the values of the ...
346 views

Why is positive definiteness necessary for kriging?

I understand from wikipedia that a variogram model must be positive definite to be used for kriging: Note that the experimental variogram is an empirical ...
30 views

Is there a way to find optimal sensor node locations in a domain when all data is known?

I have an x-y domain where I know the snow depth everywhere. (i.e. the granularity of the values is such that I can assume I know it everywhere). I want to use this information to inform where to ...
81 views

Universal Kriging

In universal kriging the general model is $$y(x)=h(x)^{t}\boldsymbol\beta + f(x)+\epsilon(x)$$ where $h(x)^{t}$ is a regression function such as $\left [1,x,x^{2} \right ]$ and $\boldsymbol\beta$ are ...
158 views

Universal kriging with non-zero mean vector

Its is known that in ordinary Kriging (Gaussian Process) the mean and variance at any new point is given as \\\begin{pmatrix} \mathbf{y}\\y^{*} \end{pmatrix} \sim N(\mathbf{0},\begin{bmatrix} K &...
226 views

Kriging: fit.method and dX argument in variogram

I would like to krige residuals (from multiple linear regression) of yearly precipitation totals from a 50 years time series. Every year has been regressed individually. The residuals will be added to ...
855 views

Transformation of non-negative data including zeros

at first I want to mention that I am fully aware of this question here as well as the answers. Still, things won't work out as intended (using R). I have a lot of hourly rainfall data that include ...
814 views

Why we krige residuals in regression-kriging?

Wikipedia says: In applied statistics, regression-kriging (RK) is a spatial prediction technique that combines a regression of the dependent variable on auxiliary variables (such as parameters ...
51 views

Why likelihood based geostatistical modelling slower than non likelihood based counterpart

Likelihood based geo-statistics (geoR etc.) are usually slower than non-likelihood based geo-statistics (i.e. those based on just least square fitting, for example <...
254 views

Spatial prediction on surface: very fine grid vs coarse grid + quick interpolation

Once I have fitted a spatial model (point-referenced data), I need to make a prediction map. A natural approach is to make prediction over a fine grid over the region. However, the required ...
263 views

Is there a recursive version of Kriging or Inverse Distance spatial interpolation?

Classic use-case of Kriging: you have a 2d space, you have $n$ observations, each of them representing an exploratory dig. It has a $x$ and $y$ coordinate, and a $V$ representing the value discovered ...
383 views

Help understanding a kriging variation for bare earth extraction

Problem: The authors of a paper (http://www.isprs.org/proceedings/XXXIV/part3/papers/paper106.pdf) develop a bare earth extraction algorithm for LiDAR that is based on kriging. What I don't ...
443 views

Gstat: Modelled semivariogram values not matching plotted model using the variogramLine function

I am trying to extract the semivariance values associated with a given semivariogram model developed in gstat, the end goal being to compare modelled semivariance with observed semivariance at defined ...
337 views

Kriging with individual errors

When evaluating a Kriging prediction, it is possible to include a hyperparameter $\lambda$ to account for noise in the data. $\lambda$ can be estimated as a parameter in a maximum likelihood ...
130 views

How to implement Kriging without x and y coordinates

I have some sample measurements taken from different locations. I want to predict the measurements in some surrounding locations that were not studies, for this I want to use Kriging. However, My ...
According to Dubrule's Cross validation of kriging in a unique neighborhood, it is possible to compute leave-one-out the gaussian process prediction $\hat{Y}_{-i}(x_i)$ at a point $x_i$ from the ...