IDF.plot: Plotting of IDF curves at a chosen station

In IDF: Estimation and Plotting of IDF Curves

Description

Plotting of IDF curves at a chosen station

Usage

IDF.plot(
  durations,
  fitparams,
  probs = c(0.5, 0.9, 0.99),
  cols = 4:2,
  add = FALSE,
  legend = TRUE,
  ...
)

Arguments

durations	vector of durations for which to calculate the quantiles.
fitparams	vector containing parameters mut, sigma0, xi, theta, eta (modified location, scale offset, shape, duration offset, duration exponent) for chosen station as obtained from gev.d.fit (or gev.d.params for model with covariates).
probs	vector of non-exceedance probabilities for which to plot IDF curves (p = 1-1/(Return Period))
cols	vector of colors for IDF curves. Should have same length as probs
add	logical indicating if plot should be added to existing plot, default is FALSE
legend	logical indicating if legend should be plotted (TRUE, the default)
...	additional parameters passed on to the plot function

Examples

data('example',package = 'IDF')
# fit d-gev
fit <- gev.d.fit(example$dat,example$d,ydat = as.matrix(example[,c("cov1","cov2")])
                 ,mutl = c(1,2),sigma0l = 1)
# get parameters for cov1 = 1, cov2 = 1
par <- gev.d.params(fit = fit, ydat = matrix(1,1,2))
# plot quantiles
IDF.plot(durations = seq(0.5,35,0.2),fitparams = par)
# add data points
points(example[example$cov1==1,]$d,example[example$cov1==1,]$dat)

IDF documentation built on April 9, 2021, 5:07 p.m.