rlplot: Return Level Plot for a Fitted Object of an Extreme-Value Distribution

Description

Plots several return levels against the return period for a fitted object from one of the ismev functions: gev.fit and gpd.fit.

Usage

rlplot(z, ci = 0.05, add.ci = FALSE)

Arguments

A list object as returned by one of gev.fit or gpd.fit (with appropriate class attribute added).

The (1-ci)*100 confidence value.

add.ci

logical if true will add confidence bounds to plot.

Value

A plot is created. If assigned to an object, a list will be returned with the following items.
periodThe return periods used for calculating the return levels.
levelThe estimated return level for each return period.
lowerIf add.ci is TRUE, then this is a vector of lower 1-ci confidence bounds. Otherwise the value is NULL.
upperIf add.ci is TRUE, then this is a vector of upper 1-ci confidence bounds. Otherwise the value is NULL.

Details

Given a fitted list object from gev.fit or gpd.fit--attributed with the class gev.fit or gpd.fit, respectively--the return level plot is generated. Confidence bounds, if included, are found by the delta method, which is generally appropriate for shorter return periods, but not for longer return periods because the return level distribution is generally skewed. Therefore, if a plot with better estimates of the confidence bounds are desired, use add.ci=FALSE, and use the R function lines to add different bounds (e.g., using values obtained from the gev.parameterCI or gpd.parameterCI functions).

This function is simply a modification of ismev functions gpd.rl and gev.rl (Coles, 2001).

References

Coles, S. (2001) An introduction to statistical modeling of extreme values. London: Springer-Verlag, 208pp.

Examples

Run this code

data(ftcanmax)
fit <- gev.fit( ftcanmax[,"Prec"])
class( fit) <- "gev.fit"
rlplot( fit)

Run the code above in your browser using DataLab