confInt_sTL: Estimated confidence intervals using sTL

In flood: Statistical Methods for the (Regional) Analysis of Flood Frequency

Description

Estimated regional (or local) (1-alpha)-confidence intervals for an estimated quantile by using seasonal TL(0,1)-moments (sTL).

Usage

confInt_sTL(x1, x2, p, j = 1, alpha = 0.05)

Arguments

x1	vector or matrix of observations from season 1 (rows: observations, columns: stations).
x2	vector or matrix of observations from season 2 (rows: observations, columns: stations).
p	a probability.
j	quantile and parameter estimation for the jth station (jth column of x). Irrelevant if is x1 and x2 are vectors.
alpha	confidence level for confidence interval.

Value

List of

• ci confidence interval.

• quant estimated quantile from a two-component GEV using trimmed L-moments (leftrim=0, rightrim=1).

• var variance of the estimated quantile.

Examples

library("evd")
# Seasonal observations of 80 years at one station:
x1 <- rgev(80, 2, 1, 0.2) # observations from season 1
x2 <- rgev(80, 3, 1, 0.3) # observations from season 2
confInt_sTL(x1=x1, x2=x2, p=0.95, alpha=0.05)

# Seasonal observations of 100 years at 4 stations:
x1 <- matrix(rgev(400, 2, 1, 0.3), ncol=4)
x2 <- matrix(rgev(400, 4, 1, 0.2), ncol=4)
confInt_sTL(x1=x1, x2=x2, j=2, p=0.95, alpha=0.05)

Example output

$ci
[1]  7.103154 12.812078

$quant
[1] 9.957616

$var
[1] 169.6846

$ci
[1]  7.749313 10.121882

$quant
[1] 8.935598

$var
[1] 36.63377

flood documentation built on May 2, 2019, 4:04 p.m.