Template:Probability Plotting Exponential Distribution: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
(Created page with '===Probability Plotting=== Estimation of the parameters for the exponential distribution via probability plotting is very similar to the process used when dealing with the Weibul…')
 
 
(11 intermediate revisions by 2 users not shown)
Line 1: Line 1:
===Probability Plotting===
#REDIRECT [[The_Exponential_Distribution#Probability_Plotting]]
Estimation of the parameters for the exponential distribution via probability plotting is very similar to the process used when dealing with the Weibull distribution. Recall, however, that the appearance of the probability plotting paper and the methods by which the parameters are estimated vary from distribution to distribution, so there will be some noticeable differences. In fact, due to the nature of the exponential <math>cdf</math>, the exponential probability plot is the only one with a negative slope. This is because the y-axis of the exponential probability plotting paper represents the reliability, whereas the y-axis for most of the other life distributions represents the unreliability.
 
This is illustrated in the process of linearizing the <math>cdf</math>, which is necessary to construct the exponential probability plotting paper. For the two-parameter exponential distribution the cumulative density function is given by:
 
 
::<math>F(T)=1-{{e}^{-\lambda (T-\gamma )}}</math>
 
 
Taking the natural logarithm of both sides of Eqn. (Fe) yields:
 
 
::<math>\ln \left[ 1-F(T) \right]=-\lambda (T-\gamma )</math>
 
 
:or:
 
 
::<math>\ln [1-F(T)]=\lambda \gamma -\lambda T</math>
 
 
:Now, let:
 
 
::<math>y=\ln [1-F(T)]</math>
 
 
::<math>a=\lambda \gamma </math>
 
 
:and:
 
 
::<math>b=-\lambda </math>
 
 
which results in the linear equation of:
 
 
::<math>y=a+bT</math>
 
 
Note that with the exponential probability plotting paper, the y-axis scale is logarithmic and the x-axis scale is linear. This means that the zero value is present only on the x-axis. For <math>t=0</math>, <math>R=1</math> and <math>F(t)=0</math>. So if we were to use <math>F(t)</math> for the y-axis, we would have to plot the point <math>(0,0)</math>. However, since the y-axis is logarithmic, there is no place to plot this on the exponential paper. Also, the failure rate, <math>\lambda </math>, is the negative of the slope of the line, but there is an easier way to determine the value of <math>\lambda </math> from the probability plot, as will be illustrated in the following example.
 
====Example 1====
Six units are put on a life test and tested to failure. The failure times are 7, 12, 19, 29, 41, and 67 hours.  Estimate the failure rate parameter for a one-parameter exponential distribution using the probability plotting method.
 
=====Solution to Example 1=====
In order to plot the points for the probability plot, the appropriate reliability estimate values must be obtained.  These will be equivalent to <math>100%-MR</math>, since the y-axis represents the reliability and the <math>MR</math> values represent unreliability estimates.
 
 
<center><math>\begin{matrix}
  \text{Time-to-} & \text{Reliability}  \\
  \text{failure, hr} & \text{Estimate, }%  \\
  7 & 100-10.91=89.09%  \\
  12 & 100-26.44=73.56%  \\
  19 & 100-42.14=57.86%  \\
  29 & 100-57.86=42.14%  \\
  41 & 100-73.56=26.44%  \\
  67 & 100-89.09=10.91%  \\
\end{matrix}</math></center>
 
 
Next, these points are plotted on exponential probability plotting paper. A sample of this type of plotting paper is shown next, with the sample points in place. Notice how these points describe a line with a negative slope.
Once the points are plotted, draw the best possible straight line through these points. The time value at which this line intersects with a horizontal line drawn at the 36.8% reliability mark is the mean life, and the reciprocal of this is the failure rate <math>\lambda </math>.
This is because at <math>t=m=\tfrac{1}{\lambda }</math>:
 
::<math>\begin{align}
  R(t)= & {{e}^{-\lambda \cdot t}} \\
  R(t)= & {{e}^{-\lambda \cdot \tfrac{1}{\lambda }}} \\
  R(t)= & {{e}^{-1}}=0.368=36.8%. 
\end{align}</math>
 
These steps are shown graphically in the next pages.
 
[[Image:weibullEPP.png|thumb|center|300px|]]
 
As can be seen in the plot below, the best-fit line through the data points crosses the <math>R=36.8%</math> line at <math>t=33</math> hours.
 
[[Image:weibullEPP2.png|thumb|center|300px|]]
 
Since <math>\tfrac{1}{\lambda }=33</math> hours, <math>\lambda =0.0303</math> failures/hour.

Latest revision as of 09:10, 6 August 2012