Weibull Distribution Characteristics: Difference between revisions
Line 29: | Line 29: | ||
[[Image:WB.8 effect of weibull.png|center|400px| Effect on <math>\beta</math> on the <math>cdf</math> on the Weibull probability plot with a fixed value of <math>\eta</math> ]] | [[Image:WB.8 effect of weibull.png|center|400px| Effect on <math>\beta</math> on the <math>cdf</math> on the Weibull probability plot with a fixed value of <math>\eta</math> ]] | ||
The above figure shows the effect of the value of < | The above figure shows the effect of the value of <math>\beta\,\!</math> on the <math>cdf</math>, as manifested in the Weibull probability plot. It is easy to see why this parameter is sometimes referred to as the slope. Note that the models represented by the three lines all have the same value of <math>\eta\,\!</math>. The following figure shows the effects of these varied values of <math>\beta\,\!</math> on the reliability plot, which is a linear analog of the probability plot. | ||
[[Image:WB.8 weibull reliability.png|center|450px| The effect of values of <math>\beta</math> on the Weibull reliability plot. ]] | [[Image:WB.8 weibull reliability.png|center|450px| The effect of values of <math>\beta</math> on the Weibull reliability plot. ]] | ||
:*< | :*<math>R(t)\,\!</math> decreases sharply and monotonically for <math>0 < \beta < 1\,\!</math> and is convex. | ||
:*For < | :*For <math>\beta = 1\,\!</math>, <math>R(t)\,\!</math> decreases monotonically but less sharply than for <math>0 < \beta < 1\,\!</math> and is convex. | ||
:*For < | :*For <math>\beta > 1\,\!</math>, <math>R(t)\,\!</math> decreases as increases. As wear-out sets in, the curve goes through an inflection point and decreases sharply. | ||
'''The Effect of <span class="texhtml">β</span> on the Weibull Failure Rate''' | '''The Effect of <span class="texhtml">β</span> on the Weibull Failure Rate''' |
Revision as of 02:22, 22 August 2012
This article also appears in the Life Data Analysis Reference and Accelerated Life Testing Data Analysis Reference books.
The Weibull distribution is widely used in reliability and life data analysis due to its versatility. Depending on the values of the parameters, the Weibull distribution can be used to model a variety of life behaviors. We will now examine how the values of the shape parameter, [math]\displaystyle{ \beta\,\! }[/math], and the scale parameter, [math]\displaystyle{ \eta\,\! }[/math], affect such distribution characteristics as the shape of the curve, the reliability and the failure rate. Note that in the rest of this section we will assume the most general form of the Weibull distribution, (i.e., the 3-parameter form). The appropriate substitutions to obtain the other forms, such as the 2-parameter form where [math]\displaystyle{ \gamma = 0,\,\! }[/math] or the 1-parameter form where [math]\displaystyle{ \beta = C = \,\! }[/math] constant, can easily be made.
Effects of the Shape Parameter, β
The Weibull shape parameter, [math]\displaystyle{ \beta\,\! }[/math], is also known as the slope. This is because the value of [math]\displaystyle{ \beta\,\! }[/math] is equal to the slope of the regressed line in a probability plot. Different values of the shape parameter can have marked effects on the behavior of the distribution. In fact, some values of the shape parameter will cause the distribution equations to reduce to those of other distributions. For example, when [math]\displaystyle{ \beta = 1\,\! }[/math], the of the 3-parameter Weibull reduces to that of the 2-parameter exponential distribution or:
- [math]\displaystyle{ f(t)={\frac{1}{\eta }}e^{-{\frac{t-\gamma }{\eta }}} }[/math]
where [math]\displaystyle{ \frac{1}{\eta }=\lambda = \,\! }[/math] failure rate. The parameter [math]\displaystyle{ \beta\,\! }[/math] is a pure number, (i.e., it is dimensionless). The following figure shows the effect of different values of the shape parameter, [math]\displaystyle{ \beta\,\! }[/math], on the shape of the [math]\displaystyle{ pdf }[/math]. As you can see, the shape can take on a variety of forms based on the value of [math]\displaystyle{ \beta\,\! }[/math].
For [math]\displaystyle{ 0\lt \beta \leq 1 \,\! }[/math]:
- As [math]\displaystyle{ t \rightarrow 0\,\! }[/math] (or [math]\displaystyle{ \gamma\,\! }[/math]), [math]\displaystyle{ f(t)\rightarrow \infty.\,\! }[/math]
- As [math]\displaystyle{ t\rightarrow \infty\,\! }[/math], [math]\displaystyle{ f(t)\rightarrow 0 }[/math].
- [math]\displaystyle{ f(t)\,\! }[/math] decreases monotonically and is convex as it increases beyond the value of [math]\displaystyle{ \gamma\,\! }[/math].
- The mode is non-existent.
For [math]\displaystyle{ \beta \gt 1 \,\! }[/math]:
- [math]\displaystyle{ f(t) = 0\,\! }[/math] at [math]\displaystyle{ t = 0\,\! }[/math] (or [math]\displaystyle{ \gamma\,\! }[/math]).
- [math]\displaystyle{ f(t)\,\! }[/math] increases as [math]\displaystyle{ t\rightarrow \tilde{T} \,\! }[/math] (the mode) and decreases thereafter.
- For [math]\displaystyle{ \beta \lt 2.6\,\! }[/math] the Weibull [math]\displaystyle{ pdf }[/math] is positively skewed (has a right tail), for [math]\displaystyle{ 2.6 \lt \beta \lt 3.7\,\! }[/math] its coefficient of skewness approaches zero (no tail). Consequently, it may approximate the normal [math]\displaystyle{ pdf }[/math] , and for [math]\displaystyle{ \beta \gt 3.7\,\! }[/math] it is negatively skewed (left tail). The way the value of [math]\displaystyle{ \beta\,\! }[/math] relates to the physical behavior of the items being modeled becomes more apparent when we observe how its different values affect the reliability and failure rate functions. Note that for [math]\displaystyle{ \beta = 0.999\,\! }[/math], [math]\displaystyle{ f(0) = \infty\,\! }[/math], but for [math]\displaystyle{ \beta = 1.001\,\! }[/math], [math]\displaystyle{ f(0) = 0.\,\! }[/math] This abrupt shift is what complicates MLE estimation when [math]\displaystyle{ \beta\,\! }[/math] is close to 1.
The Effect of β on the [math]\displaystyle{ cdf }[/math] and Reliability Function
The above figure shows the effect of the value of [math]\displaystyle{ \beta\,\! }[/math] on the [math]\displaystyle{ cdf }[/math], as manifested in the Weibull probability plot. It is easy to see why this parameter is sometimes referred to as the slope. Note that the models represented by the three lines all have the same value of [math]\displaystyle{ \eta\,\! }[/math]. The following figure shows the effects of these varied values of [math]\displaystyle{ \beta\,\! }[/math] on the reliability plot, which is a linear analog of the probability plot.
- [math]\displaystyle{ R(t)\,\! }[/math] decreases sharply and monotonically for [math]\displaystyle{ 0 \lt \beta \lt 1\,\! }[/math] and is convex.
- For [math]\displaystyle{ \beta = 1\,\! }[/math], [math]\displaystyle{ R(t)\,\! }[/math] decreases monotonically but less sharply than for [math]\displaystyle{ 0 \lt \beta \lt 1\,\! }[/math] and is convex.
- For [math]\displaystyle{ \beta \gt 1\,\! }[/math], [math]\displaystyle{ R(t)\,\! }[/math] decreases as increases. As wear-out sets in, the curve goes through an inflection point and decreases sharply.
The Effect of β on the Weibull Failure Rate
The value of β has a marked effect on the failure rate of the Weibull distribution and inferences can be drawn about a population's failure characteristics just by considering whether the value of β is less than, equal to, or greater than one.
As indicated by above figure, populations with β < 1 exhibit a failure rate that decreases with time, populations with β = 1 have a constant failure rate (consistent with the exponential distribution) and populations with β > 1 have a failure rate that increases with time. All three life stages of the bathtub curve can be modeled with the Weibull distribution and varying values of β. The Weibull failure rate for 0 < β < 1 is unbounded at (or γ). The failure rate, λ(t), decreases thereafter monotonically and is convex, approaching the value of zero as t→∞ or λ(∞) = 0. 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 encountering such behavior in a manufactured product, it may be indicative of problems in the production process, inadequate burn-in, substandard parts and components, or problems with packaging and shipping. For β = 1, λ(t) yields a constant value of [math]\displaystyle{ { \frac{1}{\eta }} }[/math] or:
- [math]\displaystyle{ \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 β > 1, λ(t) increases as t increases and becomes suitable for representing the failure rate of units exhibiting wear-out type failures. For 1 < β < 2, the λ(t) curve is concave, consequently the failure rate increases at a decreasing rate as t increases.
For β = 2 there emerges a straight line relationship between λ(t) and t, starting at a value of λ(t) = 0 at t = γ, and increasing thereafter with a slope of [math]\displaystyle{ { \frac{2}{\eta ^{2}}} }[/math]. Consequently, the failure rate increases at a constant rate as t increases. Furthermore, if η = 1 the slope becomes equal to 2, and when γ = 0, λ(t) becomes a straight line which passes through the origin with a slope of 2. Note that at β = 2, the Weibull distribution equations reduce to that of the Rayleigh distribution.
When β > 2, the λ(t) curve is convex, with its slope increasing as t increases. Consequently, the failure rate increases at an increasing rate as t increases, indicating wearout life.
Effects of the Scale Parameter, η
A change in the scale parameter η has the same effect on the distribution as a change of the abscissa scale. Increasing the value of η while holding β constant has the effect of stretching out the [math]\displaystyle{ pdf }[/math]. Since the area under a [math]\displaystyle{ pdf }[/math] curve is a constant value of one, the "peak" of the [math]\displaystyle{ pdf }[/math] curve will also decrease with the increase of η, as indicated in the above figure.
- If η is increased while β and γ are kept the same, the distribution gets stretched out to the right and its height decreases, while maintaining its shape and location.
- If η is decreased while β and γ are kept the same, the distribution gets pushed in towards the left (i.e., towards its beginning or towards 0 or γ), and its height increases.
- η has the same units as [math]\displaystyle{ t }[/math], such as hours, miles, cycles, actuations, etc.
Effects of the Location Parameter, γ
The location parameter, γ, as the name implies, locates the distribution along the abscissa. Changing the value of γ has the effect of sliding the distribution and its associated function either to the right (if γ > 0) or to the left (if γ < 0).
- When γ = 0, the distribution starts at [math]\displaystyle{ t=0 }[/math] or at the origin.
- If γ > 0, the distribution starts at the location γ to the right of the origin.
- If γ < 0, the distribution starts at the location γ to the left of the origin.
- γ provides an estimate of the earliest time-to-failure of such units.
- The life period 0 to + γ is a failure free operating period of such units.
- The parameter γ may assume all values and provides an estimate of the earliest time a failure may be observed. A negative γ may indicate that failures have occurred prior to the beginning of the test, namely during production, in storage, in transit, during checkout prior to the start of a mission, or prior to actual use.
- γ has the same units as [math]\displaystyle{ t }[/math], such as hours, miles, cycles, actuations, etc.