The Logistic Distribution: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
No edit summary
 
(58 intermediate revisions by 7 users not shown)
Line 1: Line 1:
{{template:LDABOOK|14|The Logistic Distribution}}
{{template:LDABOOK|14|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 [[Appendix:_Life_Data_Analysis_References|[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.


===The Logistic Distribution===
===Logistic Probability Density Function===
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 [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.
The logistic ''pdf'' is given by:  
====Logistic Probability Density Function====
The logistic <math>pdf</math>  is given by:  


::<math>\begin{matrix}
::<math>\begin{align}
   f(T)=\tfrac{{{e}^{z}}}{\sigma {{(1+{{e}^{z}})}^{2}}}  \\
   f(t)=\tfrac{{{e}^{z}}}{\sigma {{(1+{{e}^{z}})}^{2}}}  \\ \\
   z=\tfrac{t-\mu }{\sigma }  \\
 
   -\infty <T<\infty ,\ \ -\infty <\mu <\infty ,\sigma >0  \\
   z=\tfrac{t-\mu }{\sigma }  \\ \\
\end{matrix}</math>
 
   -\infty <t<\infty ,\ \ -\infty <\mu <\infty ,\sigma >0  \\ \\
\end{align}\,\!</math>


where:   
where:   


::<math>\begin{align}
::<math>\begin{align}
   \mu = & \text{location parameter (also denoted as }\overline{T)} \\  
   \mu = & \text{location parameter (also denoted as }\overline{T}) \\  
   \sigma = & \text{scale parameter}   
   \sigma = & \text{scale parameter}   
\end{align}</math>
\end{align}\,\!</math>


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


[[Image:ldaLD10.1.gif|thumb|center|600px| ]]
[[Image:WB.14 logistic pd.png|center|400px| ]]


====The Logistic Standard Deviation====
===The Logistic Standard Deviation===
The standard deviation of the logistic distribution,   
The standard deviation of the logistic distribution,   
<math>{{\sigma }_{T}}</math> , is given by:  
<math>{{\sigma }_{T}}\,\!</math>, is given by:  


::<math>{{\sigma }_{T}}=\sigma \pi \frac{\sqrt{3}}{3}</math>
::<math>{{\sigma }_{T}}=\sigma \pi \frac{\sqrt{3}}{3}\,\!</math>


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


====The Logistic Reliability Function====
::<math>R(t)=\int_{t}^{\infty }f(t)dt\,\!</math>
The reliability for a mission of time  <math>T</math> , starting at age 0, for the logistic distribution is determined by:
 
 
::<math>R(T)=\int_{T}^{\infty }f(t)dt</math>


or:
or:


 
::<math>R(t)=\frac{1}{1+{{e}^{z}}}\,\!</math>
::<math>R(T)=\frac{1}{1+{{e}^{z}}}</math>
 


The unreliability function is:
The unreliability function is:


 
::<math>F=\frac{{{e}^{z}}}{1+{{e}^{z}}}\,\!</math>
::<math>F=\frac{{{e}^{z}}}{1+{{e}^{z}}}</math>


where:
where:


::<math>z=\frac{t-\mu }{\sigma }\,\!</math>


::<math>z=\frac{T-\mu }{\sigma }</math>
===The Logistic Conditional Reliability Function===
 
====The Logistic Conditional Reliability Function====
The logistic conditional reliability function is given by:  
The logistic conditional reliability function is given by:  


::<math>R(t/T)=\frac{R(T+t)}{R(T)}=\frac{1+{{e}^{\tfrac{T-\mu }{\sigma }}}}{1+{{e}^{\tfrac{t+T-\mu }{\sigma }}}}</math>
::<math>R(t|T)=\frac{R(T+t)}{R(T)}=\frac{1+{{e}^{\tfrac{T-\mu }{\sigma }}}}{1+{{e}^{\tfrac{t+T-\mu }{\sigma }}}}\,\!</math>
 


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


::<math>\begin{align}
{{T}_{R}}=\mu +\sigma[\ln (1-R)-\ln (R)]
\end{align}\,\!</math>


::<math>{{T}_{R}}=\mu +\sigma [\ln (1-R)-\ln (R)]</math>
===The Logistic Failure Rate Function===
 
====The Logistic Failure Rate Function====
The logistic failure rate function is given by:  
The logistic failure rate function is given by:  


::<math>\lambda (T)=\frac{{{e}^{z}}}{\sigma (1+{{e}^{z}})}</math>
::<math>\lambda (t)=\frac{{{e}^{z}}}{\sigma (1+{{e}^{z}})}\,\!</math>
 
 
====Characteristics of the Logistic Distribution====
 
• The logistic distribution has no shape parameter. This means that the logistic  <math>pdf</math>  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,  <math>\mu </math> , or the mean life or the  <math>MTTF</math> , is also the location parameter of the logistic  <math>pdf</math> , as it locates the  <math>pdf</math>  along the abscissa. It can assume values of  <math>-\infty <\bar{T}<\infty </math> .


As <math>\mu </math> decreases, the <math>pdf</math> is shifted to the left.
==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, <math>\mu \,\!</math>, or the mean life or the <math>MTTF\,\!</math>, is also the location parameter of the logistic ''pdf'', as it locates the  ''pdf'' along the abscissa. It can assume values of <math>-\infty <\bar{T}<\infty \,\!</math>.
:* As <math>\mu \,\!</math> decreases, the ''pdf'' is shifted to the left.
:* As <math>\mu \,\!</math> increases, the ''pdf'' is shifted to the right.


• As  <math>\mu </math>  increases, the  <math>pdf</math>  is shifted to the right.
[[Image:WB.14 effect of mu.png|center|400px| ]]


[[Image:ldaLD10.2.gif|thumb|center|600px| ]]
:* As <math>\sigma \,\!</math> decreases, the ''pdf'' gets pushed toward the mean, or it becomes narrower and taller.
:* As <math>\sigma \,\!</math> increases, the ''pdf'' spreads out away from the mean, or it becomes broader and shallower.
:* The scale parameter can assume values of <math>0<\sigma <\infty \,\!</math>.


• As  <math>\sigma </math>  decreases, the  <math>pdf</math>  gets pushed toward the mean, or it becomes narrower and taller.
[[Image:WB.14 effect of sigma.png|center|400px| ]]  
 
• As  <math>\sigma </math>  increases, the  <math>pdf</math>  spreads out away from the mean, or it becomes broader and shallower.
 
• The scale parameter can assume values of  <math>0<\sigma <\infty </math>.
 
[[Image:ldaLD10.3.gif|thumb|center|600px| ]]  
The logistic <math>pdf</math>  starts at <math>T=-\infty </math> with an <math>f(T)=0</math> . As <math>T</math> increases, <math>f(T)</math> also increases, goes through its point of inflection and reaches its maximum value at <math>T=\bar{T}</math> . Thereafter, <math>f(T)</math> decreases, goes through its point of inflection and assumes a value of <math>f(T)=0</math> at  <math>T=+\infty </math> .
:* The logistic ''pdf'' starts at <math>t=-\infty \,\!</math> with an <math>f(t)=0\,\!</math>. As <math>t\,\!</math> increases, <math>f(t)\,\!</math> also increases, goes through its point of inflection and reaches its maximum value at <math>t=\bar{T}\,\!</math>. Thereafter, <math>f(t)\,\!</math> decreases, goes through its point of inflection and assumes a value of <math>f(t)=0\,\!</math> at <math>t=+\infty \,\!</math>.
:* For <math>t=\pm \infty ,\,\!</math> the ''pdf'' equals 0. The maximum value of the ''pdf'' occurs at <math>t = \mu \,\!</math> and equals <math>\tfrac{1}{4\sigma }.\,\!</math>  
:* 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 <math>t=\mu +\sigma \ln (2\pm \sqrt{3})\,\!</math> or <math>t\approx \mu \pm \sigma 1.31696\,\!</math>.
:* If the location parameter <math>\mu \,\!</math> decreases, the reliability plot is shifted to the left. If <math>\mu \,\!</math> increases, the reliability plot is shifted to the right.
:* If <math>t=\mu \,\!</math> then <math>R=0.5\,\!</math>.    is the inflection point. If <math>t<\mu \,\!</math>  then <math>R(t)\,\!</math> is concave (concave down); if <math>t>\mu \,\!</math> then <math>R(t)\,\!</math> is convex (concave up). For <math>t<\mu ,\,\!</math> <math>\lambda (t)\,\!</math> is convex (concave up), for <math>t>\mu ;\,\!</math> <math>\lambda (t)\,\!</math> 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 <math>t\,\!</math> approaching 1/ <math>\delta .\,\!</math>
:* If location parameter <math>\mu \,\!</math> decreases, the failure rate plot is shifted to the left. Vice versa if <math>\mu \,\!</math> increases, the failure rate plot is shifted to the right.
:* <math>\lambda \,\!</math> always increases. For <math>t\to -\infty \,\!</math>    for <math>t\to \infty \,\!</math>    It is always <math>0\le \lambda (t)\le \tfrac{1}{\sigma }.\,\!</math>
:* If <math>\sigma \,\!</math> increases, then <math>\lambda (t)\,\!</math> increases more slowly and smoothly. The segment of time where <math>0<\lambda (t)<\tfrac{1}{\sigma }\,\!</math> increases, too, whereas the region where <math>\lambda (t)\,\!</math> is close to 0 or <math>\tfrac{1}{\sigma }\,\!</math> gets narrower. Conversely, if <math>\sigma \,\!</math> decreases, then <math>\lambda (t)\,\!</math> increases more quickly and sharply. The segment of time where <math>0 < \lambda (t)<\tfrac{1}{\sigma }\,\!</math> decreases, too, whereas the region where <math>\lambda (t)\,\!</math> is close to 0 or <math>\tfrac{1}{\sigma }\,\!</math> gets broader.


• For  <math>T=\pm \infty ,</math>  the  <math>pdf</math>  equals  <math>0.</math>  The maximum value of the  <math>pdf</math>  occurs at  <math>T</math> = <math>\mu </math>  and equals  <math>\tfrac{1}{4\sigma }.</math>


• The point of inflection of the  <math>pdf</math>  plot is the point where the second derivative of the  <math>pdf</math>  equals zero. The inflection point occurs at  <math>T=\mu +\sigma \ln (2\pm \sqrt{3})</math>  or  <math>T\approx \mu \pm \sigma 1.31696</math>.
'''Weibull++ Notes on Negative Time Values'''


• If the location parameter  <math>\mu </math>  decreases, the reliability plot is shifted to the left. If  <math>\mu </math>  increases, the reliability plot is shifted to the right.
• If  <math>T=\mu </math>  then  <math>R=0.5</math> .    is the inflection point. If  <math>T<\mu </math>  then  <math>R(t)</math>  is concave (concave down); if  <math>T>\mu </math>  then  <math>R(t)</math>  is convex (concave up). For  <math>T<\mu ,</math>  <math>\lambda (t)</math>  is convex (concave up), for  <math>T>\mu ;</math>  <math>\lambda (t)</math>  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  <math>t</math>  approaching 1/ <math>\delta .</math>
• If location parameter  <math>\mu </math>  decreases, the failure rate plot is shifted to the left. Vice versa if  <math>\mu </math>  increases, the failure rate plot is shifted to the right.
• <math>\lambda </math>  always increases. For  <math>T\to -\infty </math>      for  <math>T\to \infty </math>      It is always  <math>0\le \lambda (t)\le \tfrac{1}{\sigma }.</math>
• If  <math>\sigma </math>  increases, then  <math>\lambda (t)</math>  increases more slowly and smoothly. The segment of time where  <math>0<\lambda (t)<\tfrac{1}{\sigma }</math>  increases, too, whereas the region where  <math>\lambda (t)</math>  is close to  <math>0</math>  or  <math>\tfrac{1}{\sigma }</math>  gets narrower. Conversely, if  <math>\sigma </math>  decreases, then  <math>\lambda (t)</math>  increases more quickly and sharply. The segment of time where  <math>0<</math>  <math>\lambda (t)<\tfrac{1}{\sigma }</math>  decreases, too, whereas the region where  <math>\lambda (t)</math>  is close to  <math>0</math>  or  <math>\tfrac{1}{\sigma }</math>  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).
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, <math>z\,\!</math> can be calculated as a function of the ''cdf'' <math>F\,\!</math> as follows:


====Probability Paper====
::<math>\begin{align}
The form of the Logistic probability paper is based on linearizing the  <math>cdf</math> .
z=\ln (F)-\ln (1-F)
From Eqn. (UnR fcn),  <math>z</math>  can be calculated as a function of the  <math>cdf</math>  <math>F</math>  as follows:
\end{align}\,\!</math>
 
::<math>z=\ln (F)-\ln (1-F)</math>
 


or using Eqn. (z func of parameters)
or using the equation for <math>z\,\!</math>


::<math>\frac{T-\mu }{\sigma }=\ln (F)-\ln (1-F)</math>
::<math>\frac{t-\mu }{\sigma }=\ln (F)-\ln (1-F)\,\!</math>


Then:  
Then:  


::<math>\ln (F)-\ln (1-F)=-\frac{\mu }{\sigma }+\frac{1}{\sigma }T</math>
::<math>\ln (F)-\ln (1-F)=-\frac{\mu }{\sigma }+\frac{1}{\sigma }t\,\!</math>
 


Now let:  
Now let:  


::<math>y=\ln (F)-\ln (1-F)</math>
::<math>\begin{align}
 
y=\ln (F)-\ln (1-F)
\end{align}\,\!</math>
 
::<math>x=T</math>


::<math>\begin{align}
x=t
\end{align}\,\!</math>


and:  
and:  


::<math>a=-\frac{\mu }{\sigma }</math>
::<math>a=-\frac{\mu }{\sigma }\,\!</math>
 
 
::<math>b=\frac{1}{\sigma }</math>


::<math>b=\frac{1}{\sigma }\,\!</math>


which results in the following linear equation:  
which results in the following linear equation:  


::<math>y=a+bx</math>
::<math>\begin{align}
 
y=a+bx
 
\end{align}\,\!</math>
The logistic probability paper resulting from this linearized  <math>cdf</math>  function is shown next.
 
[[Image:ldaLD10.4.gif|thumb|center|600px| ]]
 
Since the logistic distribution is symmetrical, the area under the  <math>pdf</math>  curve from  <math>-\infty </math>  to  <math>\mu </math>  is  <math>0.5</math> , as is the area from  <math>\mu </math>  to  <math>+\infty </math> . Consequently, the value of  <math>\mu </math>  is said to be the point where  <math>R(t)=Q(t)=50%</math> .  This means that the estimate of  <math>\mu </math>  can be read from the point where the plotted line crosses the 50% unreliability line.
 
For  <math>z=1</math> ,  <math>\sigma =t-\mu </math>  and  <math>R(t)=\tfrac{1}{1+\exp (1)}\approx 0.2689.</math>  Therefore, <math>\sigma </math>  can be found by subtracting  <math>\mu </math> 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 Chapter 5.


=====Bounds on the Parameters=====
The logistic probability paper resulting from this linearized ''cdf'' function is shown next.
The lower and upper bounds on the location parameter  <math>\widehat{\mu }</math>  are estimated from
:


::<math>{{\mu }_{U}}=\widehat{\mu }+{{K}_{\alpha }}\sqrt{Var(\widehat{\mu })\text{ }}\text{ (upper bound)}</math>
[[Image:WB.14 logistic probability plot.png|center|400px| ]]


<math>{{\mu }_{L}}=\widehat{\mu }-{{K}_{\alpha }}\sqrt{Var(\widehat{\mu })\text{ }}\text{ (lower bound)}</math>
Since the logistic distribution is symmetrical, the area under the ''pdf'' curve from <math>-\infty \,\!</math> to <math>\mu \,\!</math> is 0.5, as is the area from <math>\mu \,\!</math> to <math>+\infty \,\!</math>. Consequently, the value of <math>\mu \,\!</math> is said to be the point where <math>R(t)=Q(t)=50%\,\!</math>. This means that the estimate of <math>\mu \,\!</math> can be read from the point where the plotted line crosses the 50% unreliability line.
 
The lower and upper bounds on the scale parameter  <math>\widehat{\sigma }</math>  are estimated from:
 
::<math>{{\sigma }_{U}}=\widehat{\sigma }{{e}^{\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{\sigma })\text{ }}}{\widehat{\sigma }}}}(\text{upper bound})</math>
 
 
::<math>{{\sigma }_{L}}=\widehat{\sigma }{{e}^{\tfrac{-{{K}_{\alpha }}\sqrt{Var(\widehat{\sigma })\text{ }}}{\widehat{\sigma }}}}\text{ (lower bound)}</math>
 
where  <math>{{K}_{\alpha }}</math>  is defined by:
 
::<math>\alpha =\frac{1}{\sqrt{2\pi }}\int_{{{K}_{\alpha }}}^{\infty }{{e}^{-\tfrac{{{t}^{2}}}{2}}}dt=1-\Phi ({{K}_{\alpha }})</math>
 
 
If  <math>\delta </math>  is the confidence level, then  <math>\alpha =\tfrac{1-\delta }{2}</math>  for the two-sided bounds, and  <math>\alpha =1-\delta </math>  for the one-sided bounds.
 
The variances and covariances of  <math>\widehat{\mu }</math>  and  <math>\widehat{\sigma }</math>  are estimated from the Fisher matrix, as follows:
 
::<math>\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}</math>
 
 
<math>\Lambda </math>  is the log-likelihood function of the normal distribution, described in Chapter 3 and Appendix C.
 
=====Bounds on Reliability=====
The reliability of the logistic distribution is:
 
::<math>\widehat{R}=\frac{1}{1+{{e}^{\widehat{z}}}}</math>
 
where:
 
::<math>\widehat{z}=\frac{T-\widehat{\mu }}{\widehat{\sigma }}</math>
 
 
Here  <math>-\infty <T<\infty </math> , <math>-\infty <\mu <\infty </math> ,  <math>0<\sigma <\infty </math> . Therefore, <math>z</math>  also is changing from  <math>-\infty </math>  to  <math>+\infty </math> . Then the bounds on  <math>z</math>  are estimated from:
 
::<math>{{z}_{U}}=\widehat{z}+{{K}_{\alpha }}\sqrt{Var(\widehat{z})\text{ }}</math>
 
 
::<math>{{z}_{L}}=\widehat{z}-{{K}_{\alpha }}\sqrt{Var(\widehat{z})\text{ }}\text{ }</math>
 
 
where:
 
::<math>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 })</math>
 
