spatio temporal anisotropy: help with stAni argument

Dear Bruno,

stAni relates spatial and temporal distances to each other in the model 
for the data set under study. A value of 100 km/day reflects that two 
stations 100 km apart have about the same strength of dependence as the 
values one day apart at the same location (This is in general only true 
for at most one pair of spatial and temporal distances and a 'global' 
compromise is sought). The function gstat::estiStAni provides a few 
heuristics to estimate the spatio-temporal anisotropy.

When a spatio-temporal vgm including a metric component is fitted, the 
stAni parameter is optimized under the objective to find the theoretical 
vgm surface that best fits the empirical one.

In krigeST, an additional and approximate stAni parameter can be 
provided to ease the (pre-)selection of the most relevant observations 
in the space-time cube for the prediction when local kriging is carried 
out. Note that the most relevant (i.e. strongest correlated) 
observations need not to form a cube, but depend on the actual vgm-model 
in use (see also the vignette "Spatio-Temporal Interpolation using 
gstat" at [1]).

The unit of stAni depends on the spatial and temporal units of your data 
set. Check the units of the empirical variogram to find out which units 
R/gstat uses for your data set.

Miss-specifying the stAni parameter will put too large or too small 
weights on the observations with a temporal difference resulting in bad 
predictions (this also holds vise-versa for spatial distances).

HTH,

  Ben

[1] 
https://cran.r-project.org/web/packages/gstat/vignettes/spatio-temporal-kriging.pdf

spatio temporal anisotropy: help with stAni argument

Thread (3 messages)