vgmST: Constructing a spatio-temporal variogram

Description

Constructs a spatio-temporal variogram of a given type checking for a minimal set of parameters.

Usage

vgmST(stModel, ..., space, time, joint, sill, nugget, stAni)

Arguments

stModel

A string indentifying the spatio-temporal variogram model.

...

unused, but ensure an exact match of the following parameters.

space

A spatial variogram.

time

A temporal variogram.

joint

A joint spatio-temporal variogram.

sill

A joint spatio-temporal sill.

nugget

A joint spatio-temporal nugget.

stAni

A spatio-temporal anisotropy; the number of space units equivalent to one time unit.

Value

Returns an S3 object of class StVariogramModel.

Details

The different implemented spatio-temporal variogram models have the follwoing required parameters (see as well the example section) [object Object],[object Object],[object Object],[object Object],[object Object]

Examples

Run this code

# separable model: spatial and temporal sill will be ignored
# and kept constant at 1-nugget respectively. A joint sill is used.
separableModel <- vgmST("separable", 
                        space=vgm(0.9,"Exp", 147, 0.1),
                        time =vgm(0.9,"Exp", 3.5, 0.1),
                        sill=40)

# product sum model: spatial and temporal nugget will be ignored and kept
# constant at 0. Only a joint nugget is used.
prodSumModel <- vgmST("productSum",
                      space=vgm(39, "Sph", 343, 0),
                      time= vgm(36, "Exp",   3, 0), 
                      sill=41, nugget=17)

# sum metric model: spatial, temporal and joint nugget will be estimated
sumMetricModel <- vgmST("sumMetric",
                        space=vgm( 6.9, "Lin", 200, 3.0),
                        time =vgm(10.3, "Lin",  15, 3.6),
                        joint=vgm(37.2, "Exp",  84,11.7),
                        stAni=77.7)
                       
# simplified sumMetric model, only a overall nugget is fitted. The spatial, 
# temporal and jont nuggets are set to 0.
simpleSumMetricModel <- vgmST("simpleSumMetric",
                              space=vgm(20,"Lin", 150, 0),
                              time =vgm(20,"Lin", 10,  0),
                              joint=vgm(20,"Exp", 150, 0),
                              nugget=1, stAni=15)

# metric model
metricModel <- vgmST("metric",
                     joint=vgm(60, "Exp", 150, 10),
                     stAni=60)