# Variance

In probability theory and statistics, the variance is a measure of how far a set of numbers is spread out. It is one of several descriptors of a probability distribution, describing how far the numbers lie from the mean (expected value). In particular, the variance is one of the moments of a distribution. In that context, it forms part of a systematic approach to distinguishing between probability distributions. While other such approaches have been developed, those based on moments are advantageous in terms of mathematical and computational simplicity.

The variance is a parameter describing in part either the actual probability distribution of an observed population of numbers, or the theoretical probability distribution of a sample (a not-fully-observed population) of numbers. In the latter case a sample of data from such a distribution can be used to construct an estimate of its variance: in the simplest cases this estimate can be the sample variance.

## Definition

The variance of a random variable X is its second central moment, the expected value of the squared deviation from the mean μ = E[X]:

$\operatorname{Var}(X) = \operatorname{E}\left[(X - \mu)^2 \right].$

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

$\operatorname{Var}(X) = \operatorname{Cov}(X, X).$

The variance is also equivalent to the second cumulant of the probability distribution for X. The variance is typically designated as Var(X), $\scriptstyle\sigma_X^2$, or simply σ2 (pronounced "sigma squared"). The expression for the variance can be expanded:

$\operatorname{Var}(X)= \operatorname{E}\left[X^2 - 2X\operatorname{E}[X] + (\operatorname{E}[X])^2\right] = \operatorname{E}\left[X^2\right] - 2\operatorname{E}[X]\operatorname{E}[X] + (\operatorname{E}[X])^2 = \operatorname{E}\left[X^2 \right] - (\operatorname{E}[X])^2$

A mnemonic for the above expression is "mean of square minus square of mean".

### Continuous random variable

If the random variable X is continuous with probability density function f(x), then the variance is given by

$\operatorname{Var}(X) =\sigma^2 =\int (x-\mu)^2 \, f(x) \, dx\, =\int x^2 \, f(x) \, dx\, - \mu^2$

where $\mu$ is the expected value,

$\mu = \int x \, f(x) \, dx\,$

and where the integrals are definite integrals taken for x ranging over the range of X.

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 k satisfies 1 < k ≤ 2.

### Discrete random variable

If the random variable X is discrete with probability mass function x1 ↦ p1, ..., xn ↦ pn, then

$\operatorname{Var}(X) = \sum_{i=1}^n p_i\cdot(x_i - \mu)^2 = \sum_{i=1}^n (p_i\cdot x_i^2) - \mu^2$

where $\mu$ is the expected value, i.e.

$\mu = \sum_{i=1}^n p_i\cdot x_i$ .

