bootstrapCI: Determine confidence intervals to rainfall intensities...

Description Usage Arguments Value Author(s) Examples

View source: R/bootstrapCI.R

Description

Conducts bootstrap to randomly sample of intensity values 'n' times for a specified distribution to estimate the confidence interval for each given non-exceedance probability. Based in code posted by Dave Hutchinson http://www.headwateranalytics.com.

Usage

1
2
bootstrapCI(Intensity = .., Parameters = .., Type = "Gumbel",
  Rsample = 1000, Return.P = ..., Conf.Inter = 0.95)

Arguments

Intensity:

a numeric vector with intensity [mm/h] values of different years for a specific time duration (e.g. 5, 15, 120 minutes, etc.)

Parameters:

list with three elements: i) type of distribution function ii) fitted parameters, and iii) source to call specfic function in the lmomco package.

Type:

a character specifying a name of the probability distribution function fitted (see selecDIST) by fitDISTRI function.

Rsample:

An integer representing number of resamples to conduct when confidence interval will be computed.

Return.P:

a numeric vector with return periods like non-exceedance probabilities.

Conf.Inter:

level of the confidence interval.

Value

A list of:

Author(s)

David Zamora <[email protected]> Water Resources Engineering Research Group - GIREH

Examples

1
2
3
4
5
6
7
8
9
data(inten)
data(Pargumbel)
Tp <- c(2, 3, 5, 10, 25, 50, 100)
FR <- 1 - 1/Tp
CI.test <- bootstrapCI(Intensity = inten, Parameters = Pargumbel, Type = "gumbel",
Rsample = 1000, Return.P =FR, Conf.Inter = 0.90)
(CI.test$Conf.Inter[,2])
# result = 125.9501, 139.6132, 154.6232, 173.6736, 197.5024, 215.0634, 232.4939
 

dazamora/IDFtool documentation built on March 20, 2018, 8:56 p.m.