Sum of normally distributed random variables

In probability theory, calculation of the sum of normally distributed random variables is an instance of the arithmetic of random variables.

This is not to be confused with the sum of normal distributions which forms a mixture distribution.

Independent random variables

edit

Let X and Y be independent random variables that are normally distributed (and therefore also jointly so), then their sum is also normally distributed. i.e., if

 
 
 

then

 

This means that the sum of two independent normally distributed random variables is normal, with its mean being the sum of the two means, and its variance being the sum of the two variances (i.e., the square of the standard deviation is the sum of the squares of the standard deviations).[1]

In order for this result to hold, the assumption that X and Y are independent cannot be dropped, although it can be weakened to the assumption that X and Y are jointly, rather than separately, normally distributed.[2] (See here for an example.)

The result about the mean holds in all cases, while the result for the variance requires uncorrelatedness, but not independence.

Proofs

edit

Proof using characteristic functions

edit

The characteristic function

 

of the sum of two independent random variables X and Y is just the product of the two separate characteristic functions:

 

of X and Y.

The characteristic function of the normal distribution with expected value μ and variance σ2 is

 

So

 

This is the characteristic function of the normal distribution with expected value   and variance  

Finally, recall that no two distinct distributions can both have the same characteristic function, so the distribution of X + Y must be just this normal distribution.

Proof using convolutions

edit

For independent random variables X and Y, the distribution fZ of Z = X + Y equals the convolution of fX and fY:

 

Given that fX and fY are normal densities,

 

Substituting into the convolution:

 

Defining  , and completing the square:

 

The expression in the integral is a normal density distribution on x, and so the integral evaluates to 1. The desired result follows:

 

It can be shown that the Fourier transform of a Gaussian,  , is[3]

 

By the convolution theorem:

 

Geometric proof

edit

First consider the normalized case when X, Y ~ N(0, 1), so that their PDFs are

 

and

 

Let Z = X + Y. Then the CDF for Z will be

 

This integral is over the half-plane which lies under the line x+y = z.

The key observation is that the function

 

is radially symmetric. So we rotate the coordinate plane about the origin, choosing new coordinates   such that the line x+y = z is described by the equation   where   is determined geometrically. Because of the radial symmetry, we have  , and the CDF for Z is

 

This is easy to integrate; we find that the CDF for Z is

 

To determine the value  , note that we rotated the plane so that the line x+y = z now runs vertically with x-intercept equal to c. So c is just the distance from the origin to the line x+y = z along the perpendicular bisector, which meets the line at its nearest point to the origin, in this case  . So the distance is  , and the CDF for Z is  , i.e.,  

Now, if a, b are any real constants (not both zero) then the probability that   is found by the same integral as above, but with the bounding line  . The same rotation method works, and in this more general case we find that the closest point on the line to the origin is located a (signed) distance

 

away, so that

 

The same argument in higher dimensions shows that if

 

then

 

Now we are essentially done, because

 

So in general, if

 

then

 

Correlated random variables

edit

In the event that the variables X and Y are jointly normally distributed random variables, then X + Y is still normally distributed (see Multivariate normal distribution) and the mean is the sum of the means. However, the variances are not additive due to the correlation. Indeed,

 

where ρ is the correlation. In particular, whenever ρ < 0, then the variance is less than the sum of the variances of X and Y.

Extensions of this result can be made for more than two random variables, using the covariance matrix.

Proof

edit

In this case (with X and Y having zero means), one needs to consider

 

As above, one makes the substitution  

This integral is more complicated to simplify analytically, but can be done easily using a symbolic mathematics program. The probability distribution fZ(z) is given in this case by

 

where

 

If one considers instead Z = X − Y, then one obtains

 

which also can be rewritten with

 

The standard deviations of each distribution are obvious by comparison with the standard normal distribution.

References

edit
  1. ^ Lemons, Don S. (2002), An Introduction to Stochastic Processes in Physics, The Johns Hopkins University Press, p. 34, ISBN 0-8018-6866-1
  2. ^ Lemons (2002) pp. 35–36
  3. ^ Derpanis, Konstantinos G. (October 20, 2005). "Fourier Transform of the Gaussian" (PDF).

See also

edit