RMcurlfree: Curlfree Covariance Model

Description Usage Arguments Details Value References See Also Examples

View source: R/RMmodels.R

Description

RMcurlfree is a multivariate covariance model which depends on a univariate stationary covariance model where the covariance function phi(h) is twice differentiable.

The corresponding matrix-valued covariance function C of the model only depends on the difference h between two points and it is given by the following components:

Usage

1
RMcurlfree(phi, which, var, scale, Aniso, proj)

Arguments

phi

a univariate stationary covariance model (2- or 3-dimensional).

which

vector of integers. If not given all components are returned; otherwise the selected components are returned.

var,scale,Aniso,proj

optional arguments; same meaning for any RMmodel. If not passed, the above covariance function remains unmodified.

Details

The model returns the potential field in the first component, the corresponding curlfree field and field of sources and sinks in the last component.

See also the models RMdivfree and RMvector.

Value

RMcurlfree returns an object of class RMmodel.

References

See Also

RMderiv, RMdivfree, RMvector, RMmodel, RFsimulate, RFfit.

Examples

1
2
3
4
5
6
7
8
9
RFoptions(seed=0) ## *ANY* simulation will have the random seed 0; set
##                   RFoptions(seed=NA) to make them all random again

model <- RMcurlfree(RMgauss(), scale=4)
plot(model, dim=2)

x.seq <- y.seq <- seq(-10, 10, 0.2)
simulated <- RFsimulate(model=model, x=x.seq, y=y.seq)
plot(simulated, select.variables=list(1, c(1, 2:3), 4))

RandomFields documentation built on Jan. 19, 2022, 1:06 a.m.