IDF.plot {IDF}

Plotting of IDF curves at a chosen station

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)