Tutorials › VEE Mathematical Statistics › Covariance in Bivariate Normal Models

VEE Mathematical Statistics · Covariance and Correlation · Tutorial 324 of 880

Covariance in Bivariate Normal Models

Connect the bivariate normal density’s parameters to covariance, correlation, and independence, then apply those relationships in worked calculations.

Intermediate 11 min read

What You'll Learn

  • Identify the means, standard deviations, covariance, and correlation in a bivariate normal model
  • Convert between the covariance parameter and the correlation parameter
  • Check whether proposed bivariate normal parameters define a nonsingular density
  • Interpret the cross-product term in the bivariate normal density
  • Explain why zero correlation implies independence for a bivariate normal pair
  • Calculate joint densities and probabilities using the model’s dependence structure

Why Covariance Matters in a Bivariate Normal Model

Two variables can each have a familiar normal distribution while their joint behavior still depends on how they move together. For a bivariate normal model, that dependence is summarized by one parameter: either the covariance or, equivalently, the correlation. Understanding how that parameter enters the density lets you interpret the model, check whether its parameters are possible, and decide when zero correlation does imply independence.

Throughout, consider two continuous random variables \(X\) and \(Y\) with means \(\mu_X,\mu_Y\), positive standard deviations \(\sigma_X,\sigma_Y\), and correlation \(\rho\). The ordinary density formula below applies when \(-1<\rho<1\). The endpoints \(\rho=-1\) and \(\rho=1\) describe singular cases, not a two-dimensional density of this form.

Definition: A nondegenerate bivariate normal pair \((X,Y)\) has joint density $$ f_{X,Y}(x,y)= \frac{1}{2\pi\sigma_X\sigma_Y\sqrt{1-\rho^2}} \exp\left\{-\frac{1}{2(1-\rho^2)} \left[ \left(\frac{x-\mu_X}{\sigma_X}\right)^2 -2\rho\left(\frac{x-\mu_X}{\sigma_X}\right)\left(\frac{y-\mu_Y}{\sigma_Y}\right) +\left(\frac{y-\mu_Y}{\sigma_Y}\right)^2 \right]\right\}, $$ where \(\sigma_X>0\), \(\sigma_Y>0\), and \(-1<\rho<1\). The marginal distributions are \(X\sim N(\mu_X,\sigma_X^2)\) and \(Y\sim N(\mu_Y,\sigma_Y^2)\), and $$ \operatorname{Cov}(X,Y)=\rho\sigma_X\sigma_Y. $$

The covariance is sometimes used directly as the dependence parameter. Writing \(c=\operatorname{Cov}(X,Y)\), the same relationship becomes \(\rho=c/(\sigma_X\sigma_Y)\). Thus the two parameters measure dependence on different scales: covariance has units equal to the product of the variables’ units, whereas correlation is dimensionless and lies between \(-1\) and \(1\).

The density’s cross-product term, \(-2\rho z_Xz_Y\), where \(z_X=(x-\mu_X)/\sigma_X\) and \(z_Y=(y-\mu_Y)/\sigma_Y\), is the part that reflects dependence. If \(\rho\) is positive, points for which both standardized values have the same sign are favored relative to points with opposite signs. If \(\rho\) is negative, the opposite-sign combinations are favored. The overall density remains a joint model; the cross-product term should not be interpreted in isolation from the other terms and the normalizing constant.

Why the Covariance Parameter Is \(\rho\sigma_X\sigma_Y\)

A useful way to verify the covariance relationship is to express the pair using two independent standard normal variables \(Z_1\) and \(Z_2\). For \(-1<\rho<1\), write

$$ X=\mu_X+\sigma_XZ_1,\qquad Y=\mu_Y+\sigma_Y\left(\rho Z_1+\sqrt{1-\rho^2}\,Z_2\right). $$

Both standardized components have mean zero. Since \(Z_1\) and \(Z_2\) are independent with variance one, the variance of the standardized component of \(Y\) is \(\rho^2+(1-\rho^2)=1\). The covariance follows by expanding the product of the centered variables:

$$ \begin{aligned} \operatorname{Cov}(X,Y) &=\sigma_X\sigma_Y \operatorname{Cov}\left(Z_1,\rho Z_1+\sqrt{1-\rho^2}\,Z_2\right)\\ &=\sigma_X\sigma_Y\left(\rho\operatorname{Var}(Z_1) +\sqrt{1-\rho^2}\operatorname{Cov}(Z_1,Z_2)\right)\\ &=\rho\sigma_X\sigma_Y. \end{aligned} $$

Here \(\operatorname{Cov}(Z_1,Z_2)=0\) because the standard normal variables are independent, and \(\operatorname{Var}(Z_1)=1\). This construction also explains the density: transforming \((Z_1,Z_2)\) to \((X,Y)\) gives the stated bivariate normal density by the change-of-variables formula. In particular, the absolute Jacobian determinant is \(\sigma_X\sigma_Y\sqrt{1-\rho^2}\), which accounts for the density’s normalizing factor.

