ramps (version 0.6-12)

corRExp2: Non-Separable Exponential Spatio-Temporal Correlation Structure

Description

This function is a constructor for the 'corRExp2' class, representing a non-separable spatial correlation structure. Letting $rs$ denote the spatial range, $rt$ the temporal range, and $lambda$ the space-time interaction, the correlation between two observations a distance $d$ apart in space and $t$ in time is $\exp(-d/rs - t/rt - \lambda (d/rs) (t/rt))$.

Usage

corRExp2(value = numeric(0), form = ~ 1,
            metric = c("euclidean", "maximum", "manhattan", "haversine"),
            radius = 3956)

Arguments

value

optional numeric vector of three parameter values for the exponential correlation structure, corresponding to the spatial range, temporal range, and space-time interaction. The range parameter values mus

form

one-sided formula of the form ~ S1+...+Sp+T, specifying spatial covariates S1 through Sp and the times T at which measurement were taken.

metric

optional character string specifying the distance metric to be used. The currently available options are "euclidean" for the root sum-of-squares of distances; "maximum" for the maximum difference; "manhattan" for t

radius

radius to be used in the haversine formula for great-circle distance. Defaults to the Earth's radius of 3,956 miles.

Value

• Object of class 'corRExp2', inheriting from class 'corSpatioTemporal', representing a non-separable spatial correlation structure.

References

Cressie, N. and Huang, H.-C. (1993) Classes of Nonseperable, Spatio-Temporal Stationary Covariance Functions, Journal of the American Statistical Association, 94, 1330-1340.

See Also

corClasses, Initialize.corStruct, summary.corStruct

Examples

sp1 <- corRExp2(form = ~ x + y + t)

spatDat <- data.frame(x = (0:4)/4, y = (0:4)/4, t=(0:4)/4)

cs1Exp <- corRExp2(c(1, 1, 1), form = ~ x + y + t)
cs1Exp <- Initialize(cs1Exp, spatDat)
corMatrix(cs1Exp)

cs2Exp <- corRExp2(c(1, 1, 1), form = ~ x + y + t, metric = "man")
cs2Exp <- Initialize(cs2Exp, spatDat)
corMatrix(cs2Exp)