SkewLaplace  R Documentation 
Density function, distribution function, quantiles and random number generation for the skewLaplace distribution.
dskewlap(x, Theta, logPars = FALSE)
pskewlap(q, Theta)
qskewlap(p, Theta)
rskewlap(n, Theta)
x, q 
Vector of quantiles. 
p 
Vector of probabilities. 
n 
Number of observations to be generated. 
Theta 
Vector of parameters of the skewLaplace distribution:

.
logPars 
Logical. If 
The central skewLaplace has mode zero, and is a mixture of a (negative)
exponential distribution with mean \beta
, and the negative of an
exponential distribution with mean \alpha
. The weights of
the positive and negative components are proportional to their means.
The general skewLaplace distribution is a shifted central skewLaplace
distribution, where the mode is given by \mu
.
The density is given by:
f(x)=\frac{1}{\alpha+\beta} e^{(x  \mu)/\alpha}
for x\leq\mu
, and
f(x)=\frac{1}{\alpha+\beta} e^{(x  \mu)/\beta}
for x\geq\mu
dskewlap
gives the density, pskewlap
gives the distribution
function, qskewlap
gives the quantile function and rskewlap
generates random variates. The distribution function is obtained by
elementary integration of the density function. Random variates are
generated from exponential observations using the characterization of
the skewLaplace as a mixture of exponential observations.
David Scott d.scott@auckland.ac.nz, AiWei Lee, Richard Trendall
Fieller, N. J., Flenley, E. C. and Olbricht, W. (1992) Statistics of particle size data. Appl. Statist., 41, 127–146.
hyperbFitStart
Theta < c(1,2,1)
par(mfrow = c(1,2))
curve(dskewlap(x, Theta), from = 5, to = 8, n = 1000)
title("Density of the\n SkewLaplace Distribution")
curve(pskewlap(x, Theta), from = 5, to = 8, n = 1000)
title("Distribution Function of the\n SkewLaplace Distribution")
dataVector < rskewlap(500, Theta)
curve(dskewlap(x, Theta), range(dataVector)[1], range(dataVector)[2],
n = 500)
hist(dataVector, freq = FALSE, add =TRUE)
title("Density and Histogram\n of the SkewLaplace Distribution")
logHist(dataVector, main =
"LogDensity and LogHistogram\n of the SkewLaplace Distribution")
curve(log(dskewlap(x, Theta)), add = TRUE,
range(dataVector)[1], range(dataVector)[2], n = 500)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.