The Logistic Distribution

The logistic distribution has been used for growth models, and is used in a certain type of regression known as the logistic regression. It has also applications in modeling life data. The shape of the logistic distribution and the normal distribution are very similar, as discussed in Meeker and Escobar [27]. There are some who argue that the logistic distribution is inappropriate for modeling lifetime data because the left-hand limit of the distribution extends to negative infinity. This could conceivably result in modeling negative times-to-failure. However, provided that the distribution in question has a relatively high mean and a relatively small location parameter, the issue of negative failure times should not present itself as a problem.

Logistic Probability Density Function
The logistic pdf is given by:


 * $$\begin{matrix}

f(t)=\tfrac{\sigma {{(1+{{e}^{z}})}^{2}}} \\ \\

z=\tfrac{t-\mu }{\sigma } \\ \\

-\infty 0 \\ \\ \end{matrix}\,\!$$

where:


 * $$\begin{align}

\mu = & \text{location parameter (also denoted as }\overline{T}) \\ \sigma = & \text{scale parameter} \end{align}\,\!$$

The Logistic Mean, Median and Mode
The logistic mean or MTTF is actually one of the parameters of the distribution, usually denoted as $$\mu \,\!$$. Since the logistic distribution is symmetrical, the median and the mode are always equal to the mean, $$\mu =\tilde{T}=\breve{T}.\,\!$$



The Logistic Standard Deviation
The standard deviation of the logistic distribution, $${{\sigma }_{T}}\,\!$$, is given by:


 * $${{\sigma }_{T}}=\sigma \pi \frac{\sqrt{3}}{3}\,\!$$

The Logistic Reliability Function
The reliability for a mission of time $$t\,\!$$, starting at age 0, for the logistic distribution is determined by:


 * $$R(t)=\int_{t}^{\infty }f(t)dt\,\!$$

or:


 * $$R(t)=\frac{1}{1+{{e}^{z}}}\,\!$$

The unreliability function is:


 * $$F=\frac{1+{{e}^{z}}}\,\!$$

where:


 * $$z=\frac{t-\mu }{\sigma }\,\!$$

The Logistic Conditional Reliability Function
The logistic conditional reliability function is given by:


 * $$R(t|T)=\frac{R(T+t)}{R(T)}=\frac{1+{{e}^{\tfrac{T-\mu }{\sigma }}}}{1+{{e}^{\tfrac{t+T-\mu }{\sigma }}}}\,\!$$

The Logistic Reliable Life
The logistic reliable life is given by:


 * $$\begin{align}

{{T}_{R}}=\mu +\sigma[\ln (1-R)-\ln (R)] \end{align}\,\!$$

The Logistic Failure Rate Function
The logistic failure rate function is given by:


 * $$\lambda (t)=\frac{\sigma (1+{{e}^{z}})}\,\!$$

