ev_return_period: Estimate the return period for given quantiles

Description Usage Arguments Value See Also Examples

View source: R/5extreme.R

Description

For discharges of interrest,

Usage

1

Arguments

x

numeric vector containting the quantiles

fit

object of class evfit describing the underlying distribution, possibly created with evfit().

Value

a numeric vector of return periods.

See Also

evfit

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
data("ngaruroro")
ng <- as.xts(ngaruroro)

# yearly minima
minima <- apply.yearly(ng$discharge, min, na.rm = TRUE)

# fit a Weibull distribution
fit <- evfit(x = as.vector(minima), distribution = "wei")

# compute return periods
minima$rp <- round(ev_return_period(minima, fit), 2)

print(minima)
plot(discharge ~ rp, data = minima,
     xlab = "Flow in m^3/s", ylab = "Return period in years")

lfstat documentation built on May 30, 2017, 1:46 a.m.