dgev.d: d-GEV probability density function

Description Usage Arguments Details Value See Also Examples

View source: R/d-gev.R

Description

Probability density function of duration-dependent GEV distribution

Usage

1
dgev.d(q, mut, sigma0, xi, theta, eta, d, eta2 = NULL, tau = 0, ...)

Arguments

q

vector of quantiles

mut, sigma0, xi

numeric value, giving modified location \tilde{μ}, scale offset σ_0 and shape parameter ξ.

theta

numeric value, giving duration offset θ (defining curvature of the IDF curve)

eta

numeric value, giving duration exponent η (defining slope of the IDF curve)

d

positive numeric value, giving duration

eta2

numeric value, giving a second duration exponent η_2 (needed for multiscaling). Default: NULL, treated as η_2=η.

tau

numeric value, giving intensity offset τ (defining flattening of the IDF curve). Default: τ=0.

...

additional parameters passed to dgev

Details

For details on the d-GEV and the parameter definitions, see IDF-package.

Value

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.

See Also

pgev.d, qgev.d, rgev.d

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
x <- seq(4,20,0.1)
# calculate probability density for one duration
dgev.d(q=x,mut=4,sigma0=2,xi=0,theta=0.1,eta=0.1,d=1)

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

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)
}
legend('topright',title = 'Duration',legend = 1:4,lty=1:4)

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