Characteristics of the Logistic Distribution

 * The logistic distribution has no shape parameter. This means that the logistic pdf has only one shape, the bell shape, and this shape does not change. The shape of the logistic distribution is very similar to that of the normal distribution.
 * The mean, $$\mu \,\!$$, or the mean life or the $$MTTF\,\!$$, is also the location parameter of the logistic pdf, as it locates the pdf along the abscissa. It can assume values of $$-\infty <\bar{T}<\infty \,\!$$.
 * As $$\mu \,\!$$ decreases, the pdf is shifted to the left.
 * As $$\mu \,\!$$ increases, the pdf is shifted to the right.




 * As $$\sigma \,\!$$ decreases, the pdf gets pushed toward the mean, or it becomes narrower and taller.
 * As $$\sigma \,\!$$ increases, the pdf spreads out away from the mean, or it becomes broader and shallower.
 * The scale parameter can assume values of $$0<\sigma <\infty \,\!$$.


 * The logistic pdf starts at $$t=-\infty \,\!$$ with an $$f(t)=0\,\!$$. As $$t\,\!$$ increases, $$f(t)\,\!$$ also increases, goes through its point of inflection and reaches its maximum value at $$t=\bar{T}\,\!$$. Thereafter, $$f(t)\,\!$$ decreases, goes through its point of inflection and assumes a value of $$f(t)=0\,\!$$ at $$t=+\infty \,\!$$.
 * For $$t=\pm \infty ,\,\!$$ the pdf equals 0. The maximum value of the pdf occurs at $$t = \mu \,\!$$ and equals $$\tfrac{1}{4\sigma }.\,\!$$
 * The point of inflection of the pdf plot is the point where the second derivative of the pdf equals zero. The inflection point occurs at $$t=\mu +\sigma \ln (2\pm \sqrt{3})\,\!$$ or $$t\approx \mu \pm \sigma 1.31696\,\!$$.
 * If the location parameter $$\mu \,\!$$ decreases, the reliability plot is shifted to the left. If $$\mu \,\!$$ increases, the reliability plot is shifted to the right.
 * If $$t=\mu \,\!$$ then $$R=0.5\,\!$$.    is the inflection point. If $$t<\mu \,\!$$  then $$R(t)\,\!$$ is concave (concave down); if $$t>\mu \,\!$$ then $$R(t)\,\!$$ is convex (concave up). For $$t<\mu ,\,\!$$ $$\lambda (t)\,\!$$ is convex (concave up), for $$t>\mu ;\,\!$$ $$\lambda (t)\,\!$$ is concave (concave down).
 * The main difference between the normal distribution and logistic distribution lies in the tails and in the behavior of the failure rate function. The logistic distribution has slightly longer tails compared to the normal distribution. Also, in the upper tail of the logistic distribution, the failure rate function levels out for large $$t\,\!$$ approaching 1/ $$\delta .\,\!$$
 * If location parameter $$\mu \,\!$$ decreases, the failure rate plot is shifted to the left. Vice versa if $$\mu \,\!$$ increases, the failure rate plot is shifted to the right.
 * $$\lambda \,\!$$ always increases. For $$t\to -\infty \,\!$$    for $$t\to \infty \,\!$$     It is always $$0\le \lambda (t)\le \tfrac{1}{\sigma }.\,\!$$
 * If $$\sigma \,\!$$ increases, then $$\lambda (t)\,\!$$ increases more slowly and smoothly. The segment of time where $$0<\lambda (t)<\tfrac{1}{\sigma }\,\!$$ increases, too, whereas the region where $$\lambda (t)\,\!$$ is close to 0 or $$\tfrac{1}{\sigma }\,\!$$ gets narrower. Conversely, if $$\sigma \,\!$$ decreases, then $$\lambda (t)\,\!$$ increases more quickly and sharply. The segment of time where $$0 < \lambda (t)<\tfrac{1}{\sigma }\,\!$$ decreases, too, whereas the region where $$\lambda (t)\,\!$$ is close to 0 or $$\tfrac{1}{\sigma }\,\!$$ gets broader.

Weibull++ Notes on Negative Time Values

One of the disadvantages of using the logistic distribution for reliability calculations is the fact that the logistic distribution starts at negative infinity. This can result in negative values for some of the results. Negative values for time are not accepted in most of the components of Weibull++, nor are they implemented. Certain components of the application reserve negative values for suspensions, or will not return negative results. For example, the Quick Calculation Pad will return a null value (zero) if the result is negative. Only the Free-Form (Probit) data sheet can accept negative values for the random variable (x-axis values).

Logistic Distribution Probability Paper
The form of the Logistic probability paper is based on linearizing the cdf. From unreliability equation, $$z\,\!$$ can be calculated as a function of the cdf $$F\,\!$$ as follows:


 * $$\begin{align}

z=\ln (F)-\ln (1-F) \end{align}\,\!$$

