In probability theory and statistics, **variance** is the expectation of the squared deviation of a random variable from its mean, and it informally measures how far a set of (random) numbers are spread out from their mean. The variance has a central role in statistics. It is used in descriptive statistics, statistical inference, hypothesis testing, goodness of fit, and Monte Carlo sampling, amongst many others. This makes it a central quantity in numerous fields such as physics, biology, chemistry, cryptography, economics, and finance. The variance is the square of the standard deviation, the second central moment of a distribution, and the covariance of the random variable with itself, and it is often represented by

## Contents

- Definition
- Continuous random variable
- Discrete random variable
- Normal distribution
- Exponential distribution
- Poisson distribution
- Binomial distribution
- Coin toss
- Fair die
- Basic properties
- Sum of uncorrelated variables Bienaym formula
- Sum of correlated variables
- Matrix notation for the variance of a linear combination
- Weighted sum of variables
- Product of independent variables
- Product of statistically dependent variables
- Decomposition
- Formulae for the variance
- Calculation from the CDF
- Characteristic property
- Units of measurement
- Approximating the variance of a function
- Population variance
- Sample variance
- Distribution of the sample variance
- Samuelsons inequality
- Relations with the harmonic and arithmetic means
- Tests of equality of variances
- History
- Moment of inertia
- Semivariance
- Generalizations
- References

## Definition

The variance of a random variable

This definition encompasses random variables that are generated by processes that are discrete, continuous, neither, or mixed. The variance can also be thought of as the covariance of a random variable with itself:

The variance is also equivalent to the second cumulant of a probability distribution that generates

A mnemonic for the above expression is "mean of square minus square of mean". On computational floating point arithmetic, this equation should not be used, because it suffers from catastrophic cancellation if the two components of the equation are similar in magnitude. There exist numerically stable alternatives.

## Continuous random variable

If the random variable

where

and where the integrals are definite integrals taken for

If a continuous distribution does not have an expected value, as is the case for the Cauchy distribution, it does not have a variance either. Many other distributions for which the expected value does exist also do not have a finite variance because the integral in the variance definition diverges. An example is a Pareto distribution whose index

## Discrete random variable

If the generator of random variable

or equivalently

where

(When such a discrete weighted variance is specified by weights whose sum is not 1, then one divides by the sum of the weights.)

The variance of a set of

where

The variance of a set of

## Normal distribution

The normal distribution with parameters

In this distribution,

The role of the normal distribution in the central limit theorem is in part responsible for the prevalence of the variance in probability and statistics.

## Exponential distribution

The exponential distribution with parameter

and it has expected value

So for an exponentially distributed random variable,

## Poisson distribution

The Poisson distribution with parameter

and it has expected value

So for a Poisson-distributed random variable,

## Binomial distribution

The binomial distribution with parameters

and it has expected value

## Coin toss

The binomial distribution with

## Fair die

A six-sided fair die can be modelled with a discrete random variable with outcomes 1 through 6, each with equal probability

The general formula for the variance of the outcome *X* of a die of n sides is

## Basic properties

Variance is non-negative because the squares are positive or zero.

The variance of a constant random variable is zero, and if the variance of a variable in a data set is 0, then all the entries have the same value.

Variance is invariant with respect to changes in a location parameter. That is, if a constant is added to all values of the variable, the variance is unchanged.

If all values are scaled by a constant, the variance is scaled by the square of that constant.

The variance of a sum of two random variables is given by:

where Cov(⋅, ⋅) is the covariance. In general we have for the sum of

These results lead to the variance of a linear combination as:

If the random variables

they are said to be uncorrelated. It follows immediately from the expression given earlier that if the random variables

Since independent random variables are always uncorrelated, the equation above holds in particular when the random variables

## Sum of uncorrelated variables (Bienaymé formula)

One reason for the use of the variance in preference to other measures of dispersion is that the variance of the sum (or the difference) of uncorrelated random variables is the sum of their variances:

This statement is called the Bienaymé formula and was discovered in 1853. It is often made with the stronger condition that the variables are independent, but being uncorrelated suffices. So if all the variables have the same variance σ^{2}, then, since division by *n* is a linear transformation, this formula immediately implies that the variance of their mean is

