|
|
(11 intermediate revisions by the same user not shown) |
Line 1: |
Line 1: |
| {{template:ALTABOOK|3}} | | {{template:ALTABOOK|3}} |
| In this chapter, we will briefly present three lifetime distributions commonly used in accelerated life test analysis: the 1-parameter exponential, the 2-parameter Weibull and the lognormal distributions. Readers who are interested in a more rigorous overview (or in different forms of these and other life distributions) can refer to ReliaSoft's [[Life Data Analysis Reference Book|Life Data Analysis Reference]]. For information about the parameter estimation methods, see [[Appendix_B:_Parameter_Estimation|Appendix B]].<br> | | In this chapter, we will briefly present three lifetime distributions commonly used in accelerated life test analysis: the exponential, the Weibull and the lognormal distributions. Note that although all forms are mentioned below, ALTA uses the 1-parameter form of the exponential distribution and the 2-parameter form of the Weibull distribution. |
| | |
| | Readers who are interested in a more rigorous overview of these distributions (or for information about other life distributions) can refer to the [https://help.reliasoft.com/reference/life_data_analysis Life data analysis reference]. For information about the parameter estimation methods, see [[Appendix_B:_Parameter_Estimation|Appendix B]].<br> |
|
| |
|
| =The Exponential Distribution= | | =The Exponential Distribution= |
| | | {{Exponential Distribution Definition}} |
| The exponential distribution is a very commonly used distribution in reliability engineering. Due to its simplicity, it has been widely employed, even in cases where it doesn't apply. The exponential distribution is used to describe units that have a constant failure rate. The 1-parameter exponential ''pdf ''is given by:
| |
| | |
| <br>
| |
| | |
| ::<math>\begin{align}
| |
| & f(T)= \lambda {{e}^{-\lambda T}}=\frac{1}{m}{{e}^{-\tfrac{1}{m}T}} \
| |
| & T\ge 0,\lambda >0,m>0
| |
| \end{align}</math>
| |
| | |
| <br> where:
| |
| | |
| ::*<span class="texhtml">λ = </span> constant failure rate, in failures per unit of measurement (e.g. failures per hour, per cycle, etc.).
| |
| ::*<math>\lambda =\tfrac{1}{m}</math>.
| |
| ::*<span class="texhtml">''m'' = </span> mean time between failures, or to a failure.
| |
| ::*<span class="texhtml">''T'' = </span> operating time, life, or age, in hours, cycles, miles, actuations, etc. This distribution requires the estimation of only one parameter, <span class="texhtml">λ</span> , for its application.
| |
| <br>
| |
| ==Exponential Distribution Functions== | | ==Exponential Distribution Functions== |
| {{:Exponential Distribution Functions}} | | {{:Exponential Distribution Functions}} |
| | | == Characteristics of the Exponential Distribution== |
| | | {{:Exponential Distribution Characteristics}} |
| == Characteristics == | |
| | |
| The characteristics of the 1-parameter exponential distribution can be exemplified by examining its parameter, lambda (<span class="texhtml">λ<span style="font-family: sans-serif;">)</span></span> and the effect lambda has on the ''pdf'', reliability and failure rate functions.<br>
| |
| | |
| ==== Effects of <span class="texhtml">λ</span> on the ''pdf'' ====
| |
| | |
| <br>
| |
| | |
| [[Image:ALTA4pdf.png|center|400px]] <br>
| |
| | |
| :*The scale parameter is <math>\tfrac{1}{\lambda }</math> .
| |
| :*As <span class="texhtml">λ</span> is decreased in value, the distribution is stretched out to the right, and as <span class="texhtml">λ</span> is increased, the distribution is pushed toward the origin. | |
| :*This distribution has no shape parameter, as it has only one shape (i.e., the exponential). Its only parameter is the failure rate, <span class="texhtml">λ</span> .
| |
| :*The distribution starts at <span class="texhtml">''T'' = 0</span> at the level of <span class="texhtml">''f''(''T'' = 0) = λ</span>. It decreases thereafter exponentially and monotonically as <span class="texhtml">''T''</span> increases, and it is convex.
| |
| :*As <math>T\to \infty </math> , <math>f(T)\to 0</math>.
| |
| :*This ''pdf'' can be thought of as a special case of the Weibull ''pdf'' with <span class="texhtml">β = 1</span>.
| |
| | |
| <br> [[Image:ALTA4reliabilityvstimeplot.png|center|400px]] <br>
| |
| | |
| ==== Effects of <span class="texhtml">λ</span> on the Reliability Function ====
| |
| | |
| :*The 1-parameter exponential reliability function starts at the value of 1 at <span class="texhtml">''T'' = 0</span> . It decreases thereafter monotonically and is convex.
| |
| :*As <math>T\to \infty </math> , <math>R(T\to \infty )\to 0</math>.
| |
| | |
| ==== Effects of <span class="texhtml">λ</span> on the Failure Rate Function ====
| |
| | |
| The failure rate function for the exponential distribution is constant and equal to the parameter <span class="texhtml">λ</span> .
| |
| | |
| <br> [[Image:ALTA4FRvsTP.png|center|400px]] <br>
| |
| | |
| | |
| == Parameter Estimation ==
| |
| The parameter of the exponential distribution can be estimated graphically on probability plotting paper or analytically using either least squares or maximum likelihood. (Parameter estimation methods are presented in detail in [[Appendix B: Parameter Estimation|Appendix B]].)
| |
| | |
| ==== Probability Plotting ====
| |
| One method of calculating the parameter of the exponential distribution is by using probability plotting. To better illustrate this procedure, consider the following example.<br>
| |
| | |
| {{:Probability_Plotting_Example}}
| |
| | |
| ==== MLE Parameter Estimation ====
| |
| | |
| The parameter of the exponential distribution can also be estimated using the maximum likelihood estimation (MLE) method. This log-likelihood function is: <br>
| |
| | |
| ::<math>\ln (L)=\Lambda =\underset{i=1}{\overset{{{F}_{e}}}{\mathop \sum }}\,{{N}_{i}}\ln \left[ \lambda {{e}^{-\lambda {{T}_{i}}}} \right]-\underset{i=1}{\overset{S}{\mathop \sum }}\,N_{i}^{\prime }\lambda T_{i}^{\prime }+\overset{FI}{\mathop{\underset{i=1}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{\prime \prime }\ln [R_{Li}^{\prime \prime }-R_{Ri}^{\prime \prime }]</math>
| |
| | |
| :where:
| |
| | |
| <br>
| |
| | |
| ::<math>R_{Li}^{\prime \prime }={{e}^{-\lambda T_{Li}^{\prime \prime }}}</math>
| |
| | |
| <br>
| |
| | |
| ::<math>R_{Ri}^{\prime \prime }={{e}^{-\lambda T_{Ri}^{\prime \prime }}}</math> <br>
| |
| | |
| <br>
| |
| | |
| :and:
| |
| | |
| <br>
| |
| | |
| :*<span class="texhtml">''F''<sub>''e''</sub></span> is the number of groups of times-to-failure data points.
| |
| :*<span class="texhtml">''N''<sub>''i''</sub></span> is the number of times-to-failure in the <span class="texhtml">''i''<sup>''th''</sup></span> time-to-failure data group.
| |
| :*<span class="texhtml">λ</span> is the failure rate parameter (unknown a priori, the only parameter to be found).
| |
| :*<span class="texhtml">''T''<sub>''i''</sub></span> is the time of the <span class="texhtml">''i''<sup>''th''</sup></span> group of time-to-failure data.
| |
| :*<span class="texhtml">''S''</span> is the number of groups of suspension data points.
| |
| :*<math>N_{i}^{\prime }</math> is the number of suspensions in the <span class="texhtml">''i''<sup>''th''</sup></span> group of suspension data points.
| |
| :*<math>T_{i}^{\prime }</math> is the time of the <span class="texhtml">''i''<sup>''t'''h '''''</sup></span> suspension data group.
| |
| :*<span class="texhtml">''FI''</span> is the number of interval data groups.
| |
| :*<math>N_{i}^{\prime \prime }</math> is the number of intervals in the <span class="texhtml">''i''<sup>''th''</sup></span> group of data intervals.
| |
| :*<math>T_{Li}^{\prime \prime }</math> is the beginning of the <span class="texhtml">''i''<sup>''th''</sup></span> interval.
| |
| :*<math>T_{Ri}^{\prime \prime }</math> is the ending of the <span class="texhtml">''i''<sup>''th''</sup></span> interval.
| |
| | |
| <br> <br> The solution will be found by solving for a parameter <math>\widehat{\lambda }</math> so that <math>\tfrac{\partial \Lambda }{\partial \lambda }=0</math> where: <br>
| |
| | |
| ::<math>\frac{\partial \Lambda }{\partial \lambda }=\underset{i=1}{\overset{{{F}_{e}}}{\mathop \sum }}\,{{N}_{i}}\left( \frac{1}{\lambda }-{{T}_{i}} \right)-\underset{i=1}{\overset{S}{\mathop \sum }}\,N_{i}^{\prime }T_{i}^{\prime }-\overset{FI}{\mathop{\underset{i=1}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{\prime \prime }\frac{T_{Li}^{\prime \prime }R_{Li}^{\prime \prime }-T_{Ri}^{\prime \prime }R_{Ri}^{\prime \prime }}{R_{Li}^{\prime \prime }-R_{Ri}^{\prime \prime }}</math>
| |
| | |
| <br>
| |
| <br>
| |
| '''Example 2'''
| |
| <br>
| |
| Using the same data as in the probability plotting example (Example 1), and assuming an exponential distribution, estimate the parameter using the MLE method. <br>
| |
| <br>
| |
| '''Solution'''
| |
| <br>
| |
| In this example we have non-grouped data without suspensions. Therefore: <br>
| |
| | |
| ::<math>\frac{\partial \Lambda }{\partial \lambda }=\underset{i=1}{\overset{{{F}_{e}}}{\mathop \sum }}\,\left[ \frac{1}{\lambda }-\left( {{T}_{i}} \right) \right]=\underset{i=1}{\overset{14}{\mathop \sum }}\,\left[ \frac{1}{\lambda }-\left( {{T}_{i}} \right) \right]=0</math>
| |
| | |
| <br> Substituting the values for <span class="texhtml">''T''</span> we get: <br>
| |
| | |
| ::<math>\begin{align}
| |
| \frac{6}{\lambda }&= 4409,\text{ or:} \
| |
| \lambda& = 0.00136\text{ failure/hr}
| |
| \end{align}</math>
| |
| | |
| <br>
| |
|
| |
|
| = The Weibull Distribution = | | = The Weibull Distribution = |
| | | {{Weibull Distribution Definition}} |
| The Weibull distribution is one of the most commonly used distributions in reliability engineering because of the many shapes it attains for various values of <span class="texhtml">β</span> (slope). It can therefore model a great variety of data and life characteristics [[Reference Appendix D: References|[18]]].
| | ==Weibull Distribution Functions== |
| | | {{:Weibull Distribution Functions}} |
| The 2-parameter Weibull ''pdf'' is given by:
| | ==Characteristics of the Weibull Distribution== |
| | | {{:Weibull Distribution Characteristics}} |
| <br>
| |
| | |
| ::<math>f(T)=\frac{\beta }{\eta }{{\left( \frac{T}{\eta } \right)}^{\beta -1}}{{e}^{-{{\left( \tfrac{T}{\eta } \right)}^{\beta }}}}</math>
| |
| | |
| <br>
| |
| | |
| where:
| |
| | |
| <br>
| |
| | |
| ::<math>f(T)\ge 0,\text{ }T\ge 0,\text{ }\beta >0,\text{ }\eta >0\text{ }</math>
| |
| | |
| <br> and:
| |
| | |
| <br>
| |
| | |
| :*<span class="texhtml">η = </span> scale parameter.
| |
| :*<span class="texhtml">β = </span> shape parameter (or slope).
| |
| | |
| <br>
| |
| ==Statistical Properties Summary==
| |
| ==== The Mean or MTTF ====
| |
| | |
| The mean, <math>\overline{T}</math>, of the 2-parameter Weibull <span class="texhtml">''pdf''</span> is given by:
| |
| | |
| <br>
| |
| | |
| ::<math>\overline{T}=\eta \cdot \Gamma \left( \frac{1}{\beta }+1 \right)</math>
| |
| | |
| <br> where <math>\Gamma \left( \tfrac{1}{\beta }+1 \right)</math> is the gamma function evaluated at the value of <math>\left( \tfrac{1}{\beta }+1 \right)</math> .
| |
| | |
| <br>
| |
| | |
| ==== The Median ====
| |
| | |
| The median, <math>\breve{T}</math>, of the 2-parameter Weibull is given by:
| |
| | |
| <br>
| |
| | |
| ::<math>\breve{T}=\eta {{\left( \ln 2 \right)}^{\tfrac{1}{\beta }}}</math>
| |
| | |
| <br>
| |
| | |
| | |
| ==== The Mode ====
| |
| | |
| The mode, <math>\tilde{T}</math>, of the 2-parameter Weibull is given by:
| |
| | |
| <br>
| |
| | |
| ::<math>\tilde{T}=\eta {{\left( 1-\frac{1}{\beta } \right)}^{\tfrac{1}{\beta }}}</math>
| |
| | |
| <br>
| |
| | |
| | |
| ==== The Standard Deviation ====
| |
| | |
| The standard deviation, <span class="texhtml">σ<sub>''T''</sub></span>, of the 2-parameter Weibull is given by:
| |
| | |
| <br>
| |
| | |
| ::<math>{{\sigma }_{T}}=\eta \cdot \sqrt{\Gamma \left( \frac{2}{\beta }+1 \right)-\Gamma {{\left( \frac{1}{\beta }+1 \right)}^{2}}}</math>
| |
| | |
| <br>
| |
| | |
| | |
| ==== The <span class="texhtml">''cdf''</span> and the Reliability Function ====
| |
| | |
| The <span class="texhtml">''cdf''</span> of the 2-parameter Weibull distribution is given by:<br>
| |
| | |
| ::<math>F(T)=1-{{e}^{-{{\left( \tfrac{T}{\eta } \right)}^{\beta }}}}</math>
| |
| | |
| <br> The Weibull reliability function is given by:<br>
| |
| | |
| ::<math>\begin{align}
| |
| R(T)&= & 1-F(t) = \ {{e}^{-{{\left( \tfrac{T}{\eta } \right)}^{\beta }}}}
| |
| \end{align}</math>
| |
| | |
| <br>
| |
| | |
| | |
| ====The Conditional Reliability Function==== | |
| The Weibull conditional reliability function is given by:
| |
| | |
| | |
| ::<math>R(T,t)=\frac{R(T+t)}{R(T)}=\frac{{{e}^{-{{\left( \tfrac{T+t}{\eta } \right)}^{\beta }}}}}{{{e}^{-{{\left( \tfrac{T}{\eta } \right)}^{\beta }}}}}</math>
| |
| | |
| | |
| or:
| |
| | |
| ::<math>R(T,t)={{e}^{-\left[ {{\left( \tfrac{T+t}{\eta } \right)}^{\beta }}-{{\left( \tfrac{T}{\eta } \right)}^{\beta }} \right]}}</math>
| |
| | |
| | |
| The above equation gives the reliability for a new mission of <math>t</math> duration, having already accumulated <math>T</math> hours of operation up to the start of this new mission, and the units are checked out to assure that they will start the next mission successfully. (It is called conditional because you can calculate the reliability of a new mission based on the fact that the unit(s) already accumulated <math>T</math> hours of operation successfully.)
| |
| | |
| <br>
| |
| | |
| | |
| ====The Reliable Life====
| |
| For the 2-parameter Weibull distribution, the reliable life, <math>{{T}_{R}}</math> , of a unit for a specified reliability, starting the mission at age zero, is given by:
| |
| | |
| | |
| ::<math>{{T}_{R}}=\eta \cdot {{\left\{ -\ln \left[ R\left( {{T}_{R}} \right) \right] \right\}}^{\tfrac{1}{\beta }}}</math>
| |
| | |
| | |
| This is the life for which the unit will function successfully with a reliability of <math>R({{T}_{R}})</math> . If <math>R({{T}_{R}})=0.50</math> then <math>{{T}_{R}}=\breve{T}</math>,
| |
| the median life, or the life by which half of the units will survive.
| |
| | |
| <br>
| |
| | |
| | |
| ==== The Failure Rate Function ====
| |
| | |
| The 2-parameter Weibull failure rate function,<span class="texhtml"> λ(''T'')</span>, is given by:
| |
| | |
| <br>
| |
| | |
| ::<math>\lambda \left( T \right)=\frac{f\left( T \right)}{R\left( T \right)}=\frac{\beta }{\eta }{{\left( \frac{T}{\eta } \right)}^{\beta -1}}</math>
| |
| | |
| <br>
| |
| ==Characteristics== | |
| The characteristics of the 2-parameter Weibull distribution can be exemplified by examining the two parameters, beta, <math>\beta ,</math> and eta, <math>\eta ,</math> and the effect they have on the <math>pdf,</math> reliability and failure rate functions.
| |
| | |
| ==== Looking at <span class="texhtml">β</span> ==== | |
| | |
| Beta (<span class="texhtml">β</span>) is called the shape parameter or slope of the Weibull distribution. Changing the value of <span class="texhtml">β</span> forces a change in the shape of the <span class="texhtml">''pdf''</span> as shown in the next figure. In addition, when the <span class="texhtml">''cdf''</span> is plotted on Weibull probability paper, a change in beta is a change in the slope of the distribution on Weibull probability paper.
| |
| | |
| <br> '''Effects of <span class="texhtml">β</span> on the ''pdf'''''
| |
| | |
| <br> [[Image:ALTA4.3.png|center|400px]] <br>
| |
| | |
| :*For <span class="texhtml">0 < β < 1</span> , the failure rate decreases with time and:
| |
| ::*As <math>T\to 0,</math> <math>f(T)\to \infty</math>.
| |
| ::*As <math>T\to \infty </math> , <math>f(T)\to 0</math>.
| |
| ::*<span class="texhtml">''f''(''T'')</span> decreases monotonically and is convex as <span class="texhtml">''T''</span> increases.
| |
| ::*The mode is non-existent.
| |
| :*For <span class="texhtml">β = 1,</span> it becomes the exponential distribution, as a special case, <br>
| |
| | |
| <br>
| |
| | |
| :or:
| |
| | |
| <br>
| |
| | |
| ::<math>f(T)=\frac{1}{\eta }{{e}^{-\tfrac{T}{\eta }}};\text{ }\eta >0,T\ge 0</math>
| |
| | |
| <br>
| |
| | |
| :where <math>\tfrac{1}{\eta }=\lambda =</math> chance, useful life, or failure rate. | |
| | |
| <br>
| |
| | |
| :*For <span class="texhtml">β > 1</span>, <span class="texhtml">''f''(''T'')</span>, the Weibull distribution assumes wear-out type shapes (i.e., the failure rate increases with time) and:
| |
| ::*<span class="texhtml">''f''(''T'') = 0</span> at <span class="texhtml">''T'' = 0</span> .
| |
| ::*<span class="texhtml">''f''(''T'')</span> increases as <math>T\to \tilde{T}</math> (mode) and decreases thereafter.
| |
| ::*For <span class="texhtml">β = 2</span> it becomes the Rayleigh distribution as a special case. For <span class="texhtml">β < 2.6</span>, the Weibull <span class="texhtml">''pdf''</span> is positively skewed (has a right tail). For <span class="texhtml">2.6 < β < 3.7</span>, its coefficient of skewness approaches zero (no tail). Consequently, it may approximate the normal <span class="texhtml">''pdf''</span>, and for <span class="texhtml">β > 3.7</span> it is negatively skewed (left tail).
| |
| :*The parameter <span class="texhtml">β</span> is a pure number (i.e., it is dimensionless).
| |
| | |
| <br> '''Effects of <span class="texhtml">β</span> on the Reliability Function and the ''cdf'''''
| |
| | |
| <br> [[Image:ALTA4.4.png|center|400px]]
| |
| | |
| <br> [[Image:ALTA4.5.png|center|400px]]
| |
| | |
| <br>
| |
| | |
| :*<span class="texhtml">''R''(''T'')</span> decreases sharply and monotonically for <span class="texhtml">0 < β < 1</span>. It is convex and decreases less sharply for the same <span class="texhtml">β</span>.
| |
| :*For <span class="texhtml">β = 1</span> and the same <span class="texhtml">η</span>, <span class="texhtml">''R''(''T'')</span> decreases monotonically but less sharply than for <span class="texhtml">0 < β < 1</span>, and it is convex.
| |
| :*For <span class="texhtml">β > 1</span>, <span class="texhtml">''R''(''T'')</span> decreases as <span class="texhtml">''T''</span> increases but less sharply than before. As wear-out sets in, it decreases sharply and goes through an inflection point.
| |
| | |
| <br>
| |
| | |
| '''Effects of <span class="texhtml">β</span> on the Failure Rate Function'''
| |
| | |
| <br> [[Image:ALTA4.6.png|center|400px]] <br>
| |
| | |
| :*The Weibull failure rate for <span class="texhtml">0 < β < 1</span> is unbounded at <span class="texhtml">''T'' = 0</span>. The failure rate, <span class="texhtml">λ(''T'')</span>, decreases thereafter monotonically and is convex, approaching the value of zero as <math>T\to \infty </math> or <math>\lambda (\infty )=0</math>. This behavior makes it suitable for representing the failure rate of units exhibiting early-type failures, for which the failure rate decreases with age. When such behavior is encountered, one or more of the following conclusions can be drawn:
| |
| ::*Burn-in testing and/or environmental stress screening are not well implemented.
| |
| ::*There are problems in the production line.
| |
| ::*There is inadequate quality control.
| |
| ::*There are packaging and transit problems.
| |
| :*For <span class="texhtml">β = 1</span>, <span class="texhtml">λ(''T'')</span> yields a constant value of <math>\tfrac{1}{\eta }</math> , or:
| |
| ::<math>\lambda (T)=\lambda =\frac{1}{\eta }</math>
| |
| | |
| This makes it suitable for representing the failure rate of chance-type failures and the useful life period failure rate of units.
| |
| | |
| :*For <span class="texhtml">β > 1</span>, <span class="texhtml">λ(''T'')</span> increases as <span class="texhtml">''T''</span> increases and becomes suitable for representing the failure rate of units exhibiting wear-out type failures. For <span class="texhtml">1 < β < 2</span>, the <span class="texhtml">λ(''T'')</span> curve is concave. Consequently, the failure rate increases at a decreasing rate as <span class="texhtml">''T''</span> increases.
| |
| :*For <span class="texhtml">β = 2</span>, or for the Rayleigh distribution case, the failure rate function is given by:
| |
| ::<math>\lambda (T)=\frac{2}{\eta }\left( \frac{T}{\eta } \right)</math>
| |
| | |
| <br> Hence there emerges a straight line relationship between <span class="texhtml">λ(''T'')</span> and <span class="texhtml">''T''</span>, starting at a value of <span class="texhtml">λ(''T'') = 0</span> at <span class="texhtml">''T'' = 0</span> and increasing thereafter with a slope of <math>\tfrac{2}{{{\eta }^{2}}}</math> . Consequently, the failure rate increases at a constant rate as <span class="texhtml">''T''</span> increases. Furthermore, if <span class="texhtml">η = 1</span> the slope becomes equal to 2, and <span class="texhtml">λ(''T'')</span> becomes a straight line which passes through the origin with a slope of 2.
| |
| | |
| :*When <span class="texhtml">β > 2</span> the <span class="texhtml">λ(''T'')</span> curve is convex, with its slope increasing as <span class="texhtml">''T''</span> increases. Consequently, the failure rate increases at an increasing rate as <span class="texhtml">''T''</span> increases, indicating wear-out life.
| |
| | |
| <br>
| |
| | |
| ==== Looking at <span class="texhtml">η</span> ====
| |
| | |
| Eta (<span class="texhtml">η</span>) is called the scale parameter of the Weibull distribution. The parameter <span class="texhtml">η</span> has the same units as <span class="texhtml">''T''</span>, such as hours, miles, cycles, actuations, etc.
| |
| | |
| <br> [[Image:ALTA4.7.png|center|400px]] <br>
| |
| | |
| :*A change in the scale parameter <span class="texhtml">η</span> has the same effect on the distribution as a change of the abscissa scale.
| |
| ::*If <span class="texhtml">η</span> is increased while <span class="texhtml">β</span> is kept the same, the distribution gets stretched out to the right and its height decreases, while maintaining its shape and location.
| |
| ::*If <span class="texhtml">η</span> is decreased while <span class="texhtml">β</span> is kept the same, the distribution gets pushed in toward the left (i.e. toward its beginning, or 0) and its height increases.
| |
| | |
| <br>
| |
| | |
| == Parameter Estimation ==
| |
| | |
| The estimates of the parameters of the Weibull distribution can be found graphically on probability plotting paper, or analytically using either least squares or maximum likelihood. (Parameter estimation methods are presented in detail in [[Appendix B: Parameter Estimation|Appendix B]].) <br>
| |
| | |
| ====Probability Plotting====
| |
| | |
| One method of calculating the parameters of the Weibull distribution is by using probability plotting. To better illustrate this procedure, consider the following example [[Reference Appendix D: References|[18]]].
| |
| <br>
| |
| <br>
| |
| '''Example 3'''
| |
| <br>
| |
| | |
| Let's assume six identical units are being reliability tested at the same application and operation stress levels. All of these units fail during the test after operating the following times (in hours), <span class="texhtml">''T''<sub>''i''</sub></span> : 93, 34, 16, 120, 53 and 75. The steps for using probability plotting to determine the parameters of the Weibull <span class="texhtml">''pdf''</span> are as follows: <br>
| |
| | |
| :*Rank the failure times in ascending order as shown next.
| |
| | |
| [[Image:failureordernumber.png|center|300px]]
| |
| | |
| :*Obtain their median rank plotting positions. The failure, with their corresponding median ranks, are shown next.
| |
| | |
| [[Image:median rank.png|center|300px]]
| |
| | |
| :*On a Weibull probability paper, plot the times and their corresponding ranks. The next figure displays an example of a Weibull probability paper.
| |
| | |
| <br> [[Image:ALTA4.8.png|center|400px]] <br>
| |
| | |
| :*Draw the best possible straight line through the plotted points.
| |
| :*Obtain the slope of this line by drawing a line, parallel to the one just obtained, through the slope indicator. This value is the estimate of the shape parameter <math>\widehat{\beta }</math>. In this case <math>\widehat{\beta }=1.4</math>.
| |
| :*At the <span class="texhtml">''Q''(''t'') = 63.2%</span> ordinate point, draw a straight horizontal line until this line intersects the fitted straight line. Draw a vertical line through this intersection until it crosses the abscissa. The value at the intersection of the abscissa is the estimate of <math>\widehat{\eta }</math> . For this case <math>\widehat{\eta }=76</math> hours. (This is always at 63.2% since <math>Q(T)=1-{{e}^{-{{(\tfrac{\eta }{\eta })}^{\beta }}}}=1-{{e}^{-1}}=0.632=63.2%</math>.)
| |
| | |
| <br> [[Image:ALTA4.9.png|center|400px]] <br> Now any reliability value for any mission time <span class="texhtml">''t''</span> can be obtained. For example, the reliability for a mission of 15 hours (or any other time) can now be obtained either from the plot or analytically.
| |
| | |
| To obtain the value from the plot, draw a vertical line from the abscissa, at <span class="texhtml">''t'' = 15</span> hours, to the fitted line. Draw a horizontal line from this intersection to the ordinate and read <span class="texhtml">''Q''(''t'')</span>, in this case <span class="texhtml">''Q''(''t'' = 15) = 9.8%</span>. Thus, <span class="texhtml">''R''(''t'' = 15) = 1 − ''Q''(''t'') = 90.2%</span>. This can also be obtained analytically from the Weibull reliability function since both of the parameters are known. <br>
| |
| | |
| ::<math>R(t=15)={{e}^{-{{\left( \tfrac{15}{\eta } \right)}^{\beta }}}}={{e}^{-{{\left( \tfrac{15}{76} \right)}^{1.4}}}}=90.2%</math>
| |
| | |
| <br>
| |
| | |
| ====MLE Parameter Estimation====
| |
| The parameters of the 2-parameter Weibull distribution can also be estimated using maximum likelihood estimation (MLE). This log-likelihood function is composed of :
| |
| | |
| <br>
| |
| ::<math>\begin{align}
| |
| & \ln (L)= & \Lambda =\underset{i=1}{\overset{{{F}_{e}}}{\mathop \sum }}\,{{N}_{i}}\ln \left[ \frac{\beta }{\eta }{{\left( \frac{{{T}_{i}}}{\eta } \right)}^{\beta -1}}{{e}^{-{{\left( \tfrac{{{T}_{i}}}{\eta } \right)}^{\beta }}}} \right]-\underset{i=1}{\overset{S}{\mathop \sum }}\,N_{i}^{\prime }{{\left( \frac{T_{i}^{\prime }}{\eta } \right)}^{\beta }}\overset{FI}{\mathop{+\underset{i=1}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{\prime \prime }\ln [R_{Li}^{\prime \prime }-R_{Ri}^{\prime \prime }]
| |
| \end{align}</math>
| |
| <br>
| |
| :where:
| |
| <br>
| |
| ::<math>R_{Li}^{\prime \prime }={{e}^{-{{(\tfrac{T_{Li}^{\prime \prime }}{\eta })}^{\beta }}}}</math>
| |
| <br>
| |
| ::<math>R_{Ri}^{\prime \prime }={{e}^{-{{(\tfrac{T_{Ri}^{\prime \prime }}{\eta })}^{\beta }}}}</math>
| |
| <br>
| |
| :* <math>{{F}_{e}}</math> is the number of groups of times-to-failure data points.
| |
| :* <math>{{N}_{i}}</math> is the number of times-to-failure in the <math>{{i}^{th}}</math> time-to-failure data group.
| |
| :* <math>\beta </math> is the Weibull shape parameter (unknown a priori, the first of two parameters to be found).
| |
| :* <math>\eta </math> is the Weibull scale parameter (unknown a priori, the second of two parameters to be found).
| |
| :* <math>{{T}_{i}}</math> is the time of the <math>{{i}^{th}}</math> group of time-to-failure data.
| |
| :* <math>S</math> is the number of groups of suspension data points.
| |
| :* <math>N_{i}^{\prime }</math> is the number of suspensions in <math>{{i}^{th}}</math> group of suspension data points.
| |
| :* <math>T_{i}^{\prime }</math> is the time of the <math>{{i}^{th}}</math> suspension data group.
| |
| :* <math>FI</math> is the number of interval data groups.
| |
| :* <math>N_{i}^{\prime \prime }</math> is the number of intervals in the <math>{{i}^{th}}</math> group of data intervals.
| |
| :* <math>T_{Li}^{\prime \prime }</math> is the beginning of the <math>{{i}^{th}}</math> interval.
| |
| :* <math>T_{Ri}^{\prime \prime }</math> is the ending of the <math>{{i}^{th}}</math> interval.
| |
| <br>
| |
| The solution is found by solving for a pair of parameters <math>\left( \widehat{\beta },\widehat{\eta } \right)</math> so that <math>\tfrac{\partial \Lambda }{\partial \beta }=0</math> and <math>\tfrac{\partial \Lambda }{\partial \eta }=0</math>. (Other methods can also be used, such as direct maximization of the likelihood function, without having to compute the derivatives.)
| |
| <br>
| |
| ::<math>\begin{align}
| |
| &\frac{\partial \Lambda }{\partial \beta }= \frac{1}{\beta }\underset{i=1}{\overset{{{F}_{e}}}{\mathop{\sum }}}\,{{N}_{i}}+\underset{i=1}{\overset{{{F}_{e}}}{\mathop{\sum }}}\,{{N}_{i}}\ln \left( \frac{{{T}_{i}}}{\eta } \right) -\underset{i=1}{\overset{{{F}_{e}}}{\mathop{\sum }}}\,{{N}_{i}}{{\left( \frac{{{T}_{i}}}{\eta } \right)}^{\beta }}\ln \left( \frac{{{T}_{i}}}{\eta } \right)-\underset{i=1}{\overset{S}{\mathop{\sum }}}\,N_{i}^{\prime }{{\left( \frac{T_{i}^{\prime }}{\eta } \right)}^{\beta }}\ln \left( \frac{T_{i}^{\prime }}{\eta } \right)+\overset{FI}{\mathop{\underset{i=1}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{\prime \prime }\frac{-{{(\tfrac{T_{Li}^{\prime \prime }}{\eta })}^{\beta }}\ln (\tfrac{T_{Li}^{\prime \prime }}{\eta })R_{Li}^{\prime \prime }+{{(\tfrac{T_{Ri}^{\prime \prime }}{\eta })}^{\beta }}\ln (\tfrac{T_{Ri}^{\prime \prime }}{\eta })R_{Ri}^{\prime \prime }}{R_{Li}^{\prime \prime }-R_{Ri}^{\prime \prime }} \
| |
| &\frac{\partial \Lambda }{\partial \eta }= \frac{-\beta }{\eta }\underset{i=1}{\overset{{{F}_{e}}}{\mathop{\sum }}}\,{{N}_{i}}+\frac{\beta }{\eta }\underset{i=1}{\overset{{{F}_{e}}}{\mathop{\sum }}}\,{{N}_{i}}{{\left( \frac{{{T}_{i}}}{\eta } \right)}^{\beta }}+\frac{\beta }{\eta }\underset{i=1}{\overset{S}{\mathop{\sum }}}\,N_{i}^{\prime }{{\left( \frac{T_{i}^{\prime }}{\eta } \right)}^{\beta }}+\overset{FI}{\mathop{\underset{i=1}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{\prime \prime }\frac{\beta }{\eta }\frac{{{(\tfrac{T_{Li}^{\prime \prime }}{\eta })}^{\beta }}R_{Li}^{\prime \prime }-{{(\tfrac{T_{Ri}^{\prime \prime }}{\eta })}^{\beta }}R_{Ri}^{\prime \prime }}{R_{Li}^{\prime \prime }-R_{Ri}^{\prime \prime }}
| |
| \end{align}</math>.
| |
| <br>
| |
| <br>
| |
| '''Example 4'''
| |
| Using the same data as in the probability plotting example (Example 3), and assuming a 2-parameter Weibull distribution, estimate the parameter using the MLE method.
| |
| <br>
| |
| | |
| '''Solution'''
| |
| <br>
| |
| In this case we have non-grouped data with no suspensions, therefore the above equations become:
| |
| | |
| <br>
| |
| ::<math>\frac{\partial \Lambda }{\partial \beta }=\frac{6}{\beta }+\underset{i=1}{\overset{6}{\mathop{\sum }}}\,\ln \left( \frac{{{T}_{i}}}{\eta } \right)-\underset{i=1}{\overset{6}{\mathop{\sum }}}\,{{\left( \frac{{{T}_{i}}}{\eta } \right)}^{\beta }}\ln \left( \frac{{{T}_{i}}}{\eta } \right)=0</math>
| |
| <br>
| |
| :and:
| |
| <br>
| |
| ::<math>\frac{\partial \Lambda }{\partial \eta }=\frac{-\beta }{\eta }\cdot 6+\frac{\beta }{\eta }\underset{i=1}{\overset{6}{\mathop \sum }}\,{{\left( \frac{{{T}_{i}}}{\eta } \right)}^{\beta }}=0</math>
| |
| <br>
| |
| Solving the above equations simultaneously we get:
| |
| | |
| <br>
| |
| ::<math>\begin{matrix}
| |
| \widehat{\beta }=1.933 \
| |
| \widehat{\eta }=73.526 \
| |
| \end{matrix}</math>
| |
| <br>
| |
|
| |
|
| = The Lognormal Distribution = | | = The Lognormal Distribution = |
| | | {{Lognormal Distribution Definition}} |
| The lognormal distribution is commonly used for general reliability analysis, cycles-to-failure in fatigue, material strengths and loading variables in probabilistic design. A random variable is lognormally distributed if the logarithm of the random variable is normally distributed. Since the logarithms of a lognormally distributed random variable are normally distributed, the lognormal distribution is given by: <br>
| | ==Lognormal Distribution Functions== |
| | | {{:Lognormal Distribution Functions}} |
| <br>
| | ==Characteristics of the Lognormal Distribution== |
| | | {{:Lognormal Distribution Characteristics}} |
| ::<math>f({T}')=\frac{1}{{{\sigma }_{{{T}'}}}\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{\left( \tfrac{{T}'-\bar{{T}'}}{{{\sigma }_{{{T}'}}}} \right)}^{2}}}}</math>
| |
| | |
| <br>
| |
| | |
| where:
| |
| | |
| :* <span class="texhtml">''T''' = ln ''T''</span>, and where the <span class="texhtml">''T''</span>s are the failure times.
| |
| :* <math>\bar{{T}'}=</math> mean of the natural logarithms of the times to failure.
| |
| :* <span class="texhtml">σ<sub>''T'''</sub> = </span> standard deviation of the natural logarithms of the failure times.
| |
| | |
| <br> The lognormal <span class="texhtml">''pdf''</span> can be obtained, realizing that for equal probabilities under the normal and lognormal <span class="texhtml">''pdf''</span> s incremental areas should also be equal, or:
| |
| | |
| <br>
| |
| | |
| ::<math>\begin{align}
| |
| f(T)dT = f(T')dT'
| |
| \end{align}</math>
| |
| | |
| <br> Taking the derivative yields:
| |
| | |
| <br>
| |
| | |
| ::<math>d{T}'=\frac{dT}{T}</math>
| |
| | |
| <br> Substitution yields:
| |
| | |
| <br>
| |
| | |
| ::<math>\begin{align}
| |
| f(T)= \frac{f({T}')}{T}= \frac{1}{T\cdot {{\sigma }_{{{T}'}}}\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{\left( \tfrac{{T}'-\bar{{T}'}}{{{\sigma }_{{{T}'}}}} \right)}^{2}}}}
| |
| \end{align}</math>
| |
| | |
| <br>
| |
| | |
| where: <br>
| |
| | |
| ::<math>f(T)\ge 0,T>0,-\infty <\bar{{T}'}<\infty ,{{\sigma }_{{{T}'}}}>0</math>
| |
| | |
| <br>
| |
| | |
| ==Statistical Properties Summary==
| |
| ====The Mean or MTTF====
| |
| :* The mean of the lognormal distribution, <math>\bar{T}</math>, is given by:
| |
| | |
| | |
| ::<math>\bar{T}={{e}^{\bar{{T}'}+\tfrac{1}{2}\sigma _{{{T}'}}^{2}}}</math>
| |
| | |
| | |
| :* The mean of the natural logarithms of the failure times, <math>{{\bar{T}}^{^{\prime }}}</math>, in terms of <math>\bar{T}</math> and <math>{{\sigma }_{T}}</math> is given by:
| |
| <br>
| |
| ::<math>{{\bar{T}}^{\prime }}=\ln \left( {\bar{T}} \right)-\frac{1}{2}\ln \left( \frac{\sigma _{T}^{2}}{{{{\bar{T}}}^{2}}}+1 \right)</math>
| |
| | |
| <br>
| |
| | |
| | |
| ====The Standard Deviation====
| |
| :* The standard deviation of the lognormal distribution, <math>{{\sigma }_{T}}</math>, is given by:
| |
| <br>
| |
| | |
| ::<math>{{\sigma }_{T}}=\sqrt{\left( {{e}^{2\bar{{T}'}+\sigma _{{{T}'}}^{2}}} \right)\left( {{e}^{\sigma _{{{T}'}}^{2}}}-1 \right)}</math>
| |
| <br>
| |
| | |
| :* The standard deviation of the natural logarithms of the times-to-failure, <math>{{\sigma }_{{{T}'}}}</math>, in terms of <math>\bar{T}</math> and <math>{{\sigma }_{T}}</math> is given by:
| |
| <br>
| |
| | |
| ::<math>{{\sigma }_{{{T}'}}}=\sqrt{\ln \left( \frac{\sigma _{T}^{2}}{{{{\bar{T}}}^{2}}}+1 \right)}</math>
| |
| <br>
| |
| | |
| | |
| ====The Median====
| |
| :* The median of the lognormal distribution is given by:
| |
| <br>
| |
| | |
| ::<math>\breve{T}={{e}^{{{\bar{T}}^{\prime }}}}</math>
| |
| <br>
| |
| | |
| | |
| ====The Mode====
| |
| :* The mode of the lognormal distribution is given by:
| |
| <br>
| |
| | |
| ::<math>\tilde{T}={{e}^{{{\bar{T}}^{\prime }}-\sigma _{{{T}'}}^{2}}}</math>
| |
| <br>
| |
| | |
| | |
| ====Reliability Function====
| |
| For the lognormal distribution, the reliability for a mission of time <math>T</math>, starting at age 0, is given by:
| |
| | |
| <br>
| |
| ::<math>R(T)=\int_{T}^{\infty }f(t)dt</math>
| |
| | |
| or:
| |
| | |
| ::<math>R(T)=\int_{{{T}^{^{\prime }}}}^{\infty }\frac{1}{{{\sigma }_{{{T}'}}}\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{\left( \tfrac{t-\overline{{{T}'}}}{{{\sigma }_{{{T}'}}}} \right)}^{2}}}}dt</math>
| |
| <br>
| |
| | |
| There is no closed form solution for the lognormal reliability function. Solutions can be obtained via the use of standard normal tables.
| |
| | |
| <br>
| |
| | |
| | |
| ====Lognormal Failure Rate====
| |
| The lognormal failure rate is given by:
| |
| <br>
| |
| ::<math>\lambda (T)=\frac{f(T)}{R(T)}=\frac{\tfrac{1}{{T}'{{\sigma }_{{{T}'}}}\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{(\tfrac{{T}'-\overline{{{T}'}}}{{{\sigma }_{{{T}'}}}})}^{2}}}}}{\mathop{}_{{{T}'}}^{\infty }\tfrac{1}{{{\sigma }_{{{T}'}}}\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{(\tfrac{t-\overline{{{T}'}}}{{{\sigma }_{{{T}'}}}})}^{2}}}}dt}</math>
| |
| | |
| ==Characteristics==
| |
| :* The lognormal distribution is a distribution skewed to the right.
| |
| :* The <math>pdf</math> starts at zero, increases to its mode, and decreases thereafter.
| |
| <br>
| |
| [[Image:chp4pdf.png|center|400px|''Pdf'' of the lognormal distribution.]]
| |
| <br>
| |
| <br>
| |
| The characteristics of the lognormal distribution can be exemplified by examining the two parameters, the log-mean <math>({{\overline{T}}^{\prime }})</math> and the log-std (<math>{{\sigma }_{{{T}'}}}</math>), and the effect they have on the <math>pdf</math>.
| |
| <br>
| |
| '''Looking at the Log-Mean''' <math>({{\overline{T}}^{\prime }})</math>
| |
| :* The parameter, <math>\bar{{T}'}</math>, or the log-mean life, or the <math>MTT{F}'</math> in terms of the logarithm of the <math>{T}'s</math> is also the scale parameter and a unitless number.
| |
| :* For the same <math>{{\sigma }_{{{T}'}}}</math> the <math>pdf</math> 's skewness increases as <math>\bar{{T}'}</math> increases.
| |
| <br>
| |
| <br>
| |
| [[Image:chp4pdf2.png|center|400px|''Pdf'' of the lognormal distribution with different log-mean values.]]
| |
| <br>
| |
| | |
|
| |
| ====Looking at the Log-STD <math>({{\sigma }_{{{T}'}}})</math>====
| |
| :* The parameter <math>{{\sigma }_{{{T}'}}}</math>, or the standard deviation of the <math>{T}'s</math> in terms of their logarithm or of their <math>{T}'</math>, is also the shape parameter, and not the scale parameter as in the normal <math>pdf</math>. It is a unitless number and assumes only positive values.
| |
| :* The degree of skewness increases as <math>{{\sigma }_{{{T}'}}}</math> increases, for a given <math>\bar{{T}'}</math>.
| |
| :* For <math>{{\sigma }_{{{T}'}}}</math> values significantly greater than 1, the <math>pdf</math> rises very sharply in the beginning (i.e., for very small values of <math>T</math> near zero), and essentially follows the ordinate axis, peaks out early, and then decreases sharply like an exponential <math>pdf</math> or a Weibull <math>pdf</math> with <math>0<\beta <1</math>.
| |
| <br>
| |
| [[Image:chp4pdf3.png|center|400px|''Pdf'' of the lognormal distribution with different log-std values.]]
| |
| <br>
| |
| | |
| == Parameter Estimation ==
| |
| The estimate of the parameters of the lognormal distribution can be found graphically on probability plotting paper or analytically using either least squares or maximum likelihood. (Parameter estimation methods are presented in detail in [[Appendix B: Parameter Estimation|Appendix B]].) <br>
| |
| | |
| ====Probability Plotting====
| |
| One method of calculating the parameter of the lognormal distribution is by using probability plotting. To better illustrate this procedure, consider the following example.
| |
| <br>
| |
| <br>
| |
| '''Example 5'''
| |
| | |
| Let's assume six identical units are being reliability tested at the same application and operation stress levels. All of these units fail during the test after operating the following times (in hours), <span class="texhtml">''T''<sub>''i''</sub></span> : 144, 385, 747, 1,144, 1,576 and 2,616. The steps for determining the parameters of the lognormal <span class="texhtml">''pdf''</span> representing the data, using probability plotting, are as follows: <br>
| |
| | |
| :* Rank the failure times in ascending order as shown next.
| |
| <center><math>\begin{matrix}
| |
| \text{Failure} & \text{Failure Order Number} \
| |
| \text{Time (Hr)} & \text{out of a Sample Size of 6} \
| |
| \text{144} & \text{1} \
| |
| \text{385} & \text{2} \
| |
| \text{747} & \text{3} \
| |
| \text{1,144} & \text{4} \
| |
| \text{1,576} & \text{5} \
| |
| \text{2,616} & \text{6} \
| |
| \end{matrix}</math></center>
| |
| <br>
| |
| | |
| :* Obtain their median rank plotting positions. The failure times, with their corresponding median ranks, are shown next:
| |
| | |
| <br>
| |
| <center><math>\begin{matrix}
| |
| \text{Failure} & \text{Median} \
| |
| \text{Time (Hr)} & \text{Rank, }% \
| |
| \text{144} & \text{10}\text{.91} \
| |
| \text{385} & \text{26}\text{.44} \
| |
| \text{747} & \text{42}\text{.14} \
| |
| \text{1,144} & \text{57}\text{.86} \
| |
| \text{1,576} & \text{73}\text{.56} \
| |
| \text{2,616} & \text{89}\text{.09} \
| |
| \end{matrix}</math></center>
| |
| <br>
| |
| | |
| :* On a lognormal probability paper, plot the times and their corresponding rank value. The next figure displays an example of a lognormal probability paper. The paper is simply a log-log paper.
| |
| | |
| <br> [[Image:ALTA4.10.png|center|400px]] <br>
| |
| | |
| :* Draw the best possible straight line that goes through the <span class="texhtml">''t'' = 0</span> and <span class="texhtml">''R''(''t'') = 100%</span> point and through these points.
| |
| :* At the <span class="texhtml">''Q''(''t'') = 50%</span> ordinate point, draw a straight horizontal line until this line intersects the fitted straight line. Draw a vertical line through this intersection until it crosses the abscissa. The value at the intersection of the abscissa is the estimate of the median. For this case, <math>\breve{T}=760</math> hours which means that <math>{{\bar{T}}^{\prime }}=\ln(\breve{T})=6.633</math>.
| |
| | |
| <br>
| |
| | |
| ::
| |
| | |
| <br> [[Image:ALTA4.11.png|center|400px]] <br> <br>
| |
| | |
| :* The standard deviation, <span class="texhtml">σ<sub>''T'''</sub>,</span> can be found using the following equation:
| |
| | |
| <br>
| |
| | |
| ::<math>\begin{align}
| |
| {{\sigma }_{{{T}'}}}= & \frac{\ln \left[ T(Q=97.7%) \right]-\ln \left[ T(Q=2.3%) \right]}{4} \
| |
| = & \frac{\ln (5100)-\ln (120)}{4} \
| |
| = & 0.937376
| |
| \end{align}</math>
| |
| | |
| <br> Now any reliability value for any mission time <span class="texhtml">''t''</span> can be obtained. For example, the reliability for a mission of 200 hours, or any other time, can now be obtained either from the plot or analytically. <br> To obtain the value from the plot, draw a vertical line from the abscissa, at <span class="texhtml">''t'' = 200</span> hours, to the fitted line. Draw a horizontal line from this intersection to the ordinate and read <span class="texhtml">''Q''(''t'')</span> . In this case, <span class="texhtml">''R''(''t'' = 200) = 1 − ''Q''(''t'' = 200) = 92%</span> . This can also be obtained analytically, from the lognormal reliability function. However, standard normal tables (or the Quick Statistical Reference in ALTA) must be used.
| |
| | |
| ====MLE Parameter Estimation====
| |
| The parameters of the lognormal distribution can also be estimated using maximum likelihood estimation (MLE). This general log-likelihood function is:
| |
| | |
| <br>
| |
| ::<math>\begin{align}
| |
| & \ln (L)= \Lambda =\underset{i=1}{\overset{{{F}_{e}}}{\mathop \sum }}\,{{N}_{i}}\ln \left[ \frac{1}{{{\sigma }_{{{T}'}}}{{T}_{i}}}\phi \left( \frac{\ln \left( {{T}_{i}} \right)-{\mu }'}{{{\sigma }_{{{T}'}}}} \right) \right]\text{ }+\underset{i=1}{\overset{S}{\mathop \sum }}\,N_{i}^{\prime }\ln \left[ 1-\Phi \left( \frac{\ln \left( T_{i}^{\prime } \right)-{\mu }'}{{{\sigma }_{{{T}'}}}} \right) \right]+\overset{FI}{\mathop{\underset{i=1}{\mathop{\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{\prime \prime }\ln [\Phi (z_{Ri}^{\prime \prime })-\Phi (z_{Li}^{\prime \prime })]
| |
| \end{align}</math>
| |
| | |
| <br>
| |
| where:
| |
| | |
| <br>
| |
| ::<math>z_{Li}^{\prime \prime }=\frac{\ln T_{Li}^{\prime \prime }-{\mu }'}{\sigma _{T}^{\prime }}</math>
| |
| <br>
| |
| ::<math>z_{Ri}^{\prime \prime }=\frac{\ln T_{Ri}^{\prime \prime }-{\mu }'}{\sigma _{T}^{\prime }}</math>
| |
| <br>
| |
| and:
| |
| <br>
| |
| <br>
| |
| :* <math>{{F}_{e}}</math> is the number of groups of times-to-failure data points.
| |
| :* <math>{{N}_{i}}</math> is the number of failure times in the <math>{{i}^{th}}</math> time-to-failure data group.
| |
| :* <math>{\mu }'</math> is the mean of the natural logarithms of the failure times (unknown a priori, the first of two parameters to be found).
| |
| :* <math>{{\sigma }_{{{T}'}}}</math> is the standard deviation of the natural logarithms of the failure times (unknown a priori, the second of two parameters to be found).
| |
| :* <math>{{T}_{i}}</math> is the time of the <math>{{i}^{th}}</math> group of time-to-failure data.
| |
| :* <math>S</math> is the number of groups of suspension data points.
| |
| :* <math>N_{i}^{\prime }</math> is the number of suspensions in <math>{{i}^{th}}</math> group of suspension data points.
| |
| :* <math>T_{i}^{\prime }</math> is the time of the <math>{{i}^{th}}</math> suspension data group.
| |
| :* <math>FI</math> is the number of interval data groups.
| |
| :* <math>N_{i}^{\prime \prime }</math> is the number of intervals in the <math>i^{th}</math> group of data intervals.
| |
| :* <math>T_{Li}^{\prime \prime }</math> is the beginning of the <math>i^{th}</math> interval.
| |
| :* <math>T_{Ri}^{\prime \prime }</math> is the ending of the <math>i^{th}</math> interval.
| |
| | |
| <br>
| |
| The solution will be found by solving for a pair of parameters <math>\left( {\mu }',{{\sigma }_{{{T}'}}} \right)</math> so that <math>\tfrac{\partial \Lambda }{\partial {\mu }'}=0</math> and <math>\tfrac{\partial \Lambda }{\partial {{\sigma }_{{{T}'}}}}=0</math>, where:
| |
| <br>
| |
| ::<math>\begin{align}
| |
| & \frac{\partial \Lambda }{\partial {\mu }'}= \frac{1}{\sigma _{{{T}'}}^{2}}\underset{i=1}{\overset{{{F}_{e}}}{\mathop \sum }}\,{{N}_{i}}(\ln ({{T}_{i}})-{\mu }')+\frac{1}{{{\sigma }_{{{T}'}}}}\underset{i=1}{\overset{S}{\mathop \sum }}\,N_{i}^{\prime }\frac{\phi \left( \tfrac{\ln \left( T_{i}^{\prime } \right)-{\mu }'}{{{\sigma }_{{{T}'}}}} \right)}{1-\Phi \left( \tfrac{\ln \left( T_{i}^{\prime } \right)-{\mu }'}{{{\sigma }_{{{T}'}}}} \right)}\overset{FI}{\mathop{\underset{i=1}{\mathop{-\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{\prime \prime }\frac{\varphi (z_{Ri}^{\prime \prime })-\varphi (z_{Li}^{\prime \prime })}{\sigma _{T}^{\prime }(\Phi (z_{Ri}^{\prime \prime })-\Phi (z_{Li}^{\prime \prime }))} \
| |
| & & \
| |
| & \frac{\partial \Lambda }{\partial {{\sigma }_{{{T}'}}}}= \underset{i=1}{\overset{{{F}_{e}}}{\mathop \sum }}\,{{N}_{i}}\left( \frac{{{\left( \ln ({{T}_{i}})-{\mu }' \right)}^{2}}}{\sigma _{{{T}'}}^{3}}-\frac{1}{{{\sigma }_{{{T}'}}}} \right)+\frac{1}{{{\sigma }_{{{T}'}}}}\underset{i=1}{\overset{S}{\mathop \sum }}\,N_{i}^{\prime }\frac{\left( \tfrac{\ln \left( T_{i}^{\prime } \right)-{\mu }'}{{{\sigma }_{{{T}'}}}} \right)\phi \left( \tfrac{\ln \left( T_{i}^{\prime } \right)-{\mu }'}{{{\sigma }_{{{T}'}}}} \right)}{1-\Phi \left( \tfrac{\ln \left( T_{i}^{\prime } \right)-{\mu }'}{{{\sigma }_{{{T}'}}}} \right)}\overset{FI}{\mathop{\underset{i=1}{\mathop{-\underset{}{\overset{}{\mathop \sum }}\,}}\,}}\,N_{i}^{\prime \prime }\frac{z_{Ri}^{\prime \prime }\varphi (z_{Ri}^{\prime \prime })-z_{Li}^{\prime \prime }\varphi (z_{Li}^{\prime \prime })}{\sigma _{T}^{\prime }(\Phi (z_{Ri}^{\prime \prime })-\Phi (z_{Li}^{\prime \prime }))}
| |
| \end{align}</math>
| |
| <br>
| |
| and:
| |
| <br>
| |
| <br>
| |
| ::<math>\phi \left( x \right)=\frac{1}{\sqrt{2\pi }}\cdot {{e}^{-\tfrac{1}{2}{{\left( x \right)}^{2}}}}</math>
| |
| <br>
| |
| ::<math>\Phi (x)=\frac{1}{\sqrt{2\pi }}\int_{-\infty }^{x}{{e}^{-\tfrac{{{t}^{2}}}{2}}}dx</math>
| |
| <br>
| |
| <br>
| |
| <br>
| |
| '''Example 6'''
| |
| <br>
| |
| | |
| Using the same data as in the probability plotting example (Example 5) and assuming a lognormal distribution, estimate the parameters using the MLE method.
| |
| <br>
| |
| '''Solution '''
| |
| <br>
| |
| In this example we have non-grouped data without suspensions. Thus, the partials reduce to:<br>
| |
| | |
| <br>
| |
| | |
| ::<math>\begin{align}
| |
| & \frac{\partial \Lambda }{\partial {\mu }'}= & \frac{1}{\sigma _{{{T}'}}^{2}}\cdot \underset{i=1}{\overset{14}{\mathop \sum }}\,\ln ({{T}_{i}})-{\mu }'=0 \
| |
| & \frac{\partial \Lambda }{\partial {{\sigma }_{{{T}'}}}}= & \underset{i=1}{\overset{14}{\mathop \sum }}\,\left( \frac{\ln ({{T}_{i}})-{\mu }'}{\sigma _{{{T}'}}^{3}}-\frac{1}{{{\sigma }_{{{T}'}}}} \right)=0
| |
| \end{align}</math>
| |
| | |
| <br> Substituting the values of <span class="texhtml">''T''<sub>''i''</sub></span> and solving the above system simultaneously, we get: <br>
| |
| | |
| <br>
| |
| | |
| ::<math>\begin{align}
| |
| & {{\sigma }_{{{T}'}}}= & 0.9537 \
| |
| & {\mu }'= & 6.6356
| |
| \end{align}</math>
| |
| | |
| <br> The mean and standard deviation of the times-to-failure can be estimated by: <br>
| |
| | |
| <br>
| |
| | |
| ::<math>\overline{T}=\mu =1,200.31\text{ }hr</math>
| |
| | |
| <br> and: <br>
| |
| | |
| <br>
| |
| | |
| ::<math>\begin{align}
| |
| & {{\sigma }_{T}}= 1,461.78\text{ }hr
| |
| \end{align}</math>
| |
| | |
| <br>
| |
| | |
| [[Category:Accelerated_Life_Data_Analysis_Reference]]
| |