Worked Example: Converting Covariance to Correlation

A hypothetical model describes two quarterly claim-cost deviations, \(X\) and \(Y\), measured in thousands of dollars. Suppose their means are both zero, their standard deviations are \(8\) and \(15\) thousand dollars, and their covariance is \(-60\) thousand-dollars squared. Find the correlation, check whether a nondegenerate bivariate normal density is possible, and find the density at the pair of means.

Method. Standardize the covariance by the product of the standard deviations. A valid nondegenerate model requires the resulting correlation to be strictly between \(-1\) and \(1\).

First calculate the product of the standard deviations:

$$ \sigma_X\sigma_Y=8(15)=120. $$

Then calculate the correlation:

$$ \rho=\frac{\operatorname{Cov}(X,Y)}{\sigma_X\sigma_Y} =\frac{-60}{120}=-0.5. $$

Since \(-1<-0.5<1\), the parameters define a nondegenerate bivariate normal density. At \(x=\mu_X\) and \(y=\mu_Y\), both standardized values are zero, so the exponential factor is \(1\). Therefore,

$$ f_{X,Y}(\mu_X,\mu_Y) =\frac{1}{2\pi(8)(15)\sqrt{1-(-0.5)^2}} =\frac{1}{240\pi\sqrt{0.75}} \approx 0.001530. $$

The density is approximately \(0.001530\) per squared thousand dollars. The negative correlation indicates that deviations in opposite directions are associated in this model. The common mistake this calculation avoids is treating covariance as correlation without dividing by both standard deviations.

Reading the Density and Checking the Parameters

The covariance and correlation cannot be chosen independently of the marginal standard deviations. Since \(\rho=c/(\sigma_X\sigma_Y)\), a proposed covariance \(c\) must satisfy \(|c|<\sigma_X\sigma_Y\) for the nondegenerate density above. Equivalently, the covariance matrix

$$ \begin{pmatrix} \sigma_X^2 & c\\ c & \sigma_Y^2 \end{pmatrix} $$

must have positive determinant: \(\sigma_X^2\sigma_Y^2-c^2>0\). If equality holds, then \(|\rho|=1\), and the joint distribution lies on a straight line rather than having the displayed two-dimensional density. If \(|c|>\sigma_X\sigma_Y\), the proposed parameters are invalid.

When using the density, first calculate the standardized values \(z_X\) and \(z_Y\), then substitute them into the entire quadratic expression. In particular, preserve the signs of both standardized values in the cross-product. A negative \(z_Y\) changes the sign of \(z_Xz_Y\), which can change the exponent substantially.

Worked Example: Evaluating the Density Away from the Means

Suppose a hypothetical bivariate normal model for two continuous operational measurements has \(\mu_X=100\), \(\mu_Y=50\), \(\sigma_X=4\), \(\sigma_Y=10\), and \(\rho=0.6\). Find the covariance and evaluate the joint density at \(x=104\), \(y=40\).

Method. Use \(c=\rho\sigma_X\sigma_Y\) for the covariance, then standardize the observed values and substitute them into the bivariate normal density.

The covariance is

$$ c=0.6(4)(10)=24. $$

The standardized values are

$$ z_X=\frac{104-100}{4}=1,\qquad z_Y=\frac{40-50}{10}=-1. $$

The expression in square brackets in the exponent is

$$ z_X^2-2\rho z_Xz_Y+z_Y^2 =1-2(0.6)(1)(-1)+1 =1+1.2+1=3.2. $$

The exponent is therefore \(-3.2/[2(1-0.6^2)]=-3.2/1.28=-2.5\). The normalizing factor is

$$ \frac{1}{2\pi(4)(10)\sqrt{1-0.6^2}} =\frac{1}{64\pi} \approx 0.0049736. $$

Thus,

$$ f_{X,Y}(104,40) =\frac{1}{64\pi}e^{-2.5} \approx 0.0004083. $$

The covariance is \(24\) in the product of the measurements’ units, and the density is approximately \(0.0004083\) in reciprocal product-units. The key sign check is that \(z_Xz_Y=-1\), so the cross-product contribution is positive inside the square brackets. Using \(+1\) for that product would give the wrong exponent.

Zero Correlation and Independence

Zero covariance does not generally establish independence. In the nondegenerate bivariate normal model, however, the zero-correlation case has a special property: the joint density factors into the product of the marginal densities. This factorization is the direct reason that \(X\) and \(Y\) are independent when \(\rho=0\).

Set \(\rho=0\) in the joint density. The cross-product term disappears, the denominator contains \(\sqrt{1-\rho^2}=1\), and the exponent separates into a term involving only \(x\) and a term involving only \(y\):

