Hyperexponential distribution

Diagram showing queueing system equivalent of a hyperexponential distribution

In probability theory, a hyperexponential distribution is a continuous probability distribution whose probability density function of the random variable X is given by

$f_X(x) = \sum_{i=1}^n f_{Y_i}(x)\;p_i,$

where each Y_i is an exponentially distributed random variable with rate parameter λ_i, and p_i is the probability that X will take on the form of the exponential distribution with rate λ_i. It is named the hyperexponential distribution since its coefficient of variation is greater than that of the exponential distribution, whose coefficient of variation is 1, and the hypoexponential distribution, which has a coefficient of variation smaller than one. While the exponential distribution is the continuous analogue of the geometric distribution, the hyperexponential distribution is not analogous to the hypergeometric distribution. The hyperexponential distribution is an example of a mixture density.

An example of a hyperexponential random variable can be seen in the context of telephony, where, if someone has a modem and a phone, their phone line usage could be modeled as a hyperexponential distribution where there is probability p of them talking on the phone with rate λ₁ and probability q of them using their internet connection with rate λ₂.

Properties

Since the expected value of a sum is the sum of the expected values, the expected value of a hyperexponential random variable can be shown as

$E[X] = \int_{-\infty}^\infty x f(x) \, dx= \sum_{i=1}^n p_i\int_0^\infty x\lambda_i e^{-\lambda_ix} \, dx = \sum_{i=1}^n \frac{p_i}{\lambda_i}$

and

$E\!\left[X^2\right] = \int_{-\infty}^\infty x^2 f(x) \, dx = \sum_{i=1}^n p_i\int_0^\infty x^2\lambda_i e^{-\lambda_ix} \, dx = \sum_{i=1}^n \frac{2}{\lambda_i^2}p_i,$

from which we can derive the variance:

$\operatorname{Var}[X] = E\!\left[X^2\right] - E\!\left[X\right]^2 = \sum_{i=1}^n \frac{2}{\lambda_i^2}p_i - \left[\sum_{i=1}^n \frac{p_i}{\lambda_i}\right]^2 = \left[\sum_{i=1}^n \frac{p_i}{\lambda_i}\right]^2 + \sum_{i=1}^n \sum_{j=1}^n p_i p_j \left(\frac{1}{\lambda_i} - \frac{1}{\lambda_j} \right)^2.$

The standard deviation exceeds the mean in general (except for the degenerate case of all the λs being equal), so the coefficient of variation is greater than 1.

The moment-generating function is given by

$E\!\left[e^{tx}\right] = \int_{-\infty}^\infty e^{tx} f(x) \, dx= \sum_{i=1}^n p_i \int_0^\infty e^{tx}\lambda_i e^{-\lambda_i x} \, dx = \sum_{i=1}^n \frac{\lambda_i}{\lambda_i - t}p_i.$

Fitting

A given probability distribution, including a heavy-tailed distribution, can be approximated by a hyperexponential distribution by fitting recursively to different time scales using Prony's method.

See also

• Phase-type distribution
• Hyper-Erlang distribution