Models for stationary max-stable random fields
Here, the code of the paper on ‘Models for stationary max-stable random fields’ is given.
Martin Schlather, schlather@math.uni-mannheim.de, https://www.wim.uni-mannheim.de/schlather/
Schlather, M. (2002) Models for stationary max-stable random fields. Extremes 5, 33-44.
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
## Due to change in the handling the seeds here are different from the
## seeds in the paper.
x <- seq(0, 10, length=128)
# Fig. 1-4
## Not run: \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
}
## End(Not run)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.