Normally distributed and uncorrelated does not imply independent

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

In probability theory, two random variables being uncorrelated does not imply their independence. In some contexts, uncorrelatedness implies at least pairwise independence (as when the random variables involved have Bernoulli distributions).

It is sometimes mistakenly thought that one context in which uncorrelatedness implies independence is when the random variables involved are normally distributed. However, this is incorrect if the variables are merely marginally normally distributed but not jointly normally distributed.

Suppose two random variables X and Y are jointly normally distributed. That is the same as saying that the random vector (XY) has a multivariate normal distribution. It means that the joint probability distribution of X and Y is such that each linear combination of X and Y is normally distributed, i.e. for any two constant (i.e., non-random) scalars a and b, the random variable aX + bY is normally distributed. In that case if X and Y are uncorrelated, i.e., their covariance cov(XY) is zero, then they are independent.[1] However, it is possible for two random variables X and Y to be so distributed jointly that each one alone is marginally normally distributed, and they are uncorrelated, but they are not independent; examples are given below.

Examples

A symmetric example

File:Uncorrelated sym.png
Joint range of X and Y. Darker indicates higher value of the density function.

Suppose X has a normal distribution with expected value 0 and variance 1. Let W have the Rademacher distribution, so that W = 1 or −1, each with probability 1/2, and assume W is independent of X. Let Y = WX. Then[2]

  • X and Y are uncorrelated;
  • Both have the same normal distribution; and
  • X and Y are not independent.

Note that the distribution of the simple linear combination X + Y concentrates positive probability at 0: Pr(X + Y = 0)  = 1/2 and so is not normally distributed. By the definition above, X and Y are not jointly normally distributed.

To see that X and Y are uncorrelated, consider

 \begin{align}
\operatorname{cov}(X,Y) &{} = E(XY) - E(X)E(Y) = E(XY) - 0 = E(E(XY\mid W)) \\
& {} = E(X^2)\Pr(W=1) + E(-X^2)\Pr(W=-1) \\
& {} = 1\cdot\frac12 + (-1)\cdot\frac12 = 0.
\end{align}

To see that Y has the same normal distribution as X, consider

\begin{align}
\Pr(Y \le x) & {} = E(\Pr(Y \le x\mid W)) \\
& {} = \Pr(X \le x)\Pr(W = 1) + \Pr(-X\le x)\Pr(W = -1) \\
& {} = \Phi(x) \cdot\frac12 + \Phi(x)\cdot\frac12
\end{align}

(since X and −X both have the same normal distribution), where \Phi(x) is the cumulative distribution function of the normal distribution..

To see that X and Y are not independent, observe that |Y| = |X| or that Pr(Y > 1 | X = 1/2) = Pr(X > 1 | X = 1/2) = 0.

An asymmetric example

File:Uncorrelated asym.png
The joint density of X and Y. Darker indicates a higher value of the density.

Suppose X has a normal distribution with expected value 0 and variance 1. Let

Y=\left\{\begin{matrix} X & \text{if } \left|X\right| \leq c \\
-X & \text{if } \left|X\right|>c \end{matrix}\right.

where c is a positive number to be specified below. If c is very small, then the correlation corr(XY) is near −1; if c is very large, then corr(XY) is near 1. Since the correlation is a continuous function of c, the intermediate value theorem implies there is some particular value of c that makes the correlation 0. That value is approximately 1.54. In that case, X and Y are uncorrelated, but they are clearly not independent, since X completely determines Y.

To see that Y is normally distributed—indeed, that its distribution is the same as that of X—let us find its cumulative distribution function:

\begin{align}\Pr(Y \leq x) &= \Pr(\{|X| \leq c\text{ and }X \leq x\}\text{ or }\{|X|>c\text{ and }-X \leq x\})\\
&= \Pr(|X| \leq c\text{ and }X \leq x) + \Pr(|X|>c\text{ and }-X \leq x)\\
&= \Pr(|X| \leq c\text{ and }X \leq x) + \Pr(|X|>c\text{ and }X \leq x) \\
&= \Pr(X \leq x). \end{align}\,

where the next-to-last equality follows from the symmetry of the distribution of X and the symmetry of the condition that |X| ≤ c.

Observe that the difference X − Y is nowhere near being normally distributed, since it has a substantial probability (about 0.88) of it being equal to 0, whereas the normal distribution, being a continuous distribution, has no discrete part, i.e., does not concentrate more than zero probability at any single point. Consequently X and Y are not jointly normally distributed, even though they are separately normally distributed.[3]

References

  1. Lua error in package.lua at line 80: module 'strict' not found.
  2. UIUC, Lecture 21. The Multivariate Normal Distribution, 21.6:"Individually Gaussian Versus Jointly Gaussian".
  3. Edward L. Melnick and Aaron Tenenbein, "Misspecifications of the Normal Distribution", The American Statistician, volume 36, number 4 November 1982, pages 372–373