Why Covariance Matters in a Normal Combination
A sum or difference of random variables often represents a quantity of direct interest: a portfolio total, a difference between two measurements, or a weighted combination of risks. If the variables are correlated, their dependence affects the variability of that result. Ignoring this dependence can substantially misstate a probability, even when each variable separately has a normal distribution.
The key setting for an exact normal calculation is that the variables are jointly normal. Under that assumption, every linear combination is normal, and its variance is determined by both the individual variances and the covariances. This tutorial develops the calculation for two variables and for a vector of variables, then applies it to probability questions. It is aligned with the SOA Exam P (Probability) syllabus.
For two variables, the covariance is \(\operatorname{Cov}(X,Y)=E[(X-E[X])(Y-E[Y])]\). If the standard deviations are \(\sigma_X\) and \(\sigma_Y\), their correlation is \(\rho_{X,Y}=\operatorname{Cov}(X,Y)/(\sigma_X\sigma_Y)\), so \(\operatorname{Cov}(X,Y)=\rho_{X,Y}\sigma_X\sigma_Y\). Correlation is dimensionless; covariance carries the product of the variables’ units.
Distribution of a Linear Combination
Suppose \(X\) and \(Y\) are jointly normal, with \(X\sim N(\mu_X,\sigma_X^2)\), \(Y\sim N(\mu_Y,\sigma_Y^2)\), and covariance \(\operatorname{Cov}(X,Y)\). For fixed constants \(a\) and \(b\), the combination \(W=aX+bY\) is normal. Its mean and variance are:
The sign of the covariance contribution depends on the coefficients as well as on the covariance itself. For a difference \(X-Y\), the coefficient product is negative, so a positive covariance reduces the variance of the difference. A negative covariance instead increases it.
The final step uses the definition of covariance. Joint normality then supplies the additional fact that \(aX+bY\) has a normal distribution with this mean and variance. The variance identity itself does not require normality; the claim that the combination is normal does.
Worked Example: A Weighted Sum With Positive Dependence
In a hypothetical model, \(X\) and \(Y\) are jointly normal with \(X\sim N(8,4)\), \(Y\sim N(5,9)\), and correlation \(1/3\). Find the distribution of \(T=2X-Y\) and calculate \(P(T>15)\).
Method. First convert the correlation to covariance, then apply the linear-combination formula and standardize. The model is jointly normal, so the resulting distribution is exactly normal.
The standard deviations are \(\sigma_X=\sqrt{4}=2\) and \(\sigma_Y=\sqrt{9}=3\). Therefore,
The mean of \(T\) is \(2(8)-5=11\). Its variance is
Thus \(T\sim N(11,17)\). Standardizing the threshold gives
where \(Z\sim N(0,1)\). The probability is rounded to four decimal places. Check: the variance is positive, and because \(15\) is about \(0.97\) standard deviations above the mean, a right-tail probability near \(0.17\) is reasonable. The common mistake is to omit the covariance term, which would give variance \(16+9=25\) and the incorrect probability \(P(Z>0.8)\approx0.2119\). That calculation wrongly treats the variables as independent.
Several Variables and the Covariance Matrix
For more than two variables, the same principle applies to every pair. Let the vector of jointly normal variables be \(\mathbf{X}=(X_1,\ldots,X_n)^\mathsf{T}\), with mean vector \(\boldsymbol{\mu}\) and covariance matrix \(\boldsymbol{\Sigma}\). The diagonal entries of \(\boldsymbol{\Sigma}\) are the variances; each off-diagonal entry is the covariance of the corresponding pair. For fixed coefficient vector \(\mathbf{a}\), the combination \(L=\mathbf{a}^\mathsf{T}\mathbf{X}\) has mean \(\mathbf{a}^\mathsf{T}\boldsymbol{\mu}\) and variance \(\mathbf{a}^\mathsf{T}\boldsymbol{\Sigma}\mathbf{a}\).
The matrix expression is a compact way to include each variance and covariance contribution exactly once. In the expanded sum, the factor of \(2\) appears because each covariance would arise twice when the quadratic expression is expanded.
Worked Example: A Three-Variable Normal Combination
Suppose \((X,Y,Z)\) is jointly normal with means \(E[X]=2\), \(E[Y]=-1\), and \(E[Z]=4\), and covariance matrix
Find the distribution of \(L=3X+Y-Z\), then calculate \(P(L>1+2\sqrt{30})\).
Method. Use the coefficients \((3,1,-1)\) in the mean and variance formulas. The covariance matrix is a valid positive-definite matrix: its one-dimensional principal minors are positive, its two-dimensional principal minors are \(3\), \(8\), and \(32\), and its determinant is \(15\).
The mean is
The diagonal variance contributions are \(3^2(1)+1^2(4)+(-1)^2(9)=9+4+9=22\). The three covariance contributions are
Therefore, \(\operatorname{Var}(L)=22+6+6-4=30\), and \(L\sim N(1,30)\). The event threshold is two standard deviations above the mean:
The probability is rounded to four decimal places. Check: the threshold is exactly \(2\) standard deviations above the mean, so the upper-tail probability is the standard normal upper tail at \(2\). A common mistake is to discard the off-diagonal entries and use variance \(22\); this is incorrect because the variables are not independent. The covariance contributions here sum to \(8\), giving the correct variance \(30\).
When the Normal Result Applies—and When It Does Not
There are two separate questions to answer: what are the mean and variance of the combination, and is its distribution normal? The mean and variance formulas apply to variables with finite second moments whether or not they are normal. An exact normal distribution for the combination is guaranteed when the variables are jointly normal.
Normal marginals alone do not guarantee joint normality. Dependence can be arranged so that each variable separately has a normal distribution while a linear combination does not. Therefore, do not conclude that a sum or difference is normal merely because each component is described as normal; look for a joint normal model or an explicit statement that the variables are jointly normal.
For jointly normal variables, zero covariance between a pair implies that pair is independent. This equivalence is special to the joint normal model; in general, zero covariance does not establish independence. Conversely, independence of two variables with normal distributions implies they are jointly normal, so the same calculation applies.
For \(aX+bY\), record \(a\) and \(b\), including any negative sign. For a longer combination, list every coefficient in order.
Multiply each component mean by its coefficient and add the results.
Square each coefficient for its variance term, then include every pairwise covariance term with its signed coefficient product.
If the variables are jointly normal, state the combination’s normal distribution. Otherwise, do not claim exact normality from the mean and variance alone.
Subtract the combination’s mean and divide by its standard deviation before using the standard normal distribution.
Common Exam Traps
- Assuming independence without evidence. Independence sets the covariance to zero; it is not implied by having separate normal distributions.
- Dropping the coefficient signs. In a difference, a negative coefficient can reverse the sign of a covariance contribution.
- Using correlation as covariance. Convert using \(\operatorname{Cov}(X,Y)=\rho_{X,Y}\sigma_X\sigma_Y\) before inserting the value into a variance formula.
- Using variance where standard deviation is required. Standardization divides by the square root of the variance, not by the variance itself.
- Claiming an exact normal result too soon. The variance calculation is general, but the exact normal distribution requires joint normality or another sufficient model assumption.
Check Your Understanding
For each problem, identify the combination’s distribution where possible and use the covariance terms rather than assuming independence.
- \(U\) and \(V\) are jointly normal with means \(4\) and \(7\), variances \(1\) and \(4\), and correlation \(0.5\). Find the distribution of \(U+V\) and \(P(U+V>11+\sqrt{7})\). (Answer: \(U+V\sim N(11,7)\); probability \(0.1587\).)
- \(A\) and \(B\) are jointly normal with means \(20\) and \(18\), variances \(9\) and \(4\), and correlation \(-0.5\). Find the distribution of \(A-2B\) and \(P(A-2B<-16-\sqrt{37})\). (Answer: \(A-2B\sim N(-16,37)\); probability \(0.1587\).)
- \((A,B,C)\) are jointly normal with mean vector \((1,0,2)^\mathsf{T}\) and covariance matrix \(\begin{pmatrix}4&1&0\\1&1&0\\0&0&9\end{pmatrix}\). Find the distribution of \(A+2B-C\). (Answer: \(A+2B-C\sim N(-1,21)\).)
- \(X\) and \(Y\) are jointly normal with means \(12\) and \(8\), variances \(9\) and \(16\), and \(\operatorname{Var}(X-Y)=17\). Find their correlation and \(P(X-Y>7)\). (Answer: \(\operatorname{Cov}(X,Y)=4\), correlation \(1/3\), and probability \(0.2335\).)
- \(R\) and \(S\) are jointly normal with means \(50\) and \(40\), variances \(25\) and \(36\), and correlation \(0.4\). Find the distribution of \(2R+S\) and \(P(2R+S>140+2\sqrt{184})\). (Answer: \(2R+S\sim N(140,184)\); probability \(0.0228\).)