daep {AEP}R Documentation

Computing the probability density function (pdf) of asymmetric exponential power (AEP) distribution.

Description

The pdf of AEP distribution given by

fX(xΘ)=12σΓ(1+1α)exp{μxσ(1ϵ)α},   x<μ, f_{X}(x|\Theta)= \frac{1}{2\sigma \Gamma\bigl(1+\frac{1}{\alpha}\bigr)}\exp\biggl\{-\bigg|\frac{\mu-x}{\sigma(1-\epsilon)}\bigg|^{\alpha}\biggr\},~~{}~x < \mu,

fX(xΘ)=12σΓ(1+1α)exp{xμσ(1+ϵ)α},   xμ, f_{X}(x|\Theta)= \frac{1}{2\sigma \Gamma\bigl(1+\frac{1}{\alpha}\bigr)}\exp\biggl\{-\bigg|\frac{x-\mu}{\sigma(1+\epsilon)}\bigg|^{\alpha}\biggr\},~~{}~x \geq\mu,

where <x<+-\infty<x<+\infty, Θ=(α,σ,μ,ϵ)T\Theta=(\alpha,\sigma,\mu,\epsilon)^T with 0<α20<\alpha \leq 2, σ>0\sigma> 0, <μ<-\infty<\mu<\infty, 1<ϵ<1-1<\epsilon<1, and

Γ(u)=0+xu1exp{x}dx, u>0.\Gamma(u)=\int_{0}^{+\infty} x^{u-1}\exp\bigl\{-x\bigr\}dx,~u>0.

Usage

daep(x, alpha, sigma, mu, epsilon, log = FALSE)

Arguments

x

Vector of observation of requested random realizations.

alpha

Tail thickness parameter.

sigma

Scale parameter.

mu

Location parameter.

epsilon

Skewness parameter.

log

If TRUE, then log(fX(xΘ))\bigl(f_{X}(x|\Theta)\bigr) is returned.

Details

The AEP distribution is a special case of asymmetric exponential power distribution studied by Dongming and Zinde-Walsh (2009) when p1=p2=αp_1=p_2=\alpha. Also, note that if ϵ=0\epsilon=0, then the AEP distribution turns into a normal distribution with mean μ\mu and standard deviation 2σ\sqrt{2}\sigma. When α=2\alpha=2, the AEP distribution is a slight variant of epsilon-skew-normal distribution introduced by Mudholkar and Huston (2001).

Value

Computed pdf of AEP distribution at points of vector xx.

Author(s)

Mahdi Teimouri

References

Z. Dongming and V. Zinde-Walsh, 2009. Properties and estimation of asymmetric exponential power distribution, Journal of Econometrics, 148(1), 86-99.

G. S. Mudholkar and A. D. Huston, 2001. The epsilon-skew-normal distribution for analyzing near-normal data, Journal of Statistical Planning and Inference, 83, 291-309.

Examples

daep(x = 2, alpha = 1.5, sigma = 1, mu = 0, epsilon = 0.5, log = FALSE)

[Package AEP version 0.1.4 Index]