IDFCurve: Compute Intensity-Duration-Frequency curve

In dazamora/IDFtool: Estimation Intensity-Duration-Frequency Curves for Rainfall Series

Compute Intensity-Duration-Frequency curve

Description

An Intensity-Duration-Frequency curve (IDF Curve) is a graphical representation of the probability that a given average rainfall intensity will occur. This function allows to fit different probability distribution functions (see selecDIST) by means of four fit methods (see fitDISTRI) to determine intensity [mm/h] for different return periods and per specific time durations. Finally, it computes equations parameters of the IDF curves (see regIDF)

Usage

IDFCurve(Data, Station = "2610516", Duration = FALSE,
  Periods = FALSE, Type = "gumbel", M.fit = "lmoments",
  Plot = 1234, Strategy = 1:3, logaxe = "", CI = FALSE,
  CIpdf = TRUE, iter = 500, goodtest = FALSE, Resolution = 300,
  SAVE = FALSE, name = TRUE)

Arguments

Data	a numeric matrix with years in the first column and intensity values for each Duration time in the other columns.
Station	a string with a name to identify source of Data.
Duration	a logical value or numeric vector. If it is TRUE the durations (in minutes) by default will be used: 5, 10 ,15, 20, 30, 60, 120 y 360. In case of a numeric vector the durations must be in minutes.
Periods	a logical value or numeric vector.
Type	a character specifying the name of the distribution function that will be employed: exponencial, gamma, gev, gumbel, log.normal3, normal, pearson, log.pearson3 and wakeby (see selecDIST).
M.fit	a character specifying a name of fit method employed on pdf, just three options are available: L-moments (Lmoments), Probability-Weighted Moments (PWD), Maximum Likelihood (MLEZ) and Moments (MME) (see MME_DIST).
Plot	it is a number of one to four digits. a number (1) to plot density curves both empirical as modeled (pdf).(2) to plot if it will be plotted curves between return Periods and intensity computed by pdf fitted. Or use (12) to get both graphs. (3) to plot IDF curves for all return periods: (Durations versus Intensity) for all return Periods. (4) to plot IDF curve each for return Periods with its confidence and prediction intervals. Or use (34) to get both graphs. If other is used the graphs will not appear.
Strategy	a numeric vector used to identify Strategies to compute IDF curves with different data sets: 1 just data from HIDFUN tool, 2 just data from Ideam and 3 used this data sets. These strategies were created just for Colombian pluviograph data. For your data use strategy number 1.
logaxe	a character to plot axis in log scale: x, y or both (xy). In other case used "".
CI	a logical value specifying whether confidence and prediction intervals will be computed in IDF curves.
CIpdf	a logical value specifying whether confidence of pdf will be computed.
iter	an integer representing number of resamples to conduct when confidence interval is computed (see bootstrapCI). Use it only if CI is equal to TRUE.
goodtest	a logical value specifying whether goodness-fit tests should be cumputed to pdf fitted by means of goodFIT function.
Resolution	a number to determine the resolution that the plot function will used to save graphs. It has two options: 300 and 600 ppi. See resoPLOT.
SAVE	a logical value. TRUE will save Plot, FALSE will just show Plot.
name	a logical value. TRUE will use a default names to identify strategies: (1) "HIDFUN", (2) "IDEAM", (3) "AMBOS". In other case FALSE allows: i) selected years of data sets, and ii) to insert durations to do IDF curves (in minutes)

Value

A list of:

• Intesities a numeric matrix of intensities values per each return Periods computed by pdf fitted.

• Models a list with results of the function regIDF.

• Test.fit a list with results of the function goodFIT.

• Distribution a list with results of the function fitDISTRI.

Author(s)

David Zamora <[e-mail]> Albeiro Figueroa <[e-mail]> Water Resources Engineering Research Group - GIREH

Examples

# Meteorology station in the Farfan Airport in Tulua, Colombia.
data(inten)
Test.idftool <- IDFCurve(Data = inten, Station='2610516', Duration = FALSE,
Periods = FALSE, Type = "gumbel", M.fit = "lmoments",
Plot = 1234, Strategy = 1, logaxe = "", CI = FALSE, CIpdf = TRUE, iter = 50,
goodtest = FALSE, Resolution = 300, SAVE = FALSE, name = TRUE)

dazamora/IDFtool documentation built on Jan. 1, 2023, 3:29 p.m.