law0040.Log-Pareto-tail-normal: The Log-Pareto-tail-normal Distribution

In PoweR: Computation of Power and Level Tables for Hypothesis Tests

Description

Random generation for the Log-Pareto-tail-normal distribution with parameters alpha, mu and sigma.

This generator is called by function gensample to create random variables based on its parameters.

Details

If alpha, mu and sigma are not specified they assume the default values of 1.959964, 0.0 and 1.0 respectively.

The log-Pareto-tailed normal distribution has a symmetric and continuous density that belongs to the larger family of log-regularly varying distributions (see Desgagne, 2015). This is essentially a normal density with log-Pareto tails. Using this distribution instead of the usual normal ensures whole robustness to outliers in the estimation of location and scale parameters and in the estimation of parameters of a multiple linear regression.

The density of the log-Pareto-tailed normal distribution with parameters alpha, mu and sigma is given by

g(y|α,μ,σ) = (1/σ)φ((y-μ)(σ)) if μ - α σ <= y <= μ + α σ, and g(y|α,μ,σ) = φ(α) (α)/(|y-μ|) ((\log α)/(\log (|y-μ|/σ)))^β if |y-μ| >= α σ.

where β = 1+2\,φ(α)\,α\log(α)(1-q)^{-1} and q=Φ(α)-Φ(-α). The functions φ(α)=\frac{1}{√{2π}}\exp[-\frac{α^2}{2}] and Φ(α) are respectively the p.d.f. and the c.d.f. of the standard normal distribution. The domains of the variable and the parameters are -∞<y<∞, α>1, -∞<μ<∞ and σ>0.

Note that the normalizing constant K_{(α,β)} (see Desgagne, 2015, Definition 3) has been set to 1. The desirable consequence is that the core of the density, between μ-ασ and μ+ασ, becomes exactly the density of the N(μ,σ^2). This mass of the density corresponds to q. It follows that the parameter β is no longer free and its value depends on α as given above.

For example, if we set α=1.959964, we obtain β=4.083613 and q=0.95 of the mass is comprised between μ-ασ and μ+ασ. Note that if one is more comfortable in choosing the central mass $q$ instead of choosing directly the parameter α, then it suffices to use the equation α=Φ^{-1}((1+q)/2), with the contrainst q>0.6826895 <==> α>1.

The mean and variance of Log-Pareto-tail-normal are not defined.

Author(s)

P. Lafaye de Micheaux

References

Pierre Lafaye de Micheaux, Viet Anh Tran (2016). PoweR: A Reproducible Research Tool to Ease Monte Carlo Power Simulation Studies for Goodness-of-fit Tests in R. Journal of Statistical Software, 69(3), 1–42. doi:10.18637/jss.v069.i03

Desgagne, Alain. Robustness to outliers in location-scale parameter model using log-regularly varying distributions. Ann. Statist. 43 (2015), no. 4, 1568–1595. doi:10.1214/15-AOS1316. http://projecteuclid.org/euclid.aos/1434546215.

See Also

See Distributions for other standard distributions.

Examples

res <- gensample(40, 10000, law.pars = c(1.959964, 0.0, 1.0))
res$law
res$law.pars

PoweR documentation built on May 2, 2019, 2:09 p.m.