RMsinepower
is a isotropic covariance model. The
corresponding covariance function, the sine power function of
Soubeyrand, Enjalbert and Sache, only depends on the angle $0 \le \theta \le \pi$ between two points on the sphere and is given by
RMsinepower(alpha, var, scale, Aniso, proj)
RMmodel
. If not passed, the above
covariance function remains unmodified.RMsinepower
returns an object of class RMmodel
RMchoquet
.
RMmodel
,
RFsimulate
,
RFfit
,
spherical models
,
RMchoquet
RFoptions(seed=0) ## *ANY* simulation will have the random seed 0; set
## RFoptions(seed=NA) to make them all random again
RFoptions(coord_system="sphere")
model <- RMsinepower(alpha=1.7)
plot(model, dim=2)
## the following two pictures are the same
x <- seq(0, 0.4, 0.01)
z1 <- RFsimulate(model, x=x, y=x)
plot(z1)
x2 <- x * 180 / pi
z2 <- RFsimulate(model, x=x2, y=x2, coord_system="earth")
plot(z2)
stopifnot(all.equal(as.array(z1), as.array(z2)))
RFoptions(coord_system="auto")
Run the code above in your browser using DataLab