Engineered Software

Weibull Distribution - Hazard Plotting


Techniques

Cumulative Distribution Function

Cumulative Hazard Function

Normal Distribution

Lognormal Distribution

Exponential Distribution

Exam

Engineered Software Home Page

Hazard Plotting is tedious and computer routines are often employed. Commercial software is available for these calculations, such as the Reliability & Maintenance Analyst.

The following section describes hazard plotting for the Weibull distribution using the Reliability & Maintenance Analyst.  The manual method is located here.

Hazard plotting supports the 2-parameter and 3-parameter Weibull distribution, and is an excellent method for determining goodness-of-fit. To determine the goodness-of-fit click the "Plot" button. If the plotted points form a straight line, the distribution provides a good time to fail model for the data. The "R-Squared" value is a measure of how well the data forms a straight line. An R-Squared value of 1.0 indicates a perfectly straight line. R-Squared is also known as the coefficient of determination.  To estimate the parameters of the Weibull distribution using hazard plotting, follow these steps:

  1. Enter the data using one of the data entry grids, or connect to a database.
  2. Select the "Parameter Estimation"
  3. Select "Weibull"
  4. Select "Hazard Plot"

The figure below shows the Weibull hazard plotting screen using the data in the file "Demo2.dat".

The default for this operation is a location parameter of zero. In the Location Parameter frame, the software will determine the best value for the location parameter if "Software Estimate" is selected. The user is also given the option of entering a location parameter in the text box below the "User Entered" option. Be cautious when using a non-zero location parameter. A positive location parameter indicates a zero probability of failure for time less than the value of the location parameter, and a negative location parameter means that the population had failures before testing began. Unless there is a sound engineering for one of these conditions, it is best to use a location parameter equal to zero.

Clicking the "Plot" button gives a hazard plot If the plotted points do not follow a straight line, the Weibull distribution with the estimated parameters does not provide an adequate time to fail model. The untransformed option plots the cumulative distribution function against time. This is useful for determining the probability of failure at a given time. The title of the graphs can be changed by editing the text in the Graph Title frame. To check the spelling of the title, click the "Spell Check" button.

To predict reliability or time-to-fail using the estimated parameters use the Predicting Module.

 

Manual Hazard Plotting

The Weibull cumulative hazard function is

H(x) = -ln[1-F(x)]

Replacing F(x) and rearranging gives a linear expression.

By plotting lnH(x) versus lnx the resulting slope (censored points are not plotted) provides an estimate of b . The y-intercept of this plot is an estimate of blnq. Thus, q is estimated from the expression

where, y0 is the y-intercept of the hazard plot.

The hazard function h(x) is estimated from the inverse of the reverse rank of the ordered failures; and the cumulative hazard function, H(x), is the cumulative of the values of h(x). An alternative to plotting lnH(x) versus lnx is to directly plot H(x) versus x on specialized Weibull hazard paper. The advantage of hazard paper is that logarithmic transformations do not have to be computed. Computers have made this technique obsolete.

Example:
Determine the parameters of the Weibull distribution using the multiple censored data in the table below. A "c" following an entry indicates censoring.

Time to Fail

309 c 229
386 104 c
180 217 c
167 c 168
122 138

Solution:  The Table below is constructed to obtain the necessary plotting data.

Time to
Fail

Reverse
Rank

h(t)

H(t)

lnH(t)

lnt

104 c

10

       

122

9

0.1111

0.1111

-2.1972

4.8040

138

8

0.1250

0.2361

-1.4435

4.9273

167 c

7

       

168

6

0.1667

0.4028

-0.9094

5.1240

180

5

0.2000

0.6028

-0.5062

5.1930

217 c

4

       

229

3

0.3333

0.9361

-0.0660

5.4337

309 c

2

       

386

1

1.0000

1.9361

0.6607

5.9558

The final two columns of this table can be plotted. This plot is shown in the figure below. The slope of the best-fit straight line through the data (found by linear regression) is equal to 2.34 and provides an estimate of b. The y-intercept of the best-fit straight line through the data is -13.004. The estimated scale parameter for the Weibull distribution is

Click here to download this solution in Microsoft Excel.

BackHomeNext