or using the equation for $$z\,\!$$


 * $$\frac{t-\mu }{\sigma }=\ln (F)-\ln (1-F)\,\!$$

Then:


 * $$\ln (F)-\ln (1-F)=-\frac{\mu }{\sigma }+\frac{1}{\sigma }t\,\!$$

Now let:


 * $$\begin{align}

y=\ln (F)-\ln (1-F) \end{align}\,\!$$


 * $$\begin{align}

x=t \end{align}\,\!$$

and:


 * $$a=-\frac{\mu }{\sigma }\,\!$$


 * $$b=\frac{1}{\sigma }\,\!$$

which results in the following linear equation:


 * $$\begin{align}

y=a+bx \end{align}\,\!$$

The logistic probability paper resulting from this linearized cdf function is shown next.



Since the logistic distribution is symmetrical, the area under the pdf curve from $$-\infty \,\!$$ to $$\mu \,\!$$ is 0.5, as is the area from $$\mu \,\!$$ to $$+\infty \,\!$$. Consequently, the value of $$\mu \,\!$$ is said to be the point where $$R(t)=Q(t)=50%\,\!$$. This means that the estimate of $$\mu \,\!$$ can be read from the point where the plotted line crosses the 50% unreliability line.

For $$z=1\,\!$$, $$\sigma =t-\mu \,\!$$ and $$R(t)=\tfrac{1}{1+\exp (1)}\approx 0.2689.\,\!$$ Therefore, $$\sigma \,\!$$ can be found by subtracting $$\mu \,\!$$ from the time value where the plotted probability line crosses the 73.10% unreliability (26.89% reliability) horizontal line.

Confidence Bounds
In this section, we present the methods used in the application to estimate the different types of confidence bounds for logistically distributed data. The complete derivations were presented in detail (for a general function) in Confidence Bounds.

Bounds on the Parameters
The lower and upper bounds on the location parameter $$\widehat{\mu }\,\!$$ are estimated from


 * $${{\mu }_{U}}=\widehat{\mu }+{{K}_{\alpha }}\sqrt{Var(\widehat{\mu })\text{ }}\text{ (upper bound)}\,\!$$


 * $${{\mu }_{L}}=\widehat{\mu }-{{K}_{\alpha }}\sqrt{Var(\widehat{\mu })\text{ }}\text{ (lower bound)}\,\!$$

The lower and upper bounds on the scale parameter $$\widehat{\sigma }\,\!$$ are estimated from:


 * $${{\sigma }_{U}}=\widehat{\sigma }{{e}^{\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{\sigma })\text{ }}}{\widehat{\sigma }}}}(\text{upper bound})\,\!$$


 * $${{\sigma }_{L}}=\widehat{\sigma }{{e}^{\tfrac{-{{K}_{\alpha }}\sqrt{Var(\widehat{\sigma })\text{ }}}{\widehat{\sigma }}}}\text{ (lower bound)}\,\!$$

where $${{K}_{\alpha }}\,\!$$ is defined by:


 * $$\alpha =\frac{1}{\sqrt{2\pi }}\int_^{\infty }{{e}^{-\tfrac{2}}}dt=1-\Phi ({{K}_{\alpha }})\,\!$$

If $$\delta \,\!$$ is the confidence level, then $$\alpha =\tfrac{1-\delta }{2}\,\!$$ for the two-sided bounds, and $$\alpha =1-\delta \,\!$$ for the one-sided bounds.

The variances and covariances of $$\widehat{\mu }\,\!$$ and $$\widehat{\sigma }\,\!$$ are estimated from the Fisher matrix, as follows:


 * $$\left( \begin{matrix}

\widehat{Var}\left( \widehat{\mu } \right) & \widehat{Cov}\left( \widehat{\mu },\widehat{\sigma } \right) \\ \widehat{Cov}\left( \widehat{\mu },\widehat{\sigma } \right) & \widehat{Var}\left( \widehat{\sigma } \right) \\ \end{matrix} \right)=\left( \begin{matrix} -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{\mu }^{2}}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial \mu \partial \sigma } \\ {} & {} \\   -\tfrac{{{\partial }^{2}}\Lambda }{\partial \mu \partial \sigma } & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{\sigma }^{2}}}  \\ \end{matrix} \right)_{\mu =\widehat{\mu },\sigma =\widehat{\sigma }}^{-1}\,\!$$