or:
 
::<math>Var(\widehat{z})=\frac{1}{{{\sigma }^{2}}}(Var(\widehat{\mu })+2\widehat{z}Cov(\widehat{\mu },\widehat{\sigma })+{{\widehat{z}}^{2}}Var(\widehat{\sigma }))</math>
 
The upper and lower bounds on reliability are:
 
::<math>{{R}_{U}}=\frac{1}{1+{{e}^{{{z}_{L}}}}}\text{(upper bound)}</math>
 
::<math>{{R}_{L}}=\frac{1}{1+{{e}^{{{z}_{U}}}}}\text{(lower bound)}</math>
 
=====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:
 
::<math>\widehat{T}(\widehat{\mu },\widehat{\sigma })=\widehat{\mu }+\widehat{\sigma }z</math>
 
 
where:
 
 
::<math>z=\ln (1-R)-\ln (R)</math>
 
 
 
 
::<math>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 })</math>
 
 
or:
 
 
::<math>Var(\widehat{T})=Var(\widehat{\mu })+2\widehat{z}Cov(\widehat{\mu },\widehat{\sigma })+{{\widehat{z}}^{2}}Var(\widehat{\sigma })</math>
 
 
The upper and lower bounds are then found by:
 
::<math>{{T}_{U}}=\widehat{T}+{{K}_{\alpha }}\sqrt{Var(\widehat{T})\text{ }}(\text{upper bound})</math>
 
   
 
