RandomFields (version 3.0.10)

ExtremalGaussian: Extremal Gaussian process

Description

RPschlather defines an extremal Gaussian process.

Usage

RPschlather(phi, tcf, xi, mu, s)

Arguments

phi

an RMmodel, see Details.

tcf

an RMmodel specifying the extremal correlation function; either phi or tcf must be given.

xi,mu,s

the extreme value index, the location parameter and the scale parameter, respectively, of the generalized extreme value distribution. See Details.

Details

The parameter xi is always a number, i.e. $\xi$ is constant in space. In contrast, $\mu$ and $s$ might be constant numerical value or given a RMmodel, in particular by a RMtrend model.

The parameter phi can be any random field for which the expectation of the positive part is known at the origin.

It simulates Extremal Gaussian process $Z$ (also called Schlather model), which is defined by $$Z(x) = \max_{i=1}^\infty X_i \max(0, Y_i(x)),$$ where the $X_i$ are the points of a Poisson point process on the positive real half-axis with intensity $c x^{-2} dx$, $Y_i \sim Y$ are iid stationary Gaussian processes with a covariance function given by model, and $c$ is chosen such that $Z$ has standard Frechet margins. model must represent a stationary covariance model.

See Also

RMmodel, RPgauss, maxstable, maxstableAdvanced

Examples

RFoptions(seed=0)
model <- RMgauss()
x <- seq(0, 10, if (interactive()) 0.1 else 1)
z <- RFsimulate(RPschlather(model, xi=0, s=0.1), x)
plot(z)

if (.C("isAuthor", a=integer(1))$a) {

model <- RMgauss()
x <- seq(0,10, 0.02)
z <- RFsimulate(RPschlather(model, xi=0), x,
                n=if (interactive()) 100 else 1)
plot(z)
hist(unlist(z@data), 50, freq=FALSE)
curve(exp(-x) * exp(-exp(-x)), from=-3, to=8, add=TRUE)

## for some more sophisticated models see 'maxstableAdvanced'

}
RFoptions(seed=NA)