That is, the variance of the mean decreases when *n* increases. This formula for the variance of the mean is used in the definition of the standard error of the sample mean, which is used in the central limit theorem.

## Sum of correlated variables

In general, if the variables are correlated, then the variance of their sum is the sum of their covariances:

(Note: The second equality comes from the fact that Cov(*X*_{i},*X*_{i}) = Var(*X*_{i}).)

Here Cov(⋅, ⋅) is the covariance, which is zero for independent random variables (if it exists). The formula states that the variance of a sum is equal to the sum of all elements in the covariance matrix of the components. The next expression states equivalently that the variance of the sum is the sum of the diagonal of covariance matrix plus two times the sum of its upper triangular elements (or its lower triangular elements); this emphasizes that the covariance matrix is symmetric. This formula is used in the theory of Cronbach's alpha in classical test theory.

So if the variables have equal variance *σ*^{2} and the average correlation of distinct variables is *ρ*, then the variance of their mean is

This implies that the variance of the mean increases with the average of the correlations. In other words, additional correlated observations are not as effective as additional independent observations at reducing the uncertainty of the mean. Moreover, if the variables have unit variance, for example if they are standardized, then this simplifies to

This formula is used in the Spearman–Brown prediction formula of classical test theory. This converges to *ρ* if *n* goes to infinity, provided that the average correlation remains constant or converges too. So for the variance of the mean of standardized variables with equal correlations or converging average correlation we have

Therefore, the variance of the mean of a large number of standardized variables is approximately equal to their average correlation. This makes clear that the sample mean of correlated variables does not generally converge to the population mean, even though the Law of large numbers states that the sample mean will converge for independent variables.

## Matrix notation for the variance of a linear combination

Define

## Weighted sum of variables

The scaling property and the Bienaymé formula, along with the property of the covariance Cov(*aX*, *bY*) = *ab* Cov(*X*, *Y*) jointly imply that

This implies that in a weighted sum of variables, the variable with the largest weight will have a disproportionally large weight in the variance of the total. For example, if *X* and *Y* are uncorrelated and the weight of *X* is two times the weight of *Y*, then the weight of the variance of *X* will be four times the weight of the variance of *Y*.

The expression above can be extended to a weighted sum of multiple variables:

## Product of independent variables

If two variables X and Y are independent, the variance of their product is given by

Equivalently, using the basic properties of expectation, it is given by

## Product of statistically dependent variables

In general, if two variables are statistically dependent, the variance of their product is given by:

## Decomposition

The general formula for variance decomposition or the law of total variance is: If

where

In particular, if

where

where

A similar formula is applied in analysis of variance, where the corresponding formula is

here

This can also be derived from the additivity of variances, since the total (observed) score is the sum of the predicted score and the error score, where the latter two are uncorrelated.

