Template:MLE Example (EXP1-2): Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
(Redirected page to The Exponential Distribution)
 
(2 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 
#REDIRECT [[The_Exponential_Distribution]]
====MLE Parameter Estimation====
 
<br>
 
Let's assume six identical units are reliability tested at the same application and operation
stress levels. All of these units fail during the test after operating for the following times (in hours),  <math>{{T}_{i}}</math> : 96, 257, 498, 763, 1051 and 1744.
 
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>
<br>
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>
and
<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>\lambda </math>  is the failure rate parameter (unknown a priori, the only parameter 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 the  <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 i <math>^{th}</math>  group of data intervals.
::• <math>T_{Li}^{\prime \prime }</math>  is the beginning of the i <math>^{th}</math>  interval.
::• <math>T_{Ri}^{\prime \prime }</math>  is the ending of the i <math>^{th}</math>  interval.
 
<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>
 
====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>
''Solution''
<br>
In this example we have non-grouped data without suspensions. Thus Eqn. (exp-mle) becomes:
<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  <math>T</math>  we get:
 
<br>
::<math>\begin{align}
  \frac{6}{\lambda }= & 4409,\text{ or:} \\
  \lambda = & 0.00136\text{ failure/hr} 
\end{align}</math>
 
[[Category:Examples]]

Latest revision as of 06:52, 15 August 2012