R: On some covariance models based on normal scale mixtures

S10 {RandomFields} | R Documentation

## On some covariance models based on normal scale mixtures

### Description

Here the code of the paper on ‘On some covariance models based on normal scale mixtures’ is given.

### Author(s)

Martin Schlather, schlather@math.uni-mannheim.de http://ms.math.uni-mannheim.de/de/publications/software;

### References

- Schlather, M. (2010) On some covariance models based on normal scale mixtures. Bernoulli, 16, 780-797.

### Examples

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

    ### Example 10 in Schlather (2010).
    ## The field below has more than 80 million points. So the simulation
    ## takes a while
    y <- x <- seq(0, 10, len=256)
    T <- c(0, 0.02, 1275)
    col <- c(topo.colors(300)[1:100], cm.colors(300)[c((1:50) * 2,
             101:150)])

    model <- RMcoxisham(mu=c(1, 1), D=matrix(nr=2, c(1, 0.5, 0.5, 1)),
                        RMwhittle(nu=1))
    z <- RFsimulate(model, x, y, T=T, sp_lines=1500, every=10)
    plot(z, MARGIN.slices=3, col=col)
    plot(z, MARGIN.movie=3) # add 'file="ci.avi"' to get it stored

[Package RandomFields version 3.1.50 Index]