RMwave: Wave Covariance Model / Cardinal Sine

Description Usage Arguments Details Value References See Also Examples

View source: R/RMmodels.R

Description

RMwave is a stationary isotropic covariance model, which is valid only for dimensions d ≤ 3. The corresponding covariance function only depends on the distance r ≥ 0 between two points and is given by

C(r)=sin(r)/r 1_{r>0} + 1_{r=0} .

It is a special case of RMbessel.

Usage

1
2

Arguments

var,scale,Aniso,proj

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

Details

The model is only valid for dimensions d ≤ 3. It is a special case of RMbessel for ν = 0.5.

This covariance models a hole effect (cf. Chiles, J.-P. and Delfiner, P. (1999), p. 92).

Value

RMwave returns an object of class RMmodel.

References

See Also

RMbessel, RMmodel, RFsimulate, RFfit.

Examples

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

model <- RMwave(scale=0.1)
x <- seq(0, 10, 0.02)
plot(model)
plot(RFsimulate(model, x=x))

Example output

Loading required package: sp
Loading required package: RandomFieldsUtils

Attaching package:RandomFieldsThe following object is masked frompackage:RandomFieldsUtils:

    RFoptions

NULL
NOTE: simulation is performed with fixed random seed 0.
Set 'RFoptions(seed=NA)' to make the seed arbitrary.
New output format of RFsimulate: S4 object of class 'RFsp';
for a bare, but faster array format use 'RFoptions(spConform=FALSE)'.

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