::<math>{{T}_{L}}=\widehat{T}-{{K}_{\alpha }}\sqrt{Var(\widehat{T})\text{ }}(\text{lower bound})</math>
 
=====A Logistic Distribution Example=====
The lifetime of a mechanical valve is known to follow a logistic distribution. Ten units were tested for 28 months and the following months-to-failure data was collected.
 
 
::<math>\overset{{}}{\mathop{\text{Table 10}\text{.1 - Times-to-Failure Data with Suspensions}}}\,</math>
 
 
::<math>\begin{matrix}
  \text{Data Point Index} & \text{State F or S} & \text{State End Time}  \\
  \text{1} & \text{F} & \text{8}  \\
  \text{2} & \text{F} & \text{10}  \\
  \text{3} & \text{F} & \text{15}  \\
  \text{4} & \text{F} & \text{17}  \\
  \text{5} & \text{F} & \text{19}  \\
  \text{6} & \text{F} & \text{26}  \\
  \text{7} & \text{F} & \text{27}  \\
  \text{8} & \text{S} & \text{28}  \\
  \text{9} & \text{S} & \text{28}  \\
  \text{10} & \text{S} & \text{28}  \\
\end{matrix}</math>
 
• Determine the valve's design life if specifications call for a reliability goal of 0.90.
 
