Calculate Return Levels

Calculates return levels for specified return periods from a fitted extreme value model.

Usage

calculate_return_levels(
  fit,
  return_periods = c(10, 50, 100),
  conf_level = 0.95
)

Arguments

fit

Result from fit_gev_annual_maxima or fit_gpd_threshold

return_periods

Numeric vector of return periods in years (default: c(10, 50, 100))

conf_level

Confidence level for intervals (default: 0.95)

Value

Data frame with:

• return_period: return period in years

• return_level: estimated return level

• lower: lower confidence bound

• upper: upper confidence bound

Examples

if (FALSE) { # \dontrun{
gev_result <- fit_gev_annual_maxima(data)
levels <- calculate_return_levels(gev_result, c(10, 50, 100))
print(levels)
} # }