Logistic distribution explained

In probability theory and statistics, the logistic distribution is a continuous probability distribution. Its cumulative distribution function is the logistic function, which appears in logistic regression and feedforward neural networks. It resembles the normal distribution in shape but has heavier tails (higher kurtosis). The logistic distribution is a special case of the Tukey lambda distribution.

Specification

Probability density function

When the location parameter  is 0 and the scale parameter  is 1, then the probability density function of the logistic distribution is given by

\begin{align} f(x;0,1)&=

e-x
(1+e-x)2

\\[4pt] &=

1
(ex/2+e-x/2)2

\\[5pt] &=

1
4

\operatorname{sech}2\left(

x
2

\right). \end{align}

Thus in general the density is:

\begin{align} f(x;\mu,s)&=

e-(x-\mu)/s
s\left(1+e-(x-\mu)/s\right)2

\\[4pt] &=

1
s\left(e(x-\mu)/(2s)+e-(x-\mu)/(2s)\right)2

\\[4pt] &=

1
4s
2\left(x-\mu
2s
\operatorname{sech}

\right). \end{align}

Because this function can be expressed in terms of the square of the hyperbolic secant function "sech", it is sometimes referred to as the sech-square(d) distribution.[1] (See also: hyperbolic secant distribution).

Cumulative distribution function

The logistic distribution receives its name from its cumulative distribution function, which is an instance of the family of logistic functions. The cumulative distribution function of the logistic distribution is also a scaled version of the hyperbolic tangent.

F(x;\mu,s)=

1
1+e-(x-\mu)/s

=

12
+
12\left(
\operatorname{tanh}
x-\mu
2s

\right).

In this equation is the mean, and is a scale parameter proportional to the standard deviation.

Quantile function

The inverse cumulative distribution function (quantile function) of the logistic distribution is a generalization of the logit function. Its derivative is called the quantile density function. They are defined as follows:

Q(p;\mu,s)=\mu+sln\left(

p
1-p

\right).

Q'(p;s)=

s
p(1-p)

.

Alternative parameterization

An alternative parameterization of the logistic distribution can be derived by expressing the scale parameter,

s

, in terms of the standard deviation,

\sigma

, using the substitution

s=q\sigma

, where

q=\sqrt{3}/{\pi}=0.551328895\ldots

. The alternative forms of the above functions are reasonably straightforward.

Applications

The logistic distribution—and the S-shaped pattern of its cumulative distribution function (the logistic function) and quantile function (the logit function)—have been extensively used in many different areas.

Logistic regression

One of the most common applications is in logistic regression, which is used for modeling categorical dependent variables (e.g., yes-no choices or a choice of 3 or 4 possibilities), much as standard linear regression is used for modeling continuous variables (e.g., income or population). Specifically, logistic regression models can be phrased as latent variable models with error variables following a logistic distribution. This phrasing is common in the theory of discrete choice models, where the logistic distribution plays the same role in logistic regression as the normal distribution does in probit regression. Indeed, the logistic and normal distributions have a quite similar shape. However, the logistic distribution has heavier tails, which often increases the robustness of analyses based on it compared with using the normal distribution.

Physics

The PDF of this distribution has the same functional form as the derivative of the Fermi function. In the theory of electron properties in semiconductors and metals, this derivative sets the relative weight of the various electron energies in their contributions to electron transport. Those energy levels whose energies are closest to the distribution's "mean" (Fermi level) dominate processes such as electronic conduction, with some smearing induced by temperature.[2] Note however that the pertinent probability distribution in Fermi–Dirac statistics is actually a simple Bernoulli distribution, with the probability factor given by the Fermi function.

The logistic distribution arises as limit distribution of a finite-velocity damped random motion described by a telegraph process in which the random times between consecutive velocity changes have independent exponential distributions with linearly increasing parameters.[3]

Hydrology

In hydrology the distribution of long duration river discharge and rainfall (e.g., monthly and yearly totals, consisting of the sum of 30 respectively 360 daily values) is often thought to be almost normal according to the central limit theorem.[4] The normal distribution, however, needs a numeric approximation. As the logistic distribution, which can be solved analytically, is similar to the normal distribution, it can be used instead. The blue picture illustrates an example of fitting the logistic distribution to ranked October rainfalls—that are almost normally distributed—and it shows the 90% confidence belt based on the binomial distribution. The rainfall data are represented by plotting positions as part of the cumulative frequency analysis.

Chess ratings

The United States Chess Federation and FIDE have switched its formula for calculating chess ratings from the normal distribution to the logistic distribution; see the article on Elo rating system (itself based on the normal distribution).

Related distributions

X\simLogistic(\mu,s)

then

kX+\ell\simLogistic(k\mu+\ell,|k|s)

.

X\sim

U(0, 1) then

\mu+s(logX-log(1-X))\simLogistic(\mu,s)

.

X\simGumbel(\muX,\beta)

and

Y\simGumbel(\muY,\beta)

independently then

X-Y\simLogistic(\muX-\muY,\beta)

.

X

and

Y\simGumbel(\mu,\beta)

then

X+Y\nsimLogistic(2\mu,\beta)

(The sum is not a logistic distribution). Note that

E(X+Y)=2\mu+2\beta\gamma2\mu=E\left(Logistic(2\mu,\beta)\right)

.

\left(\alpha=e\mu,\beta=

1
s

\right)

, and exp(X) + γ ~ shifted log-logistic

\left(\alpha=e\mu,\beta=

1
s,

\gamma\right)

.

\mu+slog(eX-1)\sim\operatorname{Logistic}(\mu,s).

\mu+slog\left(X
Y

\right)\sim\operatorname{Logistic}(\mu,s).

p

are substituted for logistic parameters

\mu

and

\sigma

. The resulting metalog quantile function is highly shape flexible, has a simple closed form, and can be fit to data with linear least squares.

Derivations

Higher-order moments

The nth-order central moment can be expressed in terms of the quantile function:

\begin{align} \operatorname{E}[(X-\mu)n]&=

infty
\int
-infty

(x-\mu)ndF(x)\\ &=

1(Q(p)-\mu)
\int
0

ndp=sn

1
\int\left[ln\left(
0
p
1-p

\right)\right]ndp. \end{align}

This integral is well-known and can be expressed in terms of Bernoulli numbers:

\operatorname{E}[(X-\mu)n]=sn\pin(2

n-2)|B
n|.

See also

References

Notes and References

  1. Johnson, Kotz & Balakrishnan (1995, p.116).
  2. Book: 9780521484916 . The Physics of Low-dimensional Semiconductors: An Introduction . Davies . John H. . 1998 . Cambridge University Press .
  3. A. Di Crescenzo, B. Martinucci (2010) "A damped telegraph random process with logistic stationary distribution", J. Appl. Prob., vol. 47, pp. 84–96.
  4. Book: Ritzema. H.P.. Frequency and Regression Analysis. 1994. Chapter 6 in: Drainage Principles and Applications, Publication 16, International Institute for Land Reclamation and Improvement (ILRI), Wageningen, The Netherlands. 175–224. 90-70754-33-9.