plot.maxstab: Model checking of a fitted max-stable model

Description Usage Arguments Details Value Author(s) Examples

Description

This function produces several plots to assess the goodness of fit of a fitted max-stable model.

Usage

1
2
## S3 method for class 'maxstab'
plot(x, ..., sites)

Arguments

x

An object of class maxstab. Most often, this will be the output of fitmaxstab or lsmaxstab.

...

Here for compatibility reasons but not yet implemented.

sites

A vector of integer of length 4 specifying the locations to be used for the model checking. If missing, locations will be choosen randomly.

Details

The diagonal plots are return level plots. The lower ones are qq-plots (on the Gumbel scale) between observed pairwise maxima for each block, e.g. year, and the ones obtained by simulations from the fitted model. The upper plot compares the fitted extremal coefficient functions to semi-empirical estimates from the F-madogram - see fmadogram. The two remaining plots are the stations locations and a qq-plot of blockwise maxima where the block size is 4.

Value

Several diagnostic plots.

Author(s)

Mathieu Ribatet

Examples

1
2
3
4
5
6
7
8
9
n.site <- 20
n.obs <- 50
coord <- matrix(runif(2 * n.site, 0, 10), ncol = 2)
colnames(coord) <- c("lon", "lat")
data <- rmaxstab(n.obs, coord, "powexp", nugget = 0, range = 3, smooth =
1)
fitted <- fitmaxstab(log(data), coord, "powexp", y ~ 1, y ~ 1, y ~ 1,
nugget = 0)
plot(fitted)


Search within the SpatialExtremes package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.