Mixture Distributions and Their Actuarial Applications
Understand mixture distributions and how they model heterogeneous insurance portfolios.
Finite Mixtures
A finite mixture distribution has PDF f(x) = w1*f1(x) + w2*f2(x) + ... + wk*fk(x), where the weights wi are non-negative and sum to 1. Each fi represents a component distribution. The weights represent the probability that an observation comes from each component. In insurance, mixtures model heterogeneous risk classes: a portfolio might contain 60% low-risk policies with claim distribution f1 and 40% high-risk policies with distribution f2.
The mean of a mixture is simply the weighted average of the component means: E[X] = w1*mu1 + w2*mu2 + ... + wk*muk. This follows from the law of total expectation, where the conditioning variable is the component membership.
The variance formula is more complex and is critical for Exam P: Var(X) = sum(wi * (sigma_i^2 + mu_i^2)) - (sum(wi * mu_i))^2. This can be rewritten as E[Var(X|Class)] + Var(E[X|Class]), which is the law of total variance. The first term is the weighted average of within-class variances, and the second term captures the variability between class means.
Computing Mixture Probabilities
To compute P(X > c) for a mixture, apply the law of total probability: P(X > c) = w1 * P(X > c | Class 1) + w2 * P(X > c | Class 2) + ... + wk * P(X > c | Class k). Each term uses the CDF or survival function of the corresponding component distribution.
Example: An insurer has two risk classes. Class 1 (70% of policyholders) has Exponential(500) losses. Class 2 (30%) has Exponential(2000) losses. For a randomly selected policyholder, P(X > 1000) = 0.7 * e^(-1000/500) + 0.3 * e^(-1000/2000) = 0.7 * e^(-2) + 0.3 * e^(-0.5) = 0.7(0.1353) + 0.3(0.6065) = 0.0947 + 0.1820 = 0.2767.
Notice that the mixture probability is not simply the probability from a single exponential distribution. The heavier tail of Class 2 contributes disproportionately to the probability of large losses, even though Class 2 is only 30% of the portfolio. This is a key insight about mixtures: they produce heavier tails than any individual component.
Variance of a Two-Component Mixture
For a two-component mixture with weights w and (1-w), means mu1 and mu2, and variances sigma1^2 and sigma2^2, the total variance formula gives:
Var(X) = [w * sigma1^2 + (1-w) * sigma2^2] + [w * mu1^2 + (1-w) * mu2^2 - (w*mu1 + (1-w)*mu2)^2]
The second bracket simplifies to w(1-w)(mu1 - mu2)^2. So Var(X) = w*sigma1^2 + (1-w)*sigma2^2 + w(1-w)*(mu1 - mu2)^2. This clean formula is worth memorizing for two-component mixtures on Exam P.
The term w(1-w)*(mu1 - mu2)^2 shows that the between-class variance increases with the squared difference of the means. When the component means are far apart, the mixture variance can be much larger than either component variance alone.
Continuous Mixtures
In a continuous mixture, a parameter Theta itself has a continuous distribution. Given Theta = theta, X has conditional density f(x|theta). The unconditional (marginal) density of X is f(x) = integral of f(x|theta) * g(theta) d_theta, where g is the PDF of Theta.
The most important continuous mixture for Exam P is the Poisson-Gamma mixture. If X | Lambda ~ Poisson(Lambda) and Lambda ~ Gamma(r, beta), then the marginal distribution of X is Negative Binomial with parameters r and p = beta/(1 + beta). The marginal mean is E[X] = E[Lambda] = r/beta and the marginal variance is Var(X) = E[Lambda] + Var(Lambda) = r/beta + r/beta^2 = r(1 + beta)/beta^2.
Another common continuous mixture: if X | Theta ~ Exponential(Theta) and Theta ~ Gamma(alpha, lambda), then the marginal distribution of X is Pareto-like. Specifically, f(x) = integral from 0 to infinity of theta * e^(-theta*x) * [lambda^alpha * theta^(alpha-1) * e^(-lambda*theta) / Gamma(alpha)] d_theta = alpha * lambda^alpha / (lambda + x)^(alpha + 1), which is a (shifted) Pareto distribution.
Bayes' Theorem and Posterior Classification
A common Exam P problem type combines mixtures with Bayes' theorem. Given that an observation has a particular value (or exceeds a threshold), what is the posterior probability that it came from a specific component?
Using the risk class example: given that a loss exceeds 1000, what is the probability the policyholder is from Class 2? By Bayes' theorem: P(Class 2 | X > 1000) = P(X > 1000 | Class 2) * P(Class 2) / P(X > 1000) = 0.6065 * 0.3 / 0.2767 = 0.1820 / 0.2767 = 0.658. The posterior probability of being high-risk jumped from the prior 0.30 to 0.658 after observing a large loss.
This posterior updating is the foundation of credibility theory, which you will study in later actuarial exams. For Exam P, focus on computing posterior probabilities correctly using Bayes' theorem with the mixture structure.
Mixture Models and Insurance Pricing
Mixture distributions arise naturally in insurance because policyholders are not homogeneous. Even within a single rating class, there is unobservable heterogeneity in risk levels. The mixing distribution (e.g., Gamma for the Poisson rate) captures this heterogeneity.
A practical consequence is that mixture distributions always have a variance at least as large as the average within-class variance. The additional between-class variance represents the impact of heterogeneity. If an insurer could perfectly classify all policyholders, the between-class variance would be eliminated and premiums could be set more precisely.
On Exam P, mixture problems test your ability to apply the total expectation and total variance formulas, compute mixture probabilities, and use Bayes' theorem for posterior classification. These skills combine many of the probability tools covered in the exam syllabus.
Common Pitfalls
A frequent mistake is computing the variance of a mixture as the weighted average of the component variances: w1*sigma1^2 + w2*sigma2^2. This gives only the within-class variance E[Var(X|Class)] and omits the between-class variance Var(E[X|Class]). The correct formula always includes both terms.
Another error is confusing the mixture density with the conditional density. The mixture f(x) = w1*f1(x) + w2*f2(x) is the unconditional density. The conditional density given class membership is simply f_i(x). Using f(x) when you should use f_i(x) (or vice versa) leads to incorrect probabilities.
Students sometimes forget that the CDF of a mixture is also a weighted average: F(x) = w1*F1(x) + w2*F2(x). This follows directly from integrating the mixture density and makes computing probabilities straightforward.
Study Tips
For Exam P, practice the complete workflow: compute mixture mean and variance, find probabilities for the mixture, and apply Bayes' theorem for posterior classification. Memorize the two-component variance formula Var(X) = w*sigma1^2 + (1-w)*sigma2^2 + w(1-w)*(mu1 - mu2)^2. Know the Poisson-Gamma conjugate pair and its Negative Binomial result. Work through at least five problems that combine mixtures with Bayes' theorem, as this is a common exam pattern. These problems typically require four to five minutes and test multiple concepts simultaneously.