Learn R Programming

RandomFields (version 3.1.16)

RMbcw: Model bridging stationary and intrinsically stationary processes

Description

RMbcw is a variogram model that bridges between some intrinsically stationary isotropic processes and some stationary ones. It reunifies the RMgenfbm, RMdewijsian and RMgencauchy.

The corresponding centered semi-variogram only depends on the distance $r \ge 0$ between two points and is given by $$\gamma(r) = \frac{(r^{\alpha}+1)^{\beta/alpha}-1}{2^{\beta/alpha} -1}$$ where $0 < \alpha \le 2$ and $\beta

Usage

RMbcw(alpha, beta, var, scale, Aniso, proj)

Arguments

alpha
a numerical value; should be in the interval (0,2].
beta
a numerical value; should be in the interval (-infty,2].
var,scale,Aniso,proj
optional arguments; same meaning for any RMmodel. If not passed, the above variogram remains unmodified.

Value

RMbcw returns an object of class RMmodel

Details

For $betaa >0$, $beta<0$, $beta="0$" we="" have="" the="" generalised="" fractal="" brownian="" motion="" RMgenfbm, the generalised Cauchy model RMgencauchy, and the de Wisjian model RMdewijsian, respectively.

Hence its two arguments alpha and beta allow for modelling the smoothness and a wide range of tail behaviour, respectively.

References

  • Schlather, M (2014) A parameteric variogram model bridging between stationary and intrinsically stationary processes. arxiv 1412.1914.

See Also

RMgenfbm, RMgencauchy, RMdewijsian, RMmodel, RFsimulate, RFfit.

Examples

Run this code
RFoptions(seed=0) ## *ANY* simulation will have the random seed 0; set
##                   RFoptions(seed=NA) to make them all random again

model <- RMbcw(alpha=1, beta=0.5)
x <- seq(0, 10, 0.02)
plot(model)
plot(RFsimulate(model, x=x))

Run the code above in your browser using DataLab