Example of Parametric Distribution Analysis (Right Censoring)

A reliability engineer studies the failure rates of engine windings of turbine assemblies to determine the times at which the windings fail. At high temperatures, the windings might decompose too fast.

The engineer records failure times for the engine windings at 80° C and 100° C. However, some of the units must be removed from the test before they fail. Therefore, the data are right censored. The engineer uses Parametric Distribution Analysis (Right Censoring) to determine the following:
  • The times at which various percentages of the windings fail. The engineer is particularly interested in the 0.1th percentile
  • The percentage of windings that will survive past 70 hours
  • The survival function for the engine windings (as shown on a survival plot)
  • The fit of the lognormal distribution for the data (as shown on a probability plot)
  1. Open the sample data, EngineWindingReliability.MTW.
  2. Choose Stat > Reliability/Survival > Distribution Analysis (Right Censoring) > Parametric Distribution Analysis.
  3. In Variables, enter Temp80 Temp100.
  4. From Assumed distribution, select Lognormal.
  5. Click Censor. Under Use censoring columns, enter Cens80 Cens100.
  6. In Censoring value, type 0. Click OK.
  7. Click Estimate. In Estimate percentiles for these additional percents, enter 0.1.
  8. In Estimate probabilities for these times (values), enter 70. Click OK.
  9. Click Graphs. Select Survival plot.
  10. Click OK in each dialog box.

Interpret the results

Using the Table of Percentiles, the engineer can determine the times at which various percentages of the windings fail. At 80° C, 1% of the windings to fail by 19.3281 hours. The values for the 0.1th percentile, which the engineer requested for the analysis, are also shown in the table. At 80° C, 0.1% of the windings fail by 13.3317 hours. At 100° C, 0.1% of the windings fail by 3.93505 hours. Therefore, the increase in temperature decreases the percentile by a value of approximately 9.5 hours.

Using the Table of Survival Probabilities, the engineer can determine what proportion of windings are expected to survive for more than 70 hours. At 80° C, 37.43% of the windings are expected to survive for more than 70 hours. At 100° C, 19.82% of the windings are expected to survive for more than 70 hours.

The engineer uses the survival plot to view the survival probabilities over time, and the probability plot to check that the lognormal distribution adequately fits the data.

80° C
Variable: Temp80

Censoring

Censoring InformationCount
Uncensored value37
Right censored value13
Censoring value: Cens80 = 0
Estimation Method: Maximum Likelihood
Distribution:   Lognormal

Parameter Estimates



Standard
Error
95.0% Normal CI
ParameterEstimateLowerUpper
Location4.092670.07196813.951614.23372
Scale0.4862160.06062470.3807990.620816
Log-Likelihood = -181.625

Goodness-of-Fit

Anderson-Darling
(Adjusted)
67.800

Characteristics of Distribution



Standard
Error
95.0% Normal CI

EstimateLowerUpper
Mean(MTTF)67.41535.5524557.365679.2255
Standard Deviation34.81456.7982723.743551.0476
Median59.89954.3108552.019268.9735
First Quartile(Q1)43.15163.2952637.153150.1186
Third Quartile(Q3)83.14757.3769069.876398.9392
Interquartile Range(IQR)39.99596.3331729.324554.5505

Table of Percentiles



Standard
Error
95.0% Normal CI
PercentPercentileLowerUpper
0.113.33172.515599.2102619.2975
119.32812.8375014.495325.7722
222.06742.9255917.017828.6154
324.00342.9726118.830430.5975
425.57093.0035520.312632.1906
526.92123.0262121.597833.5566
628.12653.0440322.750634.7727
729.22763.0588123.807435.8819
830.25013.0716524.791036.9113
931.21103.0832625.717037.8788
1032.12253.0940926.596238.7970
2039.78373.2099733.964646.5999
3046.41843.4101540.193653.6073
4052.95733.7566946.083360.8568
5059.89954.3108552.019268.9735
6067.75175.1591058.358478.6569
7077.29586.4592065.618491.0514
8090.18638.5821174.8412108.678
90111.69612.810389.2100139.849
91114.95813.511291.3052144.738
92118.61014.312093.6288150.255
93122.75915.241796.2426156.581
94127.56516.343799.2372163.979
95133.27617.6863102.753172.866
96140.31419.3873107.026183.955
97149.47721.6739112.500198.608
98162.59025.0764120.175219.977
99185.63431.3868133.271258.570

Table of Survival Probabilities



95.0% Normal CI
TimeProbabilityLowerUpper
700.3742990.2631020.497141
100° C
Variable: Temp100

Censoring

Censoring InformationCount
Uncensored value34
Right censored value6
Censoring value: Cens100 = 0
Estimation Method: Maximum Likelihood
Distribution:   Lognormal

Parameter Estimates



Standard
Error
95.0% Normal CI
ParameterEstimateLowerUpper
Location3.628690.1177853.397843.85955
Scale0.7309390.09198080.5711720.935397
Log-Likelihood = -160.688

Goodness-of-Fit

Anderson-Darling
(Adjusted)
17.253

Characteristics of Distribution



Standard
Error
95.0% Normal CI

EstimateLowerUpper
Mean(MTTF)49.19696.9176137.346564.8076
Standard Deviation41.343111.041624.494769.7806
Median37.66364.4362029.899547.4439
First Quartile(Q1)23.00442.9505517.891029.5791
Third Quartile(Q3)61.66438.4984347.067780.7876
Interquartile Range(IQR)38.66007.2449526.775955.8185

Table of Percentiles



Standard
Error
95.0% Normal CI
PercentPercentileLowerUpper
0.13.935051.172892.194017.05767
16.877641.616984.3382710.9034
28.394101.794205.5212112.7619
39.525281.911136.4282714.1144
410.47562.001467.2036015.2338
511.31812.076587.8995416.2162
612.08842.141878.5418417.1076
712.80692.200319.1453517.9343
813.48632.253769.7194918.7129
914.13542.3034410.270719.4544
1014.76062.3502510.803620.1667
2020.35892.7525615.619726.5362
3025.67173.1661920.159232.6916
4031.29673.6949624.831639.4451
5037.66364.4362029.899547.4439
6045.32585.5315835.683257.5740
7055.25727.2444742.735971.4473
8069.676910.205452.289692.8456
9096.104016.696868.3686135.091
91100.35417.842070.8271142.191
92105.18519.172773.5864150.351
93110.76520.746476.7308159.894
94117.34722.650280.3853171.305
95125.33425.024284.7457185.362
96135.41428.114190.1452203.417
97148.92532.405097.2189228.130
98168.99339.0628107.427265.843
99206.25552.1976125.600338.704

Table of Survival Probabilities



95.0% Normal CI
TimeProbabilityLowerUpper
700.1982330.1071870.324816