Shifted log-logistic distribution

Shifted log-logistic

Probability density function ${\displaystyle \mu =0,\sigma =1,}$ values of ${\displaystyle \xi }$ as shown in legend
Cumulative distribution function ${\displaystyle \mu =0,\sigma =1,}$ values of ${\displaystyle \xi }$ as shown in legend
Parameters	${\displaystyle \mu \in (-\infty ,+\infty )\,}$ location (real) ${\displaystyle \sigma \in (0,+\infty )\,}$ scale (real) ${\displaystyle \xi \in (-\infty ,+\infty )\,}$ shape (real)
Support	${\displaystyle x\geqslant \mu -\sigma /\xi \,\;(\xi >0)}$ ${\displaystyle x\leqslant \mu -\sigma /\xi \,\;(\xi <0)}$ ${\displaystyle x\in (-\infty ,+\infty )\,\;(\xi =0)}$
PDF	${\displaystyle {\frac {(1+\xi z)^{-(1/\xi +1)}}{\sigma \left(1+(1+\xi z)^{-1/\xi }\right)^{2}}}}$ where ${\displaystyle z=(x-\mu )/\sigma \,}$
CDF	${\displaystyle \left(1+(1+\xi z)^{-1/\xi }\right)^{-1}\,}$ where ${\displaystyle z=(x-\mu )/\sigma \,}$
Mean	${\displaystyle \mu +{\frac {\sigma }{\xi }}(\alpha \csc(\alpha )-1)}$ where ${\displaystyle \alpha =\pi \xi \,}$
Median	${\displaystyle \mu \,}$
Mode	${\displaystyle \mu +{\frac {\sigma }{\xi }}\left[\left({\frac {1-\xi }{1+\xi }}\right)^{\xi }-1\right]}$
Variance	${\displaystyle {\frac {\sigma ^{2}}{\xi ^{2}}}[2\alpha \csc(2\alpha )-(\alpha \csc(\alpha ))^{2}]}$ where ${\displaystyle \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 ${\displaystyle \delta }$. Thus if ${\displaystyle X}$ has a log-logistic distribution then ${\displaystyle X+\delta }$ has a shifted log-logistic distribution. So ${\displaystyle Y}$ has a shifted log-logistic distribution if ${\displaystyle \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 (CDF) of the shifted log-logistic distribution is

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

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

The probability density function (PDF) is

${\displaystyle 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 ${\displaystyle 1+\xi (x-\mu )/\sigma \geqslant 0.}$

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

Related distributions

• When ${\displaystyle \mu =\sigma /\xi ,}$ the shifted log-logistic reduces to the log-logistic distribution.
• When ${\displaystyle \xi }$ → 0, the shifted log-logistic reduces to the logistic distribution.
• The shifted log-logistic with shape parameter ${\displaystyle \xi =1}$ is the same as the generalized Pareto distribution with shape parameter ${\displaystyle \xi =1.}$

Applications

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

Alternate parameterization

An alternate parameterization with simpler expressions for the PDF and CDF is as follows. For the shape parameter ${\displaystyle \alpha }$, scale parameter ${\displaystyle \beta }$ and location parameter ${\displaystyle \gamma }$, the PDF is given by ^{[6] [7]}

${\displaystyle f(x)={\frac {\alpha }{\beta }}{\bigg (}{\frac {x-\gamma }{\beta }}{\bigg )}^{\alpha -1}{\bigg (}1+{\bigg (}{\frac {x-\gamma }{\beta }}{\bigg )}^{\alpha }{\bigg )}^{-2}}$

The CDF is given by

${\displaystyle F(x)={\bigg (}1+{\bigg (}{\frac {\beta }{x-\gamma }}{\bigg )}^{\alpha }{\bigg )}^{-1}}$

The mean is ${\displaystyle \beta \theta \csc(\theta )+\gamma }$ and the variance is ${\displaystyle \beta ^{2}\theta [2\csc(2\theta )-\theta \csc ^{2}(\theta )]}$, where ${\displaystyle \theta ={\frac {\pi }{\alpha }}}$. ^[7]

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. 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. 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  
6. "EasyFit - Log-Logistic Distribution" . Retrieved 1 October 2016. 
7. "Guide to Using - RISK7_EN.pdf" (PDF). Retrieved 1 October 2016.