8.1.6.3. Extreme value distributions

8. Assessing Product Reliability
8.1. Introduction
8.1.6. What are the basic lifetime distribution models used for non-repairable populations?

8.1.6.3. Extreme value distributions

The Extreme Value Distribution usually refers to the distribution of the minimum of a large number of unbounded random observations

Description, Formulas and Plots

We have already referred to Extreme Value Distributions when describing the uses of the Weibull distribution. Extreme value distributions are the limiting distributions for the minimum or the maximum of a very large collection of random observations from the same arbitrary distribution. Gumbel (1958) showed that for any well-behaved initial distribution (i.e., F(x) is continuous and has an inverse), only a few models are needed, depending on whether you are interested in the maximum or the minimum, and also if the observations are bounded above or below.

In the context of reliability modeling, extreme value distributions for the minimum are frequently encountered. For example, if a system consists of n identical components in series, and the system fails when the first of these components fails, then system failure times are the minimum of n random component failure times. Extreme value theory says that, independent of the choice of component model, the system model will approach a Weibull as n becomes large. The same reasoning can also be applied at a component level, if the component failure occurs when the first of many similar competing failure processes reaches a critical level.

The distribution often referred to as the Extreme Value Distribution (Type I) is the limiting distribution of the minimum of a large number of unbounded identically distributed random variables. The PDF and CDF are given by:

Extreme Value Distribution formulas and PDF shapes

f(x) = (1/beta)*EXP((x-mu)/beta)*EXP(-EXP((x-mu)/beta), -inf < x < inf; F(x) = 1 - EXP(-EXP((x-mu)/beta), -inf < x < inf

If the x values are bounded below (as is the case with times of failure) then the limiting distribution is the Weibull. Formulas and uses of the Weibull have already been discussed.

PDF Shapes for the (minimum) Extreme Value Distribution (Type I) are shown in the following figure.

Sample plots for the extreme value PDF

The natural log of Weibull data is extreme value data

Uses of the Extreme Value Distribution Model

In any modeling application for which the variable of interest is the minimum of many random factors, all of which can take positive or negative values, try the extreme value distribution as a likely candidate model. For lifetime distribution modeling, since failure times are bounded below by zero, the Weibull distribution is a better choice.
The Weibull distribution and the extreme value distribution have a useful mathematical relationship. If t₁, t₂, ...,t_n are a sample of random times of fail from a Weibull distribution, then ln t₁, ln t₂, ...,ln t_n are random observations from the extreme value distribution. In other words, the natural log of a Weibull random time is an extreme value random observation.

If the Weibull has the shape parameter gamma and the characteristic life alpha, then the extreme value distribution (after taking natural logarithms) has mu = LN(alpha), beta=1/gamma

Because of this relationship, computer programs and graph papers designed for the extreme value distribution can be used to analyze Weibull data. The situation exactly parallels using normal distribution programs to analyze lognormal data, after first taking natural logarithms of the data points.

Dataplot commands for the extreme value distribution

DATAPLOT for the Extreme Value Distribution

Assume µ = ln 200,000 = 12.206 and = 1/2 = .5. The extreme value distribution with these parameters could be obtained by taking natural logarithms of data from a Weibull population with characteristic life = 200,000 and shape = 2. We will use Dataplot to evaluate PDF's, CDF's and generate random numbers from this distribution. Note that you must first set MINMAX to 1 in order to do (minimum) extreme value type I calculations.

SET MINMAX 1
LET BET = .5
LET M = LOG(200000)
LET X = DATA 5 8 10 12 12.8
LET PD = EV1PDF(X, M, BET)
LET CD = EV1CDF(X, M, BET)

Dataplot will calculate PDF and CDF values corresponding to the points 5, 8, 10, 12, 12.8. The PDF's are .110E^-5, .444E^-3, .024, .683 and .247. The CDF's are .551E^-6, .222E^-3, .012, .484 and .962.

Finally, we generate 100 random numbers from this distribution and construct an extreme value distribution probability plot as follows:

LET SAM = EXTREME VALUE TYPE 1 RANDOM NUMBERS FOR I = 1 1 100
LET SAM = (BET*SAMPLE) + M
EXTREME VALUE TYPE 1 PROBABILITY PLOT SAM

Extreme value type I probability plot for 100 extreme value type I random numbers

Data from an extreme value distribution will line up approximately along a straight line when this kind of plot is constructed. The slope of the line is an estimate of , and the "y-axis" value on the line corresponding to the "x-axis" 0 point is an estimate of µ. For the graph above, these turn out to be very close to the actual values of and µ.