RandomFields (version 3.1.12)

S02: Models for stationary max-stable random fields

Description

Here the code of the paper on Models for stationary max-stable random fields is given.

References

• Schlather, M. (2002) Models for stationary max-stable random fields.Extremes5, 33-44.

Examples

RFoptions(seed=0, xi=1)
## seed = 0 : *ANY* simulation will have the random seed 0; set
##            RFoptions(seed=NA) to make them all random again
## xi = 0.5: Frechet margins with alpha=2
StartExample()
## Due to change in the handling the seeds here are different from the
## seed in the paper.

x <- seq(0, 10, length=128)
if(RFoptions()$internal$examples_red){warning("reduced <x>");x<-seq(0,10,length=4)}
# Fig. 1-4
\dontshow{plot(RFsimulate(RPsmith(RMgauss(s=1.5)), x, x))   # < 1 sec
plot(RFsimulate(RPsmith(RMball(s=RRspheric(2, 3,
R=3.3))), x, x)) # 30  sec
plot(RFsimulate(RPschlather(RMexp()), x, x))      #   1 sec
plot(RFsimulate(RPschlather(RMgauss()), x, x))    #  17 sec
}

FinalizeExample()