$$\Lambda \,\!$$ is the log-likelihood function of the normal distribution, described in Parameter Estimation and Appendix D.

Bounds on Reliability
The reliability of the logistic distribution is:


 * $$\widehat{R}=\frac{1}{1+{{e}^{\widehat{z}}}}\,\!$$

where:


 * $$\widehat{z}=\frac{t-\widehat{\mu }}{\widehat{\sigma }}\,\!$$

Here $$-\infty <t<\infty \,\!$$, $$-\infty <\mu <\infty \,\!$$, $$0<\sigma <\infty \,\!$$. Therefore, $$z\,\!$$ also is changing from $$-\infty \,\!$$ to $$+\infty \,\!$$. Then the bounds on $$z\,\!$$ are estimated from:


 * $${{z}_{U}}=\widehat{z}+{{K}_{\alpha }}\sqrt{Var(\widehat{z})\text{ }}\,\!$$


 * $${{z}_{L}}=\widehat{z}-{{K}_{\alpha }}\sqrt{Var(\widehat{z})\text{ }}\text{ }\,\!$$

where:


 * $$Var(\widehat{z})={{(\frac{\partial z}{\partial \mu })}^{2}}Var(\widehat{\mu })+2(\frac{\partial z}{\partial \mu })(\frac{\partial z}{\partial \sigma })Cov(\widehat{\mu },\widehat{\sigma })+{{(\frac{\partial z}{\partial \sigma })}^{2}}Var(\widehat{\sigma })\,\!$$

or:


 * $$Var(\widehat{z})=\frac{1}(Var(\widehat{\mu })+2\widehat{z}Cov(\widehat{\mu },\widehat{\sigma })+{{\widehat{z}}^{2}}Var(\widehat{\sigma }))\,\!$$

The upper and lower bounds on reliability are:


 * $${{R}_{U}}=\frac{1}{1+{{e}^}}\text{(upper bound)}\,\!$$


 * $${{R}_{L}}=\frac{1}{1+{{e}^}}\text{(lower bound)}\,\!$$

Bounds on Time
The bounds around time for a given logistic percentile (unreliability) are estimated by first solving the reliability equation with respect to time as follows:


 * $$\widehat{T}(\widehat{\mu },\widehat{\sigma })=\widehat{\mu }+\widehat{\sigma }z\,\!$$

where:


 * $$\begin{align}

z=\ln (1-R)-\ln (R) \end{align}\,\!$$


 * $$Var(\widehat{T})={{(\frac{\partial T}{\partial \mu })}^{2}}Var(\widehat{\mu })+2(\frac{\partial T}{\partial \mu })(\frac{\partial T}{\partial \sigma })Cov(\widehat{\mu },\widehat{\sigma })+{{(\frac{\partial T}{\partial \sigma })}^{2}}Var(\widehat{\sigma })\,\!$$

or:


 * $$Var(\widehat{T})=Var(\widehat{\mu })+2\widehat{z}Cov(\widehat{\mu },\widehat{\sigma })+{{\widehat{z}}^{2}}Var(\widehat{\sigma })\,\!$$

The upper and lower bounds are then found by:


 * $${{T}_{U}}=\widehat{T}+{{K}_{\alpha }}\sqrt{Var(\widehat{T})\text{ }}(\text{upper bound})\,\!$$


 * $${{T}_{L}}=\widehat{T}-{{K}_{\alpha }}\sqrt{Var(\widehat{T})\text{ }}(\text{lower bound})\,\!$$