Triangular Distribution
Triangular Distribution
Carl Sorensen
April 15, 2011
It is sometimes difficult to get good data on setup and repair times. Ideally,
a large number of events (either setup or repair) would be recorded, and an
average and standard deviation could be calculated from the recorded data.
Additionally, a histogram of the times could be plotted and the general shape
of the data could be recorded.
Rare events like breakdowns, however, are difficult to measure in this way.
The distributions tend not to be normal, and the effects of rare but very long
repairs are likely to either not show up in the sample, or to overly bias the
sample if they do show up.
One way around this difficulty is to use the triangular distribution to model
the repair times. The triangular distribution has three characteristic values:
a minimum, a maximum, and a mode. The probability density increases in a
straight line from the minimum to the mode, then decreases in a straight line
to the maximum. The height of the mode is adjusted to make the area under
the curve be 1.
Probability Density 6
2 2
xmax −xmin = δ3
QQ
Q
Q
Q
Q
Q
Q
Q -
xmin xmode xmax
Random Value
δ1 = xmode − xmin
δ2 = xmax − xmode (1)
δ3 = xmax − xmin
1
A graph of the triangular distribution is shown in Figure 1. The probability
density of the triangular distribution is given by
( 2(x−xmin ) 2(x−xmin )
(xmode −xmin )(xmax −xmin ) = δ1 δ3 for xmin < x < xmode
y= 2(xmax −x) 2(xmax −x) (2)
(xmax −xmode )(xmax −xmin ) = δ2 δ3 for xmode < x < xmax
h(x − min)
p(x) = (7)
δ1
To the right of the mode, the density function is given by
h(max − x)
p(x) = (8)
δ2
2
This gives us the value of the density function as
(
2(x−xmin )
δ1 δ3 ; xmin ≤ x ≤ xmode
p(x) = 2(xmax −x) (9)
δ2 δ3 ; xmode ≤ x ≤ xmax
The cumulative distribution function G(x) can be determined by integrating
the density function p(x):
x τ Rx
; xmin ≤ x ≤ xmode
Z Z
xRmin
p(τ )dτ
G(x) = p(τ )dτ = p(τ )dτ = xmax
−∞ xmin 1− x
; xmode ≤ x ≤ xmax
p(τ )dτ
(10)
When combined with the expression for p(x) given above, we obtain:
(x−xmin )2
(
δ1 δ3 ; xmin ≤ x ≤ xmode
G(x) = (xmax −x)2
(11)
1− δ2 δ3 ; xmode ≤ x ≤ xmax
x4mode x4 x3
4
x4 x3
2 xmin 2 xmode xmax
µ2 = + min − mode + + max − mode
δ1 δ3 4 12 3 δ2 δ3 4 12 3
(13)
An easier way to evaluate this integral is to recognize that
µ2 = E(x2 ) = σ 2 + µ2 (15)
and the previous equations for the mean and standard deviation can be used.