return.levels.mev: Return Levels for the MEVD/SMEV/TMEV extreme value...

In mevr: Fitting the Metastatistical Extreme Value Distribution MEVD

Return Levels for the MEVD/SMEV/TMEV extreme value distributions

Description

Calculate return levels for a MEVD, SMEV or TMEV extreme value distributions from an object of class mevr.

Usage

return.levels.mev(
  x,
  return.periods = c(2, 10, 20, 30, 50, 75, 100, 150, 200),
  ci = FALSE,
  alpha = 0.05,
  method = "boot",
  R = 502,
  ncores = 2L
)

Arguments

x	An object of class mevr, either fitted with the MEVD, SMEV or TMEV
return.periods	A vector of return periods in years, excluding 1.
ci	If ci=TRUE, confidence intervals are calculated depending on the type of distribution (only for MEVD or SMEV).
alpha	Number between zero and one giving the 1 - alpha confidence level. Defaults to alpha=0.05.
method	Character string giving the method for confidence interval calculation. Option method='boot' employs a parametric bootstrap that simulates data from the fitted model, and then fits the chosen MEVD type to each simulated data set to obtain a sample of parameters or return levels (very slow).
R	The number of bootstrap iterations.
ncores	Number of cores used for parallel computing of confidence intervals. Defaults to 2.

Details

Note that bootstraping the confidence intervals is very slow.

Value

A list with return levels, chosen return periods and, if ci=TRUE, alpha/2 and 1 - alpha/2 confidence intervals.

Examples

data(dailyrainfall)

fit <- fmev(dailyrainfall)
return.levels.mev(fit)
plot(fit)
 

mevr documentation built on June 30, 2024, 9:06 a.m.