In this example, we will determine the median rank value used for plotting the sixth failure from a sample size of ten. This will be used to illustrate two of the built-in functions in Weibull++'s Quick Statistical Reference.
First, open the Quick Statistical Reference by clicking its icon.
or by selecting Quick Statistical Reference from the Tools menu.
In this example N = 10, j = 6, m = 2(10 - 6 + 1) = 10, and n = 2 x 6 = 12.
Thus, from the F-distribution rank equation:
Calculate the value of F0.50:10:12 by using the Inverse F-Distribution Values option from the Quick Statistical Reference, or F0.50;10;12 = 0.9886 as shown next:
Consequently:
Another method is to use the Median Ranks option directly, which yields MR(%) = 54.8305%, as shown next:
What is the unreliability of the units in Example 1 for a mission duration of 30 hours, starting the mission at age zero? To replicate the results in this reference with Weibull++, choose RRX (Rank Regression on X) as the calculation method.
First, use Weibull++ to obtain the parameters using RRX.
Then, we investigate several methods of solution for this problem. The first, and more laborious, method is to extract the information directly from the plot. You may do this with either the screen plot in RS Draw or the printed copy of the plot. (When extracting information from the screen plot in RS Draw, note that the translated axis position of your mouse is always shown on the bottom right corner.)
Using this first method, enter either the screen plot or the printed plot with T = 30 hours, go up vertically to the straight line fitted to the data, then go horizontally to the ordinate, and read off . Then, a good estimate of the sought unreliability is 23%. (Also, the reliability estimate is 1.0 - 0.23 = 0.77 or 77%.)
The second method involves the use of the Quick Calculation Pad (QCP). To activate the Quick Calculation Pad either click its icon on the Data Analysis toolbar, or Data Folio Control Panel
or select it from the Tools menu. Once activated, the QCP will appear on your desktop. Switch to the Basic Calculations page of the QCP, select Std. Prob. Calculations from the Options For Calculations box, and enter 30 hours in the Mission End Time box, making sure that the Results as Probability of Failure option is selected in the Results Option box.
Click Calculate to get the result. The result is shown in the figure above, or
Note that the results in QCP vary according to the parameter estimation method used. The above results are obtained using RRX.
What is the reliability for a new mission of t = 10 hours duration, starting the new mission at the age of T = 30 hours, for the same data as Example 7?
The conditional reliability is given by:
or:
Again, the Quick Calculation Pad can provide this result directly and more accurately than the plot.
What is the longest mission that this equipment should undertake for a reliability of 90%?
Using the QCP again, choose Warranty (Time) Information and enter the Required Reliability, 0.90, and click Calculate. The result is 15.9933 hours.
Assume that ten identical units (N = 10) are being reliability tested at the same application and operation stress levels. Six of these units fail during this test after operating the following numbers of hours, Tj: 150, 105, 83, 123, 64 and 46. The test is stopped at the sixth failure. Find the parameters of the Weibull pdf that represents these data.
Open a new Data Folio choosing Times-to-failure data, My data set contains suspensions (right censored data) and I want to enter data in groups.
Enter the data in the appropriate columns. Note that there are four suspensions, as only six of the ten units were tested to failure (the next figure shows the data as entered). Use the three-parameter Weibull and MLE for the calculations.
Plot the data.
Note that the original data points, on the curved line, were adjusted by subtracting 30.92 hours to yield a straight line as shown above.
Suppose that we have run an experiment with eight units being tested and the following is a table of their last inspection times and times-to-failure:
Table 6.5 - The test data for Example 16
|
Data point index |
Last Inspection |
Time-to-failure |
|
1 |
30 |
32 |
|
2 |
32 |
35 |
|
3 |
35 |
37 |
|
4 |
37 |
40 |
|
5 |
42 |
42 |
|
6 |
45 |
45 |
|
7 |
50 |
50 |
|
8 |
55 |
55 |
Analyze the data using several different parameter estimation techniques and compare the results.
This data set can be entered into Weibull++ by opening a new Data Folio and choosing Times-to-failure and My data set contains interval and/or left censored data.
The data is entered as follows,
The computed parameters using maximum likelihood are:
using RRX or rank regression on X:
and using RRY or rank regression on Y:
The plot of the MLE solution with the two-sided 90% confidence bounds is:
ACME company manufactures widgets, and is currently engaged in reliability testing a new widget design. Nineteen units are being reliability tested, but due to the tremendous demand for widgets, units are removed from the test whenever the production cannot cover the demand. The test is terminated at the 67th day when the last widget is removed from the test. Table 6.6 contains the collected data.
Table 6.6 - Widget test data
|
Data point index |
State (S or F) |
Time-to-failure |
|
1 |
F |
2 |
|
2 |
S |
3 |
|
3 |
F |
5 |
|
4 |
S |
7 |
|
5 |
F |
11 |
|
6 |
S |
13 |
|
7 |
S |
17 |
|
8 |
S |
19 |
|
9 |
F |
23 |
|
10 |
F |
29 |
|
11 |
S |
31 |
|
12 |
F |
37 |
|
13 |
S |
41 |
|
14 |
F |
43 |
|
15 |
S |
47 |
|
16 |
S |
53 |
|
17 |
F |
59 |
|
18 |
S |
61 |
|
19 |
S |
67 |
In this example, we see that the number of failures is less than the number of suspensions. This is a very common situation, since reliability tests are often terminated before all units fail due to financial or time constraints. Further, some suspensions will be recorded when a failure occurs that is not due to a legitimate failure mode, such as operator error. In cases such as this, a suspension is recorded, since the unit under test cannot be said to have had a legitimate failure.
This data set can be entered into Weibull++ using Times-to-failure and My data set contains suspensions (right censored data).
We will use the two-parameter Weibull to solve this problem. The parameters using maximum likelihood are:
using RRX:
and using RRY:
Suppose we want to model a left censored, right censored, interval, and complete data set, consisting of 274 units under test of which 185 units fail. Table 6.8 contains the data.
Table 6.8 - The test data for Example 13
|
Data point index |
Number in State |
Last Inspection |
State |
State End Time |
|
1 |
2 |
5 |
F |
5 |
|
2 |
23 |
5 |
S |
5 |
|
3 |
28 |
0 |
F |
7 |
|
4 |
4 |
10 |
F |
10 |
|
5 |
7 |
15 |
F |
15 |
|
6 |
8 |
20 |
F |
20 |
|
7 |
29 |
20 |
S |
20 |
|
8 |
32 |
0 |
F |
22 |
|
9 |
6 |
25 |
F |
25 |
|
10 |
4 |
27 |
F |
30 |
|
11 |
8 |
30 |
F |
35 |
|
12 |
5 |
30 |
F |
40 |
|
13 |
9 |
27 |
F |
45 |
|
14 |
7 |
25 |
F |
50 |
|
15 |
5 |
20 |
F |
55 |
|
16 |
3 |
15 |
F |
60 |
|
17 |
6 |
10 |
F |
65 |
|
18 |
3 |
5 |
F |
70 |
|
19 |
37 |
100 |
S |
100 |
|
20 |
48 |
0 |
F |
102 |
This data set can be entered into Weibull++ by selecting the Times-to-failure and My data set contains suspensions (right censored data), My data set contains interval and/or left censored data and I want to enter data in groups options.
Since standard ranking methods for dealing with these different data types are inadequate, we will want to use the ReliaSoft ranking method. This option is the default in Weibull++ when dealing with interval data.
The computed parameters using MLE are:
using RRX:
and using RRY:
The plot with the two-sided 90% confidence bounds for the rank regression on X solution is:
See Also:
The Weibull Distribution
Go to weibull.com
Go to ReliaSoft.com
©1996-2006. ReliaSoft Corporation. ALL RIGHTS RESERVED.