IDF.plot.Rd 1.49 KB
Newer Older
Christoph Ritschel's avatar
Christoph Ritschel committed
1
% Generated by roxygen2: do not edit by hand
Rust Henning's avatar
Rust Henning committed
2
3
4
% Please edit documentation in R/IDF.R
\name{IDF.plot}
\alias{IDF.plot}
5
\title{Plotting of IDF curves at a chosen station}
Rust Henning's avatar
Rust Henning committed
6
\usage{
7
8
IDF.plot(durations, fitparams, probs = c(0.5, 0.9, 0.99), cols = 4:2,
  add = FALSE, legend = TRUE, ...)
Rust Henning's avatar
Rust Henning committed
9
10
}
\arguments{
11
\item{durations}{vector of durations for which to calculate the quantiles.}
Rust Henning's avatar
Rust Henning committed
12

13
14
\item{fitparams}{vector containing parameters mut, sigma0, xi, theta, eta
(modified location, scale, shape, duration offset, duration exponent) for chosen station
15
16
as obtained from \code{\link{gev.d.fit}}
(or \code{\link{gev.d.params}} for model with covariates).}
Rust Henning's avatar
Rust Henning committed
17

18
\item{probs}{vector of exeedance probabilities for which to plot IDF curves (p = 1-1/ReturnPeriod)}
Rust Henning's avatar
Rust Henning committed
19

20
\item{cols}{vector of colors for IDF curves. Should have same length as \code{probs}}
Christoph Ritschel's avatar
Christoph Ritschel committed
21

22
\item{add}{logical indicating if plot should be added to existing plot}
Rust Henning's avatar
Rust Henning committed
23

24
25
26
\item{legend}{logical indicating if legend should be plotted}

\item{...}{additional parameters passed on to the \code{plot} function}
Rust Henning's avatar
Rust Henning committed
27
28
}
\description{
29
Plotting of IDF curves at a chosen station
Rust Henning's avatar
Rust Henning committed
30
}
Christoph Ritschel's avatar
Christoph Ritschel committed
31
\examples{
32
data('example',package = 'IDF')
33
# fit d-gev
34
35
fit <- gev.d.fit(example$dat,example$d,ydat = as.matrix(example[,c("cov1","cov2")])
                 ,mul = c(1,2),sigl = 1)
36
# get parameters for cov1 = 1, cov2 = 1
37
par <- gev.d.params(fit = fit, ydat = matrix(1,1,2))
38
39
40
41
# 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)
Christoph Ritschel's avatar
Christoph Ritschel committed
42
}