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_^n f_(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 ''hyper''exponential 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 = \int_^\infty x f(x) \, dx= \sum_^n p_i\int_0^\infty x\lambda_i e^ \, dx = \sum_^n \frac

and

: E\!\left ^2\right= \int_^\infty x^2 f(x) \, dx = \sum_^n p_i\int_0^\infty x^2\lambda_i e^ \, dx = \sum_^n \fracp_i,

from which we can derive the variance:

:\operatorname = E\!\left ^2\right- E\!\left \right2 = \sum_^n \fracp_i - \left sum_^n \frac\right2
= \left sum_^n \frac\right2 + \sum_^n \sum_^n p_i p_j \left(\frac - \frac \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 ^\right= \int_^\infty e^ f(x) \, dx= \sum_^n p_i \int_0^\infty e^\lambda_i e^ \, dx = \sum_^n \fracp_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
* Lomax distribution (continuous mixture of exponentials)

References

{{DEFAULTSORT:Hyperexponential Distribution
Continuous distributions