• The valve is to be used in a pumping device that requires 1 month of continuous operation. What is the probability of the pump failing due to the valve?
 
This data set can be entered into Weibull++ as follows:
 
[[File:ldaLD10.5.gif|center]]
 
The computed parameters for maximum likelihood are:
 
::<math>\begin{align}
  & \widehat{\mu }= & 22.34 \\
& \hat{\sigma }= & 6.15 
\end{align}</math>


• The valve's design life, along with 90% two sided confidence bounds, can be obtained using the QCP as follows:
For  <math>z=1\,\!</math>, <math>\sigma =t-\mu \,\!</math> and <math>R(t)=\tfrac{1}{1+\exp (1)}\approx 0.2689.\,\!</math> Therefore, <math>\sigma \,\!</math> can be found by subtracting <math>\mu \,\!</math> from the time value where the plotted probability line crosses the 73.10% unreliability (26.89% reliability) horizontal line.


[[File:ldaLD10.6.gif|center]]


• The probability, along with 90% two sided confidence bounds, that the pump fails due to a valve failure during the first month is obtained as follows:


==General Example==


[[File:ldaLD10.7.gif|center]]
{{:Logistic Distribution Example}}

Latest revision as of 00:59, 10 March 2023

New format available! This reference is now available in a new format that offers faster page load, improved display for calculations and images, more targeted search and the latest content available as a PDF. As of September 2023, this Reliawiki page will not continue to be updated. Please update all links and bookmarks to the latest reference at help.reliasoft.com/reference/life_data_analysis