Similar decompositions are possible for the sum of squared deviations (sum of squares,

## Formulae for the variance

A formula often used for deriving the variance of a theoretical distribution is as follows:

This will be useful when it is possible to derive formulae for the expected value and for the expected value of the square.

This formula is also sometimes used in connection with the sample variance. While useful for hand calculations, it is not advised for computer calculations as it suffers from catastrophic cancellation if the two components of the equation are similar in magnitude and floating point arithmetic is used. This is discussed in the article Algorithms for calculating variance.

## Calculation from the CDF

The population variance for a non-negative random variable can be expressed in terms of the cumulative distribution function *F* using

This expression can be used to calculate the variance in situations where the CDF, but not the density, can be conveniently expressed.

## Characteristic property

The second moment of a random variable attains the minimum value when taken around the first moment (i.e., mean) of the random variable, i.e. *X*, then it is necessarily of the form *a* > 0. This also holds in the multidimensional case.

## Units of measurement

Unlike expected absolute deviation, the variance of a variable has units that are the square of the units of the variable itself. For example, a variable measured in meters will have a variance measured in meters squared. For this reason, describing data sets via their standard deviation or root mean square deviation is often preferred over using the variance. In the dice example the standard deviation is √2.9 ≈ 1.7, slightly larger than the expected absolute deviation of 1.5.

The standard deviation and the expected absolute deviation can both be used as an indicator of the "spread" of a distribution. The standard deviation is more amenable to algebraic manipulation than the expected absolute deviation, and, together with variance and its generalization covariance, is used frequently in theoretical statistics; however the expected absolute deviation tends to be more robust as it is less sensitive to outliers arising from measurement anomalies or an unduly heavy-tailed distribution.

## Approximating the variance of a function

The delta method uses second-order Taylor expansions to approximate the variance of a function of one or more random variables: see Taylor expansions for the moments of functions of random variables. For example, the approximate variance of a function of one variable is given by

provided that *f* is twice differentiable and that the mean and variance of *X* are finite.

Real-world observations such as the measurements of yesterday's rain throughout the day typically cannot be complete sets of all possible observations that could be made. As such, the variance calculated from the finite set will in general not match the variance that would have been calculated from the full population of possible observations. This means that one estimates the mean and variance that would have been calculated from an omniscient set of observations by using an estimator equation. The estimator is a function of the sample of *n* observations drawn without observational bias from the whole population of potential observations. In this example that sample would be the set of actual measurements of yesterday's rainfall from available rain gauges within the geography of interest.

The simplest estimators for population mean and population variance are simply the mean and variance of the sample, the **sample mean** and **(uncorrected) sample variance** – these are consistent estimators (they converge to the correct value as the number of samples increases), but can be improved. Estimating the population variance by taking the sample's variance is close to optimal in general, but can be improved in two ways. Most simply, the sample variance is computed as an average of squared deviations about the (sample) mean, by dividing by *n.* However, using values other than *n* improves the estimator in various ways. Four common values for the denominator are *n,* *n* − 1, *n* + 1, and *n* − 1.5: *n* is the simplest (population variance of the sample), *n* − 1 eliminates bias, *n* + 1 minimizes mean squared error for the normal distribution, and *n* − 1.5 mostly eliminates bias in unbiased estimation of standard deviation for the normal distribution.

Firstly, if the omniscient mean is unknown (and is computed as the sample mean), then the sample variance is a biased estimator: it underestimates the variance by a factor of (*n* − 1) / *n*; correcting by this factor (dividing by *n* − 1 instead of *n*) is called Bessel's correction. The resulting estimator is unbiased, and is called the **(corrected) sample variance** or **unbiased sample variance**. For example, when *n* = 1 the variance of a single observation about the sample mean (itself) is obviously zero regardless of the population variance. If the mean is determined in some other way than from the same samples used to estimate the variance then this bias does not arise and the variance can safely be estimated as that of the samples about the (independently known) mean.

Secondly, the sample variance does not generally minimize mean squared error between sample variance and population variance. Correcting for bias often makes this worse: one can always choose a scale factor that performs better than the corrected sample variance, though the optimal scale factor depends on the excess kurtosis of the population (see mean squared error: variance), and introduces bias. This always consists of scaling down the unbiased estimator (dividing by a number larger than *n* − 1), and is a simple example of a shrinkage estimator: one "shrinks" the unbiased estimator towards zero. For the normal distribution, dividing by *n* + 1 (instead of *n* − 1 or *n*) minimizes mean squared error. The resulting estimator is biased, however, and is known as the **biased sample variation**.

## Population variance

In general, the * population variance* of a

*finite*population of size

*N*with values

*x*

_{i}is given by

where the population mean is

The population variance can also be computed using

This is true because

The population variance matches the variance of the generating probability distribution. In this sense, the concept of population can be extended to continuous random variables with infinite populations.

## Sample variance

In many practical situations, the true variance of a population is not known *a priori* and must be computed somehow. When dealing with extremely large populations, it is not possible to count every object in the population, so the computation must be performed on a sample of the population. Sample variance can also be applied to the estimation of the variance of a continuous distribution from a sample of that distribution.

We take a sample with replacement of *n* values *y*_{1}, ..., *y*_{n} from the population, where *n* < *N*, and estimate the variance on the basis of this sample. Directly taking the variance of the sample data gives the average of the squared deviations:

Here,

Since the *y*_{i} are selected randomly, both *y*_{i}} from the population. For

Hence *biased sample variance*. Correcting for this bias yields the *unbiased sample variance*:

Either estimator may be simply referred to as the *sample variance* when the version can be determined by context. The same proof is also applicable for samples taken from a continuous probability distribution.

