dgev.d.Rd 1.68 KB
Newer Older
Christoph Ritschel's avatar
Christoph Ritschel committed
1
% Generated by roxygen2: do not edit by hand
2
% Please edit documentation in R/d-gev.R
Rust Henning's avatar
Rust Henning committed
3
4
\name{dgev.d}
\alias{dgev.d}
5
\title{d-GEV probability density function}
Rust Henning's avatar
Rust Henning committed
6
\usage{
7
dgev.d(q, mut, sigma0, xi, theta, eta, tau, d, ...)
Rust Henning's avatar
Rust Henning committed
8
9
}
\arguments{
10
\item{q}{vector of quantiles}
Rust Henning's avatar
Rust Henning committed
11

12
13
\item{mut, sigma0, xi}{numeric value, giving modified location \eqn{\tilde{\mu}}, scale offset \eqn{\tilde{\sigma_0}} and 
shape parameter \eqn{\xi}.}
Rust Henning's avatar
Rust Henning committed
14

15
\item{theta}{numeric value, giving duration offset \eqn{\theta} (defining curvature of the IDF curve)}
Rust Henning's avatar
Rust Henning committed
16

17
\item{eta}{numeric value, giving duration exponent \eqn{\eta} (defining slope of the IDF curve)}
Rust Henning's avatar
Rust Henning committed
18

19
20
\item{tau}{numeric value, giving intensity offset \eqn{\tau} (defining flattening of the IDF curve)}

Laura Mack's avatar
Laura Mack committed
21
\item{d}{positive numeric value, giving duration}
Rust Henning's avatar
Rust Henning committed
22

23
\item{...}{additional parameters passed to \code{\link[evd]{dgev}}}
Rust Henning's avatar
Rust Henning committed
24
25
}
\value{
26
27
list containing vectors of density values for given quantiles.
The first element of the list are the density values for the first given duration etc.
Rust Henning's avatar
Rust Henning committed
28
29
}
\description{
30
Probability density function of duration-dependent GEV distribution
Rust Henning's avatar
Rust Henning committed
31
}
32
\details{
33
For details on the d-GEV and the parameter definitions, see \link{IDF-package}.
34
35
36
}
\examples{
x <- seq(4,20,0.1)
37
# calculate probability density for one duration
38
dgev.d(q=x,mut=4,sigma0=2,xi=0,theta=0.1,eta=0.1,tau=0.1,d=1)
39
40
41

# calculate probability density for different durations
ds <- 1:4
42
dens <- lapply(ds,dgev.d,q=x,mut=4,sigma0=2,xi=0,theta=0.1,eta=0.1, tau=0.1)
43
44
45
46
47

plot(x,dens[[1]],type='l',ylim = c(0,0.21),ylab = 'Probability Density')
for(i in 2:4){
  lines(x,dens[[i]],lty=i)
}
48
legend('topright',title = 'Duration',legend = 1:4,lty=1:4)
49
}
Rust Henning's avatar
Rust Henning committed
50
\seealso{
51
\code{\link{pgev.d}}, \code{\link{qgev.d}}, \code{\link{rgev.d}}
Rust Henning's avatar
Rust Henning committed
52
}