$$ \begin{aligned} f_{X,Y}(x,y) &=\frac{1}{2\pi\sigma_X\sigma_Y} \exp\left\{-\frac{1}{2} \left[ \left(\frac{x-\mu_X}{\sigma_X}\right)^2+ \left(\frac{y-\mu_Y}{\sigma_Y}\right)^2 \right]\right\}\\ &=\left[ \frac{1}{\sqrt{2\pi}\sigma_X} \exp\left\{-\frac{(x-\mu_X)^2}{2\sigma_X^2}\right\} \right] \left[ \frac{1}{\sqrt{2\pi}\sigma_Y} \exp\left\{-\frac{(y-\mu_Y)^2}{2\sigma_Y^2}\right\} \right]\\ &=f_X(x)f_Y(y). \end{aligned} $$

A joint density that equals the product of its marginal densities gives independence: probabilities over rectangular regions factor into the corresponding marginal probabilities. Conversely, if \(X\) and \(Y\) are independent and have finite variances, then their covariance is zero. Because the standard deviations in this model are positive, \(\operatorname{Cov}(X,Y)=\rho\sigma_X\sigma_Y=0\) implies \(\rho=0\). Thus, for a nondegenerate bivariate normal pair, independence and zero correlation are equivalent.

Worked Example: A Joint Probability When Correlation Is Zero

Suppose a hypothetical bivariate normal model describes two continuous measurements \(X\) and \(Y\), with means \(\mu_X=120\) and \(\mu_Y=90\), standard deviations \(12\) and \(20\), and correlation \(\rho=0\). Find \(P(X>120,\;Y>90)\).

Method. First use the bivariate normal factorization at zero correlation to establish independence. Then use the symmetry of each normal marginal about its mean.

Because \(\rho=0\), the joint density factors, so \(X\) and \(Y\) are independent. A normal distribution is symmetric about its mean, and it assigns probability zero to any single point. Therefore,

$$ P(X>120)=0.5,\qquad P(Y>90)=0.5. $$

Independence gives

$$ P(X>120,\;Y>90) =P(X>120)P(Y>90) =(0.5)(0.5)=0.25. $$

The answer is \(0.25\). The check is that each event has probability one-half and, under independence, their intersection has probability one-quarter. The zero-correlation conclusion is justified here by the bivariate normal model, not by a general rule for arbitrary pairs of random variables.

Common Exam Traps and Practical Checks

  • Confusing covariance with correlation. The correct conversion is \(\rho=c/(\sigma_X\sigma_Y)\). For example, if the standard deviations are \(4\) and \(10\) and the covariance is \(24\), then \(\rho=24/(4\cdot10)=0.6\), not \(24\).
  • Ignoring the allowed parameter range. A nondegenerate density requires \(|c|<\sigma_X\sigma_Y\). A proposed covariance of \(45\) with standard deviations \(5\) and \(8\) is impossible for this model because \(45>5(8)=40\).
  • Assuming zero correlation always means independence. That implication is valid for the bivariate normal model, but it is not valid for every joint distribution. Use the model assumption before making the independence conclusion.
  • Losing a sign in the density. Compute both standardized values before evaluating the cross-product. If one is negative and the other positive, their product is negative.
  • Treating the endpoints as ordinary density cases. At \(\rho=\pm1\), the covariance matrix is singular, so the given two-dimensional density formula does not apply.
Key takeaway: In a nondegenerate bivariate normal model, \(\operatorname{Cov}(X,Y)=\rho\sigma_X\sigma_Y\). The covariance must be consistent with the marginal standard deviations, and zero correlation makes the density factor into its marginals, proving independence for this model.

Check Your Understanding

For each problem, identify the relevant model relationship before calculating.

  1. A bivariate normal pair has standard deviations \(6\) and \(9\), with covariance \(27\). Find its correlation and state whether the parameters allow a nondegenerate density. (Answer: \(\rho=0.5\); yes.)
  2. A proposed bivariate normal pair has standard deviations \(3\) and \(7\), with covariance \(-14\). Find the correlation and state whether the density is nondegenerate. (Answer: \(\rho=-\frac{2}{3}\approx-0.6667\); yes.)
  3. A bivariate normal pair has standard deviations \(5\) and \(8\) and correlation \(-0.4\). Calculate its covariance, including its sign. (Answer: \(-16\).)
  4. A bivariate normal pair has correlation zero. Find \(P(X\leq\mu_X,\;Y>\mu_Y)\), explaining which model property permits multiplication. (Answer: \(0.25\), by independence and symmetry of the normal marginals.)
  5. A bivariate normal pair has standard deviations \(5\) and \(12\), and covariance \(45\). Find its correlation and determine whether the variables are independent. (Answer: \(\rho=0.75\); no, since the nonzero correlation implies dependence in this model.)
Clinical Biostats