Chapter 14: The Logistic Distribution


Weibullbox.png

Chapter 14  
The Logistic Distribution  

Synthesis-icon.png

Available Software:
Weibull++

Examples icon.png

More Resources:
Weibull++ Examples Collection

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:

[math]\displaystyle{ \begin{align} f(t)=\tfrac{{{e}^{z}}}{\sigma {{(1+{{e}^{z}})}^{2}}} \\ \\ z=\tfrac{t-\mu }{\sigma } \\ \\ -\infty \lt t\lt \infty ,\ \ -\infty \lt \mu \lt \infty ,\sigma \gt 0 \\ \\ \end{align}\,\! }[/math]

where:

[math]\displaystyle{ \begin{align} \mu = & \text{location parameter (also denoted as }\overline{T}) \\ \sigma = & \text{scale parameter} \end{align}\,\! }[/math]

The Logistic Mean, Median and Mode

The logistic mean or MTTF is actually one of the parameters of the distribution, usually denoted as [math]\displaystyle{ \mu \,\! }[/math]. Since the logistic distribution is symmetrical, the median and the mode are always equal to the mean, [math]\displaystyle{ \mu =\tilde{T}=\breve{T}.\,\! }[/math]

WB.14 logistic pd.png

The Logistic Standard Deviation

The standard deviation of the logistic distribution, [math]\displaystyle{ {{\sigma }_{T}}\,\! }[/math], is given by:

[math]\displaystyle{ {{\sigma }_{T}}=\sigma \pi \frac{\sqrt{3}}{3}\,\! }[/math]

The Logistic Reliability Function

The reliability for a mission of time [math]\displaystyle{ t\,\! }[/math], starting at age 0, for the logistic distribution is determined by:

[math]\displaystyle{ R(t)=\int_{t}^{\infty }f(t)dt\,\! }[/math]

or:

[math]\displaystyle{ R(t)=\frac{1}{1+{{e}^{z}}}\,\! }[/math]

The unreliability function is:

[math]\displaystyle{ F=\frac{{{e}^{z}}}{1+{{e}^{z}}}\,\! }[/math]

where:

[math]\displaystyle{ z=\frac{t-\mu }{\sigma }\,\! }[/math]

The Logistic Conditional Reliability Function

The logistic conditional reliability function is given by:

[math]\displaystyle{ R(t|T)=\frac{R(T+t)}{R(T)}=\frac{1+{{e}^{\tfrac{T-\mu }{\sigma }}}}{1+{{e}^{\tfrac{t+T-\mu }{\sigma }}}}\,\! }[/math]

The Logistic Reliable Life

The logistic reliable life is given by:

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

The Logistic Failure Rate Function

The logistic failure rate function is given by:

