Example of Parametric Distribution Analysis with right-censored data

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 do not fail before the test ends. Therefore, the data are right censored. The engineer uses Parametric Distribution Analysis 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.MWX.
  2. Choose Stat > Reliability/Survival > Distribution Analysis > Parametric Distribution Analysis.
  3. From Reliability Data Definition, select Custom Definition.
  4. From Censoring, select Right censored data with censoring column.
  5. From Data, select One or more samples, each in a column.
  6. From Failure modes, select None.
  7. In Sample column, enter Temp80.
  8. In Censoring column, enter Cens80.
  9. In Censoring value, select 0.
  10. Select Add Row.
  11. In the second row of Sample column, enter Temp100.
  12. In the second row of Censoring column, enter Cens100.
  13. In the second row of Censoring value, select 0.
  14. From Assumed distribution, select Lognormal.
  15. Select Percentiles/Probabilities. In Estimate percentiles for these additional percents, enter 0.1.
  16. In Estimate, select Survival probabilities for the following times. In Time values, enter 70.
  17. Select OK.
  18. Select Graphs. Select Survival plot.
  19. Select 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 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 time to 0.1% of failures by 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.

Tip

The analysis stores a reliability data definition that includes the sample columns and the censoring scheme. Select the definition to quickly run Distribution ID Plot, Distribution Overview Plot, Parametric Distribution Analysis, and Nonparametric Distribution Analysis on the same 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