Shifted log-logistic distribution

Shifted log-logistic
Probability density function
\mu=0, \sigma=1, values of \xi as shown in legend
Cumulative distribution function
\mu=0, \sigma=1, values of \xi as shown in legend
Parameters

\mu \in (-\infty,+\infty) \, location (real)
\sigma \in (0,+\infty)    \, scale (real)

\xi\in (-\infty,+\infty)  \, shape (real)
Support

x \geqslant \mu -\sigma/\xi\,\;(\xi > 0)
x \leqslant \mu -\sigma/\xi\,\;(\xi < 0)

x \in (-\infty, +\infty)   \,\;(\xi = 0)
PDF

\frac{(1+\xi z)^{-(1/\xi +1)}}{\sigma\left(1 + (1+\xi z)^{-1/\xi}\right)^2}

where z=(x-\mu)/\sigma\,
CDF

\left(1+(1 + \xi z)^{-1/\xi}\right)^{-1} \,

where z=(x-\mu)/\sigma\,
Mean

\mu + \frac{\sigma}{\xi}(\alpha \csc(\alpha)-1)

where \alpha= \pi \xi\,
Median \mu \,
Mode \mu + \frac{\sigma}{\xi}\left[\left(\frac{1-\xi}{1+\xi}\right)^\xi - 1 \right]
Variance

 \frac{\sigma^2}{\xi^2}[2\alpha \csc(2 \alpha) - (\alpha \csc(\alpha))^2]

where \alpha= \pi \xi\,

The shifted log-logistic distribution is a probability distribution also known as the generalized log-logistic or the three-parameter log-logistic distribution.[1][2] It has also been called the generalized logistic distribution,[3] but this conflicts with other uses of the term: see generalized logistic distribution.

Definition

The shifted log-logistic distribution can be obtained from the log-logistic distribution by addition of a shift parameter \delta. Thus if X has a log-logistic distribution then X+\delta has a shifted log-logistic distribution. So Y has a shifted log-logistic distribution if \log(Y-\delta) has a logistic distribution. The shift parameter adds a location parameter to the scale and shape parameters of the (unshifted) log-logistic.

The properties of this distribution are straightforward to derive from those of the log-logistic distribution. However, an alternative parameterisation, similar to that used for the generalized Pareto distribution and the generalized extreme value distribution, gives more interpretable parameters and also aids their estimation.

In this parameterisation, the cumulative distribution function of the shifted log-logistic distribution is

F(x; \mu,\sigma,\xi) = \frac{1}{ 1 + \left(1+ \frac{\xi(x-\mu)}{\sigma}\right)^{-1/\xi}}

for  1 + \xi(x-\mu)/\sigma \geqslant 0, where \mu\in\mathbb R is the location parameter, \sigma>0\, the scale parameter and \xi\in\mathbb R the shape parameter. Note that some references use  \kappa =  - \xi\,\! to parameterise the shape.[3][4]

The probability density function is

 f(x; \mu,\sigma,\xi) = \frac{\left(1+\frac{\xi(x-\mu)}{\sigma}\right)^{-(1/\xi +1)}}
{\sigma\left[1 + \left(1+\frac{\xi(x-\mu)}{\sigma}\right)^{-1/\xi}\right]^2},

again, for  1 + \xi(x-\mu)/\sigma \geqslant 0.

The shape parameter \xi is often restricted to lie in [-1,1], when the probability density function is bounded. When |\xi|>1, it has an asymptote at x = \mu - \sigma/\xi. Reversing the sign of \xi reflects the pdf and the cdf about x=\mu..

Related distributions

Applications

The three-parameter log-logistic distribution is used in hydrology for modelling flood frequency.[3][4][5]

References

  1. Venter, Gary G. (Spring 1994), "Introduction to selected papers from the variability in reserves prize program" (PDF), Casualty Actuarial Society Forum 1: 91–101
  2. Geskus, Ronald B. (2001), "Methods for estimating the AIDS incubation time distribution when date of seroconversion is censored", Statistics in Medicine 20 (5): 795–812, doi:10.1002/sim.700, PMID 11241577
  3. 1 2 3 Hosking, Jonathan R. M.; Wallis, James R (1997), Regional Frequency Analysis: An Approach Based on L-Moments, Cambridge University Press, ISBN 0-521-43045-3
  4. 1 2 Robson, A.; Reed, D. (1999), Flood Estimation Handbook, 3: "Statistical Procedures for Flood Frequency Estimation", Wallingford, UK: Institute of Hydrology, ISBN 0-948540-89-3
  5. Ahmad, M. I.; Sinclair, C. D.; Werritty, A. (1988), "Log-logistic flood frequency analysis", Journal of Hydrology 98 (3–4): 205–224, doi:10.1016/0022-1694(88)90015-7
This article is issued from Wikipedia - version of the Tuesday, December 10, 2013. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.