[math]\displaystyle{ \lambda (t)=\frac{{{e}^{z}}}{\sigma (1+{{e}^{z}})}\,\! }[/math]

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, [math]\displaystyle{ \mu \,\! }[/math], or the mean life or the [math]\displaystyle{ MTTF\,\! }[/math], is also the location parameter of the logistic pdf, as it locates the pdf along the abscissa. It can assume values of [math]\displaystyle{ -\infty \lt \bar{T}\lt \infty \,\! }[/math].
  • As [math]\displaystyle{ \mu \,\! }[/math] decreases, the pdf is shifted to the left.
  • As [math]\displaystyle{ \mu \,\! }[/math] increases, the pdf is shifted to the right.
WB.14 effect of mu.png
  • As [math]\displaystyle{ \sigma \,\! }[/math] decreases, the pdf gets pushed toward the mean, or it becomes narrower and taller.
  • As [math]\displaystyle{ \sigma \,\! }[/math] increases, the pdf spreads out away from the mean, or it becomes broader and shallower.
  • The scale parameter can assume values of [math]\displaystyle{ 0\lt \sigma \lt \infty \,\! }[/math].
WB.14 effect of sigma.png
  • The logistic pdf starts at [math]\displaystyle{ t=-\infty \,\! }[/math] with an [math]\displaystyle{ f(t)=0\,\! }[/math]. As [math]\displaystyle{ t\,\! }[/math] increases, [math]\displaystyle{ f(t)\,\! }[/math] also increases, goes through its point of inflection and reaches its maximum value at [math]\displaystyle{ t=\bar{T}\,\! }[/math]. Thereafter, [math]\displaystyle{ f(t)\,\! }[/math] decreases, goes through its point of inflection and assumes a value of [math]\displaystyle{ f(t)=0\,\! }[/math] at [math]\displaystyle{ t=+\infty \,\! }[/math].
  • For [math]\displaystyle{ t=\pm \infty ,\,\! }[/math] the pdf equals 0. The maximum value of the pdf occurs at [math]\displaystyle{ t = \mu \,\! }[/math] and equals [math]\displaystyle{ \tfrac{1}{4\sigma }.\,\! }[/math]
  • 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 [math]\displaystyle{ t=\mu +\sigma \ln (2\pm \sqrt{3})\,\! }[/math] or [math]\displaystyle{ t\approx \mu \pm \sigma 1.31696\,\! }[/math].
  • If the location parameter [math]\displaystyle{ \mu \,\! }[/math] decreases, the reliability plot is shifted to the left. If [math]\displaystyle{ \mu \,\! }[/math] increases, the reliability plot is shifted to the right.
  • If [math]\displaystyle{ t=\mu \,\! }[/math] then [math]\displaystyle{ R=0.5\,\! }[/math]. is the inflection point. If [math]\displaystyle{ t\lt \mu \,\! }[/math] then [math]\displaystyle{ R(t)\,\! }[/math] is concave (concave down); if [math]\displaystyle{ t\gt \mu \,\! }[/math] then [math]\displaystyle{ R(t)\,\! }[/math] is convex (concave up). For [math]\displaystyle{ t\lt \mu ,\,\! }[/math] [math]\displaystyle{ \lambda (t)\,\! }[/math] is convex (concave up), for [math]\displaystyle{ t\gt \mu ;\,\! }[/math] [math]\displaystyle{ \lambda (t)\,\! }[/math] 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 [math]\displaystyle{ t\,\! }[/math] approaching 1/ [math]\displaystyle{ \delta .\,\! }[/math]
  • If location parameter [math]\displaystyle{ \mu \,\! }[/math] decreases, the failure rate plot is shifted to the left. Vice versa if [math]\displaystyle{ \mu \,\! }[/math] increases, the failure rate plot is shifted to the right.
  • [math]\displaystyle{ \lambda \,\! }[/math] always increases. For [math]\displaystyle{ t\to -\infty \,\! }[/math] for [math]\displaystyle{ t\to \infty \,\! }[/math] It is always [math]\displaystyle{ 0\le \lambda (t)\le \tfrac{1}{\sigma }.\,\! }[/math]
  • If [math]\displaystyle{ \sigma \,\! }[/math] increases, then [math]\displaystyle{ \lambda (t)\,\! }[/math] increases more slowly and smoothly. The segment of time where [math]\displaystyle{ 0\lt \lambda (t)\lt \tfrac{1}{\sigma }\,\! }[/math] increases, too, whereas the region where [math]\displaystyle{ \lambda (t)\,\! }[/math] is close to 0 or [math]\displaystyle{ \tfrac{1}{\sigma }\,\! }[/math] gets narrower. Conversely, if [math]\displaystyle{ \sigma \,\! }[/math] decreases, then [math]\displaystyle{ \lambda (t)\,\! }[/math] increases more quickly and sharply. The segment of time where [math]\displaystyle{ 0 \lt \lambda (t)\lt \tfrac{1}{\sigma }\,\! }[/math] decreases, too, whereas the region where [math]\displaystyle{ \lambda (t)\,\! }[/math] is close to 0 or [math]\displaystyle{ \tfrac{1}{\sigma }\,\! }[/math] 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, [math]\displaystyle{ z\,\! }[/math] can be calculated as a function of the cdf [math]\displaystyle{ F\,\! }[/math] as follows:

