Engineered Software

Weibull Distribution - Probability Plotting


Techniques

Cumulative Distribution Function

Cumulative Hazard Function

Normal Distribution

Lognormal Distribution

Exponential Distribution

Exam

Engineered Software Home Page

Probability 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 probability plotting for the Weibull distribution using the Reliability & Maintenance Analyst.  The manual method is located here.

Probability 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, select the "Transformed" option in the Plot Type frame, and 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 probability 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 "Probability Plot"

The figure below shows the Weibull probability 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 probability plot with 2 options - transformed data and untransformed data. The transformed data produced a plot of the data that has been transformed to produce a straight line with a 95% confidence interval surrounding it. 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 Probability Plotting

By taking the logarithm of the Weibull cumulative distribution function twice and rearranging,

By plotting versus ln x, and fitting a straight line to the points, the parameters of the Weibull distribution can be estimated. The slope of the plot provides an estimate of b , and the y-intercept can be used to estimate q .

The cumulative distribution function, F(x), is usually estimated from the median rank, but other estimates such as the mean rank and the Kaplan-Meier product limit estimator are also used. Specialized probability paper is available for probability plotting. Using probability paper eliminates the need to transform the data prior to plotting.

Example: Determine the parameters of the Weibull distribution using probability plotting for the data given below.

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

Time to
Fail

Median
Rank,

104 c

         

122

1.1000

1.1000

0.0769

-2.5252

4.8040

138

1.1000

2.2000

0.1827

-1.6008

4.9273

167 c

         

168

1.2571

3.4571

0.3036

-1.0167

5.1240

180

1.2571

4.7143

0.4245

-0.5934

5.1930

217 c

         

229

1.5714

6.2857

0.5755

-0.1544

5.4337

309 c

         

386

2.3571

8.6429

0.8022

0.4827

5.9558

The last two columns of this table are plotted in the figure below. The slope of the best-fit straight line through the data (found using linear regression) is 2.41, which is the estimated value of b. The y-intercept of the best-fit straight line through the data is -13.55. The estimated shape parameter is

Click here to download this solution in Microsoft Excel.

Confidence limits can be added to this plot using 5% and 95% ranks.   Plotting position for the 5% and 95% ranks are found from the expression

BackHomeNext