(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 n equally likely values can be written as

$\operatorname{Var}(X) = \frac{1}{n} \sum_{i=1}^n (x_i - \mu)^2.$

The variance of a set of n equally likely values can be equivalently expressed, without directly referring to the mean, in terms of squared deviations of all points from each other:

$\operatorname{Var}(X) = \frac{1}{n^2} \sum_{i=1}^n \sum_{j=1}^n \frac{1}{2}(x_i - x_j)^2.$
↑Jump back a section

## Examples

### Normal distribution

The normal distribution with parameters μ and σ is a continuous distribution whose probability density function is given by:

$f(x) = \frac{1}{\sqrt{2\pi \sigma^2}} e^{ -\frac{(x-\mu)^2}{2\sigma^2} }.$

It has mean μ and variance equal to:

$\operatorname{Var}(X) = \int_{-\infty}^\infty \frac{(x - \mu)^2}{\sqrt{2\pi \sigma^2}} e^{ -\frac{(x-\mu)^2}{2\sigma^2} } \, dx = \sigma^2.$

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 λ is a continuous distribution whose support is the semi-infinite interval [0,∞). Its probability density function is given by:

$f(x) = \lambda e^{-\lambda x},\,$

and it has expected value μ = λ−1. The variance is equal to:

$\operatorname{Var}(X) = \int_0^\infty (x - \lambda^{-1})^2 \, \lambda e^{-\lambda x} dx = \lambda^{-2}.\,$

So for an exponentially distributed random variable σ2 = μ2.

### Poisson distribution

The Poisson distribution with parameter λ is a discrete distribution for k = 0, 1, 2, ... Its probability mass function is given by:

$p(k) = \frac{\lambda^k}{k!} e^{-\lambda k},$

and it has expected value μ = λ. The variance is equal to:

$\operatorname{Var}(X) = \sum_{k=1}^{n} \frac{\lambda^k}{k!} e^{-\lambda k} (k-\lambda)^2 = \lambda,$

So for a Poisson-distributed random variable σ2 = μ.

### Binomial distribution

The binomial distribution with parameters n and p is a discrete distribution for k = 0, 1, 2, ... Its probability mass function is given by:

$p(k) = {n\choose k}p^k(1-p)^{n-k},$

and it has expected value μ = np. The variance is equal to:

$\operatorname{Var}(X) = \sum_{k=1}^{n} {n\choose k}p^k(1-p)^{n-k} (k-np)^2 = np(1-p),$

#### Coin Toss

The binomial distribution with p = 0.5 describes the probability of getting k heads in n tosses. Thus the expected value of the number of heads is n/2, and the variance is n/4.

### Fair die

A six-sided fair die can be modelled with a discrete random variable with outcomes 1 through 6, each with equal probability $\textstyle\frac{1}{6}$. The expected value is (1 + 2 + 3 + 4 + 5 + 6)/6 = 3.5. Therefore the variance can be computed to be:

\begin{align} \sum_{i=1}^6 \tfrac{1}{6}(i - 3.5)^2 = \tfrac{1}{6}\sum_{i=1}^6 (i - 3.5)^2 & = \tfrac{1}{6}\left((-2.5)^2{+}(-1.5)^2{+}(-0.5)^2{+}0.5^2{+}1.5^2{+}2.5^2\right) \\ & = \tfrac{1}{6} \cdot 17.50 = \tfrac{35}{12} \approx 2.92. \end{align}

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

\begin{align} \sigma^2=E(X^2)-(E(X))^2 &=\frac{1}{n}\sum_{i=1}^n i^2-\left(\frac{1}{n}\sum_{i=1}^n i\right)^2 \\ &=\tfrac 16 (n+1)(2n+1) - \tfrac 14 (n+1)^2\\ &=\frac{ n^2-1 }{12}. \end{align}
↑Jump back a section

## Properties

### Basic properties

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

$\operatorname{Var}(X)\ge 0.$

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.

$P(X=a) = 1\Leftrightarrow \operatorname{Var}(X)= 0.$

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.

$\operatorname{Var}(X+a)=\operatorname{Var}(X).$

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

$\operatorname{Var}(aX)=a^2\operatorname{Var}(X).$

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

$\operatorname{Var}(aX+bY)=a^2\operatorname{Var}(X)+b^2\operatorname{Var}(Y)+2ab\, \operatorname{Cov}(X,Y),$
$\operatorname{Var}(X-Y)=\operatorname{Var}(X)+\operatorname{Var}(Y)-2\, \operatorname{Cov}(X,Y),$

In general we have for the sum of $N$ random variables:

$\operatorname{Var}\left(\sum_{i=1}^N X_i\right)=\sum_{i,j=1}^N\operatorname{Cov}(X_i,X_j)=\sum_{i=1}^N\operatorname{Var}(X_i)+\sum_{i\ne j}\operatorname{Cov}(X_i,X_j).$

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

\begin{align} \operatorname{Var}\left( \sum_{i=1}^{N} a_iX_i\right) &=\sum_{i,j=1}^{N} a_ia_j\operatorname{Cov}(X_i,X_j) \\ &=\sum_{i=1}^{N}a_i^2\operatorname{Var}(X_i)+\sum_{i\not=j}a_ia_j\operatorname{Cov}(X_i,X_j)\\ & =\sum_{i=1}^{N}a_i^2\operatorname{Var}(X_i)+2\sum_{1\le i

The variance of a finite sum of uncorrelated random variables is equal to the sum of their variances. This stems from the above identity and the fact that for uncorrelated variables the covariance is zero; that is, if

$\operatorname{Cov}(X_i,X_j)=0\ (i\ne j) ,$

then

$\operatorname{Var}\left(\sum_{i=1}^N X_i\right)=\sum_{i=1}^N\operatorname{Var}(X_i).$

### 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:

$\operatorname{Var}\Big(\sum_{i=1}^n X_i\Big) = \sum_{i=1}^n \operatorname{Var}(X_i).$

This statement is called the Bienaymé formula.[1] and was discovered in 1853.[citation needed] It is often made with the stronger condition that the variables are independent, but uncorrelatedness 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

$\operatorname{Var}\left(\overline{X}\right) = \operatorname{Var}\left(\frac {1} {n}\sum_{i=1}^n X_i\right) = \frac {1} {n^2}\sum_{i=1}^n \operatorname{Var}\left(X_i\right) = \frac {\sigma^2} {n}.$

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.

### Product of independent variables

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

$\operatorname{Var}(XY) = [E(X)]^{2}\operatorname{Var}(Y) + [E(Y)]^{2}\operatorname{Var}(X) + \operatorname{Var}(X)\operatorname{Var}(Y).$

### Sum of correlated variables

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

$\operatorname{Var}\left(\sum_{i=1}^n X_i\right) = \sum_{i=1}^n \sum_{j=1}^n \operatorname{Cov}(X_i, X_j).$

(Note: This by definition includes the variance of each variable, since Cov(Xi,Xi) = Var(Xi).)

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. 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

$\operatorname{Var}(\overline{X}) = \frac {\sigma^2} {n} + \frac {n-1} {n} \rho \sigma^2.$

This implies that the variance of the mean increases with the average of the correlations. Moreover, if the variables have unit variance, for example if they are standardized, then this simplifies to

$\operatorname{Var}(\overline{X}) = \frac {1} {n} + \frac {n-1} {n} \rho.$

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

$\lim_{n \to \infty} \operatorname{Var}(\overline{X}) = \rho.$

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 generally not converge to the population mean, even though the Law of large numbers states that the sample mean will converge for independent variables.

### Weighted sum of variables

The scaling property and the Bienaymé formula, along with this property from the covariance page: Cov(aXbY) = ab Cov(XY) jointly imply that

$\operatorname{Var}(aX+bY) =a^2 \operatorname{Var}(X) + b^2 \operatorname{Var}(Y) + 2ab\, \operatorname{Cov}(X, Y).$

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:

$\operatorname{Var}\left(\sum_{i}^n a_iX_i\right) = \sum_{i=1}^na_i^2 \operatorname{Var}(X_i) + 2\sum_{1\le i}\sum_{

### Decomposition

The general formula for variance decomposition or the law of total variance is: If $X$ and $Y$ are two random variables and the variance of $X$ exists, then

$\operatorname{Var}(X) = \operatorname{Var}(\operatorname{E}(X|Y))+ \operatorname{E}(\operatorname{Var}(X|Y)).$

Here, $\operatorname E(X|Y)$ is the conditional expectation of $X$ given $Y$, and $\operatorname{Var}(X|Y)$ is the conditional variance of $X$ given $Y$. (A more intuitive explanation is that given a particular value of $Y$, then $X$ follows a distribution with mean $\operatorname E(X|Y)$ and variance $\operatorname{Var}(X|Y)$. The above formula tells how to find $\operatorname{Var}(X)$ based on the distributions of these two quantities when $Y$ is allowed to vary.) This formula is often applied in analysis of variance, where the corresponding formula is

$\mathit{MS}_\mathrm{Total} = \mathit{MS}_\mathrm{Between} + \mathit{MS}_\mathrm{Within};$

here $\mathit{MS}$ refers to the Mean of the Squares. It is also used in linear regression analysis, where the corresponding formula is

$\mathit{MS}_\mathrm{Total} = \mathit{MS}_\mathrm{Regression} + \mathit{MS}_\mathrm{Residual}.$

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, $\mathit{SS}$):

$\mathit{SS}_\mathrm{Total} = \mathit{SS}_\mathrm{Between} + \mathit{SS}_\mathrm{Within},$
$\mathit{SS}_\mathrm{Total} = \mathit{SS}_\mathrm{Regression} + \mathit{SS}_\mathrm{Residual}.$

### Formulae for the variance

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

$\operatorname{Var}(X) =\operatorname{E}(X^2) - (\operatorname{E}(X))^2.$

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.[citation needed] This is discussed below.

### 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

$2\int_0^\infty uH(u)\,du - \Big(\int_0^\infty H(u)\,du\Big)^2.$

where H(u) = 1 − F(u) is the right tail function. 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. $\mathrm{argmin}_m\,\mathrm{E}((X - m)^2) = \mathrm{E}(X)\,$. Conversely, if a continuous function $\varphi$ satisfies $\mathrm{argmin}_m\,\mathrm{E}(\varphi(X - m)) = \mathrm{E}(X)\,$ for all random variables X, then it is necessarily of the form $\varphi(x) = a x^2 + b$, where a > 0. This also holds in the multidimensional case.[4]

### Matrix notation for the variance of a linear combination

Let's define $X$ as a column vector of n random variables $X_1, ...,X_n$, and c as a column vector of N scalars $c_1, ...,c_n$. Therefore $c^T X$ is a linear combination of these random variables, where $c^T$ denotes the transpose of vector $c$. Let also be $\Sigma$ the variance-covariance matrix of the vector X. The variance of $c^TX$ is given by:[5]

$\operatorname{Var}(c^T X) = c^T \Sigma c .$

### 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 inches will have a variance measured in square inches. 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.

↑Jump back a section

## 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

$\operatorname{Var}\left[f(X)\right]\approx \left(f'(\operatorname{E}\left[X\right])\right)^2\operatorname{Var}\left[X\right]$

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

↑Jump back a section

## Population variance and sample variance

Real-world distributions such as the distribution of yesterday's rain throughout the day are typically not fully known, unlike the behavior of perfect dice or an ideal distribution such as the normal distribution, because it is impractical to account for every raindrop. Instead one estimates the mean and variance of the whole distribution as the computed mean and variance of a sample of n observations drawn suitably randomly from the whole sample space, in this example the set of all measurements of yesterday's rainfall in all available rain gauges.

This method of estimation is close to optimal, with the caveat that it underestimates the variance by a factor of (n − 1) / n. (For example, when n = 1 the variance of a single observation is obviously zero regardless of the true variance). This gives a bias which should be corrected for when n is small by multiplying by n / (n − 1). 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.

### Population Variance

In general, the population variance of a finite population of size N with values xi is given by

$\sigma^2 = \frac 1N \sum_{i=1}^N \left(x_i - \mu \right)^2 = \left(\frac 1N \sum_{i=1}^N x_i^2\right) - \mu^2$

where

$\mu = \frac 1N \sum_{i=1}^N x_i$

is the population mean. The population variance therefore is the variance of the underlying 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.[6] 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 y1, ..., yn from the population, where n < N, and estimate the variance on the basis of this sample.[7] Directly taking the variance of the sample gives:

$\sigma_y^2 = \frac 1n \sum_{i=1}^n \left(y_i - \overline{y} \right)^2$

Here, $\overline{y}$ denotes the sample mean:

$\overline{y}=\frac 1n \sum_{i=1}^n y_i .$

Since the yi are selected randomly, both $\scriptstyle\overline{y}$ and $\scriptstyle\sigma_y^2$ are random variables. Their expected values can be evaluated by summing over the ensemble of all possible samples {yi} from the population. For $\scriptstyle\sigma_y^2$ this gives:

\begin{align} E[\sigma_y^2] & = E\left[ \frac 1n \sum_{i=1}^n \left(y_i - \frac 1n \sum_{j=1}^n y_j \right)^2 \right] \\ & = \frac 1n \sum_{i=1}^n E\left[ y_i^2 - \frac 2n y_i \sum_{j=1}^n y_j + \frac{1}{n^2} \sum_{j=1}^n y_j \sum_{k=1}^n y_k \right] \\ & = \frac 1n \sum_{i=1}^n \left[ \frac{n-2}{n} E[y_i^2] - \frac 2n \sum_{j \neq i} E[y_i y_j] + \frac{1}{n^2} \sum_{j=1}^n \sum_{k \neq j} E[y_j y_k] +\frac{1}{n^2} \sum_{j=1}^n E[y_j^2] \right] \\ & = \frac 1n \sum_{i=1}^n \left[ \frac{n-2}{n} (\sigma^2+\mu^2) - \frac 2n (n-1) \mu^2 + \frac{1}{n^2} n (n-1) \mu^2 + \frac 1n (\sigma^2+\mu^2) \right] \\ & = \frac{n-1}{n} \sigma^2. \end{align}

Hence $\scriptstyle\sigma_y^2$ gives an estimate of the population variance that is biased by a factor of (n-1)/n. For this reason, $\scriptstyle\sigma_y^2$ is referred to as the biased sample variance. Correcting for this bias yields the unbiased sample variance:

$s^2 = \frac{1}{n-1} \sum_{i=1}^n \left(y_i - \overline{y} \right)^2$

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 ƒ(y1y2) = (y1 − y2)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

Distribution and cumulative distribution of s22, for various values of ν = n-1, when the yi are independent normally distributed.

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 yi are independent observations from a normal distribution, Cochran's theorem shows that s2 follows a scaled chi-squared distribution:[8]

$(n-1)\frac{s^2}{\sigma^2}\sim\chi^2_{n-1}.$

As a direct consequence, it follows that

$\operatorname{E}(s^2)=\operatorname{E}\left(\frac{\sigma^2}{n-1} \chi^2_{n-1}\right)=\sigma^2 ,$

and[9]

$\operatorname{Var}[s^2] =\operatorname{Var}\left(\frac{\sigma^2}{n-1} \chi^2_{n-1}\right)=\frac{\sigma^4}{(n-1)^2}\operatorname{Var}\left( \chi^2_{n-1}\right)=\frac{2\sigma^4 }{n-1}.$

If the yi are independent and identically distributed, but not necessarily normally distributed, then[10]

$\operatorname{E}[s^2] = \sigma^2, \quad \operatorname{Var}[s^2] = \sigma^4 \left (\frac{2}{n-1} + \frac{\kappa}{n} \right) = \frac{1}{n} \left(\mu_4 - \frac{n-3}{n-1}\sigma^4\right),$

where κ is the excess kurtosis of the distribution and μ4 is the fourth moment about the mean.

If the conditions of the law of large numbers hold for the squared observations, s2 is a consistent estimator of σ2.[citation needed]. One can see indeed that the variance of the estimator tends asymptotically to zero.

### 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.[11] Values must lie within the limits $\bar y \pm \sigma_y (n-1)^{1/2}.$

### Relations with the harmonic and arithmetic means

It has been shown[12] that for a sample {yi} of real numbers,

$\sigma_y^2 \le 2y_{max} (A - H)$

where ymax is the maximum of the sample, A is the arithmetic mean, H is the harmonic mean of the sample and $\sigma_y^2$ is the (biased) variance of the sample.

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

$\sigma_y^2 \le \frac{ y_{max} ( A - H )( y_{max} - A ) } { y_{max} - H },$
$\sigma_y^2 \ge \frac{ y_{min} ( A - H )( A - y_{min} ) } { H - y_{min} },$

where ymin is the minimum of the sample.[13]

↑Jump back a section

## Generalizations

If $X$ is a vector-valued random variable, with values in $\mathbb{R}^n$, and thought of as a column vector, then the natural generalization of variance is $\operatorname{E}((X - \mu)(X - \mu)^{\operatorname{T}})$, where $\mu = \operatorname{E}(X)$ and $X^{\operatorname{T}}$ is the transpose of $X$, and so is a row vector. This variance is a positive semi-definite square matrix, commonly referred to as the covariance matrix.

If $X$ is a complex-valued random variable, with values in $\mathbb{C}$, then its variance is $\operatorname{E}((X - \mu)(X - \mu)^{\dagger})$, where $X^{\dagger}$ is the conjugate transpose of $X$. This variance is also a positive semi-definite square matrix.

↑Jump back a section

## 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 sensitive to non normality and are not recommended for this purpose.

Several non parametric tests have been proposed: these include the Barton-David-Ansari-Fruend-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-Fruend-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 Lehman 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.

↑Jump back a section

## History

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

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 $\theta_1$ and $\theta_2$, it is found that the distribution, when both causes act together, has a standard deviation $\sqrt{\theta_1^2 + \theta_2^2}$. It is therefore desirable in analysing the causes of variability to deal with the square of the standard deviation as the measure of variability. We shall term this quantity the Variance...

↑Jump back a section

## 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.[citation needed] It is because of this analogy that such things as the variance are called moments of probability distributions.[citation needed] 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 $\Sigma$ is given by[citation needed]

$I=n (\mathbf{1}_{3\times 3} \operatorname{tr}(\Sigma) - \Sigma).$

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

$\Sigma=\begin{bmatrix}10 & 0 & 0\\0 & 0.1 & 0 \\ 0 & 0 & 0.1\end{bmatrix}.$

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

$I=n\begin{bmatrix}0.2 & 0 & 0\\0 & 10.1 & 0 \\ 0 & 0 & 10.1\end{bmatrix}.$
↑Jump back a section

## Notes

1. ^ Loeve, M. (1977) "Probability Theory", Graduate Texts in Mathematics, Volume 45, 4th edition, Springer-Verlag, p. 12.
2. ^ Goodman, Leo A., "On the exact variance of products," Journal of the American Statistical Association, December 1960, 708–713.
3. ^ Goodman, Leo A., "The variance of the product of K random variables," Journal of the American Statistical Association, March 1962, 54ff.
4. ^ Kagan, A.; Shepp, L. A. (1998). "Why the variance?". Statistics & Probability Letters 38 (4): 329. doi:10.1016/S0167-7152(98)00041-8. edit
5. ^ Johnson, Richard; Wichern, Dean (2001). Applied Multivariate Statistical Analysis. Prentice Hall. p. 76. ISBN 0-13-187715-1
6. ^ Navidi, William (2006) Statistics for Engineers and Scientists, McGraw-Hill, pg 14.
7. ^ Montgomery, D. C. and Runger, G. C. (1994) Applied statistics and probability for engineers, page 201. John Wiley & Sons New York
8. ^ Knight K. (2000), Mathematical Statistics, Chapman and Hall, New York. (proposition 2.11)
9. ^ Casella and Berger (2002) Statistical Inference, Example 7.3.3, p. 331[full citation needed]
10. ^ Neter, Wasserman, and Kutner (1990) Applied Linear Statistical Models, 3rd edition, pp. 622-623[full citation needed]
11. ^ Samuelson, Paul (1968)"How Deviant Can You Be?", Journal of the American Statistical Association, 63, number 324 (December, 1968), pp. 1522–1525 JSTOR 2285901
12. ^ Mercer A McD (2000) Bounds for A-G, A-H, G-H, and a family of inequalities of Ky Fan’s type, using a general method. J Math Anal Appl 243, 163–173
13. ^ Sharma R (2008) Some more inequalities for arithmetic mean, harmonic mean and variance. J Math Inequalities 2 (1) 109–114
14. ^
↑Jump back a section