[math]\displaystyle{ \begin{align} z=\ln (F)-\ln (1-F) \end{align}\,\! }[/math]

or using the equation for [math]\displaystyle{ z\,\! }[/math]

[math]\displaystyle{ \frac{t-\mu }{\sigma }=\ln (F)-\ln (1-F)\,\! }[/math]

Then:

[math]\displaystyle{ \ln (F)-\ln (1-F)=-\frac{\mu }{\sigma }+\frac{1}{\sigma }t\,\! }[/math]

Now let:

[math]\displaystyle{ \begin{align} y=\ln (F)-\ln (1-F) \end{align}\,\! }[/math]
[math]\displaystyle{ \begin{align} x=t \end{align}\,\! }[/math]

and:

[math]\displaystyle{ a=-\frac{\mu }{\sigma }\,\! }[/math]
[math]\displaystyle{ b=\frac{1}{\sigma }\,\! }[/math]

which results in the following linear equation:

[math]\displaystyle{ \begin{align} y=a+bx \end{align}\,\! }[/math]

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

WB.14 logistic probability plot.png

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

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


General Example

The lifetime of a mechanical valve is known to follow a logistic distribution. 10 units were tested for 28 months and the following months-to-failure data were collected.

Data Point Index State F or S State End Time
1 F 8
2 F 10
3 F 15
4 F 17
5 F 19
6 F 26
7 F 27
8 S 28
9 S 28
10 S 28


  • Determine the valve's design life if specifications call for a reliability goal of 0.90.
  • The valve is to be used in a pumping device that requires 1 month of continuous operation. What is the probability of the pump failing due to the valve?

Enter the data set in a Weibull++ standard folio, as follows:

Logistic Distribution Exmaple 1 Data.png

The computed parameters for maximum likelihood are:

[math]\displaystyle{ \begin{align} & \widehat{\mu }= & 22.34 \\ & \hat{\sigma }= & 6.15 \end{align}\,\! }[/math]

The valve's design life, along with 90% two sided confidence bounds, can be obtained using the QCP as follows:

Logistic Distribution Exmaple 1 QCP Reliable Life.png

The probability, along with 90% two sided confidence bounds, that the pump fails due to a valve failure during the first month is obtained as follows:

Logistic Distribution Exmaple 1 QCP Reliability.png