The use of the term *n* − 1 is called Bessel's correction, and it is also used in sample covariance and the sample standard deviation (the square root of variance). The square root is a concave function and thus introduces negative bias (by Jensen's inequality), which depends on the distribution, and thus the corrected sample standard deviation (using Bessel's correction) is biased. The unbiased estimation of standard deviation is a technically involved problem, though for the normal distribution using the term *n* − 1.5 yields an almost unbiased estimator.

The unbiased sample variance is a U-statistic for the function *ƒ*(*y*_{1}, *y*_{2}) = (*y*_{1} − *y*_{2})^{2}/2, meaning that it is obtained by averaging a 2-sample statistic over 2-element subsets of the population.

## Distribution of the sample variance

Being a function of random variables, the sample variance is itself a random variable, and it is natural to study its distribution. In the case that *y*_{i} are independent observations from a normal distribution, Cochran's theorem shows that *s*^{2} follows a scaled chi-squared distribution:

As a direct consequence, it follows that

and

If the *y*_{i} are independent and identically distributed, but not necessarily normally distributed, then

where *κ* is the kurtosis of the distribution and *μ _{4}* is the fourth central moment.

If the conditions of the law of large numbers hold for the squared observations, *s*^{2} is a consistent estimator of *σ*^{2}. One can see indeed that the variance of the estimator tends asymptotically to zero. An asymptotically equivalent formula was given in Kenney and Keeping (1951:164), Rose and Smith (2002:264), and Weisstein (n.d.).

## Samuelson's inequality

Samuelson's inequality is a result that states bounds on the values that individual observations in a sample can take, given that the sample mean and (biased) variance have been calculated. Values must lie within the limits

## Relations with the harmonic and arithmetic means

It has been shown that for a sample {*y*_{i}} of real numbers,

where *y*_{max} is the maximum of the sample, *A* is the arithmetic mean, *H* is the harmonic mean of the sample and

This bound has been improved, and it is known that variance is bounded by

where *y*_{min} is the minimum of the sample.

## Tests of equality of variances

Testing for the equality of two or more variances is difficult. The F test and chi square tests are both adversely affected by non-normality and are not recommended for this purpose.

Several non parametric tests have been proposed: these include the Barton–David–Ansari–Freund–Siegel–Tukey test, the Capon test, Mood test, the Klotz test and the Sukhatme test. The Sukhatme test applies to two variances and requires that both medians be known and equal to zero. The Mood, Klotz, Capon and Barton–David–Ansari–Freund–Siegel–Tukey tests also apply to two variances. They allow the median to be unknown but do require that the two medians are equal.

The Lehmann test is a parametric test of two variances. Of this test there are several variants known. Other tests of the equality of variances include the Box test, the Box–Anderson test and the Moses test.

Resampling methods, which include the bootstrap and the jackknife, may be used to test the equality of variances.

## History

The term *variance* was first introduced by Ronald Fisher in his 1918 paper *The Correlation Between Relatives on the Supposition of Mendelian Inheritance*:

The great body of available statistics show us that the deviations of a human measurement from its mean follow very closely the Normal Law of Errors, and, therefore, that the variability may be uniformly measured by the standard deviation corresponding to the square root of the mean square error. When there are two independent causes of variability capable of producing in an otherwise uniform population distributions with standard deviations

## Moment of inertia

The variance of a probability distribution is analogous to the moment of inertia in classical mechanics of a corresponding mass distribution along a line, with respect to rotation about its center of mass. It is because of this analogy that such things as the variance are called *moments* of probability distributions. The covariance matrix is related to the moment of inertia tensor for multivariate distributions. The moment of inertia of a cloud of *n* points with a covariance matrix of

This difference between moment of inertia in physics and in statistics is clear for points that are gathered along a line. Suppose many points are close to the *x* axis and distributed along it. The covariance matrix might look like

That is, there is the most variance in the *x* direction. Physicists would consider this to have a low moment *about* the *x* axis so the moment-of-inertia tensor is

## Semivariance

The *semivariance* is calculated in the same manner as the variance but only those observations that fall below the mean are included in the calculation. It is sometimes described as a measure of downside risk in an investments context. For skewed distributions, the semivariance can provide additional information that a variance does not.

## Generalizations

If

If

If