Article 4SDR9 Hyperexponential and hypoexponential distributions

Hyperexponential and hypoexponential distributions

by
John
from John D. Cook on (#4SDR9)

There are a couple different ways to combine random variables into a new random variable: means and mixtures. To take the mean of X and Y you average their values. To take the mixture of X and Y you average their densities. The former makes the tails thinner. The latter makes the tails thicker. When X and Y are exponential random variables, the mean has a hypoexponential distribution and the mixture has a hyperexponential distribution.

Hypoexponential distributions

Suppose X and Y are exponentially distributed with mean I1/4. Then their sum X + Y has a gamma distribution with shape 2 and scale I1/4. The sum has mean 2I1/4 and variance 2I1/4^2. The coefficient of variation, the ratio of the standard deviation to the mean, is 1/a2. The hypoexponential distribution is so-called because its coefficient of variation is less than 1, whereas an exponential distribution has coefficient of variation 1 because the mean and standard deviation are always the same.

The means of X and Y don't have to be the same. When they're different, the sum does not have a gamma distribution, and so hypoexponential distributions are more general than gamma distributions.

A hypoexponential random variable can also be the sum of more than two exponentials. If it is the sum of n exponentials, each with the same mean, then the coefficient of variation is 1/an. In general, the coefficient of variation for a hypoexponential distribution is the coefficient of variation of the means [1].

In the introduction we talked about means rather than sums, but it makes no difference to the coefficient of variation because going from sum to mean divides the mean and the standard deviation by the same amount.

Hyperexponential distributions

Hyperexponential random variables are constructed as a mixture of exponentials rather than an average. Instead selecting a value from X and a value from Y, we select a value from X or a value from Y. Given a mixture probability p, we choose a sample from X with probability p and a value from Y with probability 1 - p.

The density function for a mixture is a an average of the densities of the two components. So if X and Y have density functions fX and fY, then the mixture has density

p fX + (1 - p) fY

If you have more than two random variables, the distribution of their mixture is a convex combination of their individual densities. The coefficients in the convex combination are the probabilities of selecting each random variable.

If X and Y are exponential with means I1/4X and I1/4Y, and we have a mixture that selects X with probability p, then then mean of the mixture is the mixture of the means

I1/4 = p I1/4X + (1 - p) I1/4Y

which you might expect, but the variance

If^2 = p I1/4X ^2 + (1 - p) I1/4Y ^2 + p(1 - p)(I1/4X - I1/4Y)^2

is not quite analogous because of the extra p(1 - p)(I1/4X - I1/4Y)^2 term at the end. If I1/4X = I1/4Y this last term drops out and the coefficient of variation is 1: mixing two identically distributed random variables doesn't do anything to the distribution. But when the means are different, the coefficient of variation is greater than 1 because of the extra term in the variance of the mixture.

Example

Suppose I1/4X = 2 and I1/4Y = 8. Then the average of X and Y has mean 5, and so does an equal mixture of X and Y.

The average of X and Y has standard deviation a17, and so the coefficient of variation is a17/5 = 0.8246.

An exponential distribution with mean 5 would have standard deviation 5, and so the coefficient of variation 1.

An equal mixture of X and Y has standard deviation a43, and so the coefficient of variation is a43/5 = 1.3114.

Related posts

[1] If exponential random variables Xi have means I1/4i, then the coefficient of variation of their sum (or average) is

a(I1/41^2 + I1/42^2 + " + I1/4n^2) / (I1/41 + I1/42 + " + I1/4n)

q5SsIUNXUjY
External Content
Source RSS or Atom Feed
Feed Location http://feeds.feedburner.com/TheEndeavour?format=xml
Feed Title John D. Cook
Feed Link https://www.johndcook.com/blog
Reply 0 comments