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

Description Usage Arguments Examples

View source: R/IDF.R

Description

Plotting of IDF curves at a chosen station

Usage

1
2
3
4
5
6
7
8
9
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
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.