Variance
Variance
There are two distinct concepts that are both called "variance". One, as discussed above, is part of a
theoretical probability distribution and is defined by an equation. The other variance is a characteristic of a
set of observations. When variance is calculated from observations, those observations are typically
measured from a real world system. If all possible observations of the system are present then the calculated
variance is called the population variance. Normally, however, only a subset is available, and the variance
calculated from this is called the sample variance. The variance calculated from a sample is considered an
estimate of the full population variance. There are multiple ways to calculate an estimate of the population
variance, as discussed in the section below.
The two kinds of variance are closely related. To see how, consider that a theoretical probability distribution
can be used as a generator of hypothetical observations. If an infinite number of observations are generated
using a distribution, then the sample variance calculated from that infinite set will match the value
calculated using the distribution's equation for variance. Variance has a central role in statistics, where some
ideas that use it include descriptive statistics, statistical inference, hypothesis testing, goodness of fit, and
Monte Carlo sampling.
Etymology
The term variance was first introduced by Ronald Fisher in his 1918 paper The Correlation Between
Relatives on the Supposition of Mendelian Inheritance:[2]
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 and , it is
found that the distribution, when both causes act together, has a standard deviation
. 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...
Definition
The variance of a random variable is the expected value of the
squared deviation from the mean of , :
1. A frequency distribution is
constructed.
2. The centroid of the
distribution gives its mean.
3. A square with sides equal
In other words, the variance of X is equal to the mean of the square to the difference of each
of X minus the square of the mean of X. This equation should not value from the mean is
be used for computations using floating point arithmetic, because it formed for each value.
suffers from catastrophic cancellation if the two components of the 4. Arranging the squares into
equation are similar in magnitude. For other numerically stable a rectangle with one side
alternatives, see Algorithms for calculating variance. equal to the number of
values, n, results in the
other side being the
Discrete random variable distribution's variance, σ2.
(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 equally likely values can be equivalently expressed, without directly referring to
the mean, in terms of squared deviations of all pairwise squared distances of points from each other:[3]
If the random variable has a probability density function , and is the corresponding
cumulative distribution function, then
or equivalently,
where is the expected value of given by
In these formulas, the integrals with respect to and are Lebesgue and Lebesgue–Stieltjes
integrals, respectively.
Examples
Exponential distribution
The exponential distribution with parameter λ is a continuous distribution whose probability density
function is given by
Using integration by parts and making use of the expected value already calculated, we have:
Fair die
A fair six-sided die can be modeled as a discrete random variable, X, with outcomes 1 through 6, each with
equal probability 1/6. The expected value of X is Therefore, the
variance of X is
The general formula for the variance of the outcome, X, of an n -sided die is
The following table lists the variance for some commonly used probability distributions.
Name of the probability distribution Probability distribution function Mean Variance
Binomial distribution
Geometric distribution
Normal distribution
Exponential distribution
Poisson distribution
Properties
Basic properties
Conversely, if the variance of a random variable is 0, then it is almost surely a constant. That is, it always
has the same value:
Issues of finiteness
If a distribution does not have a finite expected value, as is the case for the Cauchy distribution, then the
variance cannot be finite either. However, some distributions may not have a finite variance, despite their
expected value being finite. An example is a Pareto distribution whose index satisfies
Decomposition
The general formula for variance decomposition or the law of total variance is: If and are two random
variables, and the variance of exists, then
The conditional expectation of given , and the conditional variance may be
understood as follows. Given any particular value y of the random variable Y, there is a conditional
expectation given the event Y = y. This quantity depends on the particular value y; it is a
function . That same function evaluated at the random variable Y is the conditional
expectation
here refers to the Mean of the Squares. In linear regression analysis the corresponding formula is
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, ):
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. . Conversely, if a continuous
function satisfies for all random variables X, then it is necessarily of
the form , where a > 0. This also holds in the multidimensional case.[4]
Units of measurement
Unlike the 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.
Propagation
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:
Linear combinations
In general, for the sum of random variables , the variance becomes:
then they are said to be uncorrelated. It follows immediately from the expression given earlier that if the
random variables are uncorrelated, then the variance of their sum is equal to the sum of their
variances, or, expressed symbolically:
Since independent random variables are always uncorrelated (see Covariance § Uncorrelatedness and
independence), the equation above holds in particular when the random variables are
independent. Thus, independence is sufficient but not necessary for the variance of the sum to equal the
sum of the variances.
This implies that the variance of the mean can be written as (with a column vector of ones)
Sum of variables
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[6] and was discovered in 1853.[7][8] 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.
The general result then follows by induction. Starting with the definition,
Using the linearity of the expectation operator and the assumption of independence (or uncorrelatedness) of
X and Y, this further simplifies as follows:
In general, the variance of the sum of n variables is the sum of their covariances:
(Note: The second equality comes from the fact that Cov(Xi,Xi) = Var(Xi).)
Here, 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.
There are cases when a sample is taken without knowing, in advance, how many observations will be
acceptable according to some criterion. In such cases, the sample size N is a random variable whose
variation adds to the variation of X, such that,
[9]
sample mean).
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.
Product of variables
If two variables X and Y are independent, the variance of their product is given by[10]
In general, if two variables are statistically dependent, then the variance of their product is given by:
Arbitrary functions
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 from a limited 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 true population mean is unknown, then the sample variance (which uses the sample mean in
place of the true mean) 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 xi is given by
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.[11] This is generally
referred to as sample variance or empirical variance. 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.[12] Directly taking the variance of the sample data gives the
average of the squared deviations:
Since the Yi are selected randomly, both and are random variables. Their expected values can be
evaluated by averaging over the ensemble of all possible samples {Yi} of size n from the population. For
this gives:
Hence gives an estimate of the population variance that is biased by a factor of . For this reason,
Correcting for this bias yields the unbiased sample variance, denoted :
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 ƒ(y1 , y2 ) = (y1 − y2 )2 /2, meaning that it is
obtained by averaging a 2-sample statistic over 2-element subsets of the population.
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 (see also: asymptotic properties and an
elementary proof):[13]
As a direct consequence, it follows that
and[14]
If the Yi are independent and identically distributed, but not necessarily normally distributed, then[15]
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, S2 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.).[16][17][18]
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.[19] Values must lie within
the limits
It has been shown[20] that for a sample {yi} of positive real numbers,
where ymax is the maximum of the sample, A is the arithmetic mean, H is the harmonic mean of the sample
and is the (biased) variance of the sample.
This bound has been improved, and it is known that variance is bounded by
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.
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 is given by
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 also described as a specific measure in different fields of application. For skewed distributions, the
semivariance can provide additional information that a variance does not.[22]
For inequalities associated with the semivariance, see Chebyshev's inequality § Semivariances.
Generalizations
As a matrix
If is a vector-valued random variable, with values in and thought of as a column vector, then a
natural generalization of variance is where and is the transpose
of and so is a row vector. The result is a positive semi-definite square matrix, commonly referred to as
the variance-covariance matrix (or simply as the covariance matrix).
If is a vector- and complex-valued random variable, with values in then the covariance matrix is
where is the conjugate transpose of This matrix is also positive semi-
definite and square.
As a scalar
Another generalization of variance for vector-valued random variables , which results in a scalar value
rather than in a matrix, is the generalized variance , the determinant of the covariance matrix. The
generalized variance can be shown to be related to the multidimensional scatter of points around their
mean.[23]
A different generalization is obtained by considering the variance of the Euclidean distance between the
random variable and its mean. This results in which is the trace of the
covariance matrix.
See also
Mathematics
portal
Bhatia–Davis inequality
Coefficient of variation
Homoscedasticity
Least-squares spectral analysis for computing a frequency spectrum with spectral
magnitudes in % of variance or in dB
Modern portfolio theory
Popoviciu's inequality on variances
Measures for statistical dispersion
Variance-stabilizing transformation
Types of variance
Correlation
Distance variance
Explained variance
Pooled variance
Pseudo-variance
References
1. Wasserman, Larry (2005). All of Statistics: a concise course in statistical inference. Springer
texts in statistics. p. 51. ISBN 9781441923226.
2. Ronald Fisher (1918) The correlation between relatives on the supposition of Mendelian
Inheritance (http://digital.library.adelaide.edu.au/dspace/bitstream/2440/15097/1/9.pdf)
3. Yuli Zhang, Huaiyu Wu, Lei Cheng (June 2012). Some new deformation formulas about
variance and covariance. Proceedings of 4th International Conference on Modelling,
Identification and Control(ICMIC2012). pp. 987–992.
4. Kagan, A.; Shepp, L. A. (1998). "Why the variance?". Statistics & Probability Letters. 38 (4):
329–333. doi:10.1016/S0167-7152(98)00041-8 (https://doi.org/10.1016%2FS0167-7152%2
898%2900041-8).
5. Johnson, Richard; Wichern, Dean (2001). Applied Multivariate Statistical Analysis (https://ar
chive.org/details/appliedmultivari00john_130). Prentice Hall. p. 76 (https://archive.org/detail
s/appliedmultivari00john_130/page/n96). ISBN 0-13-187715-1.
6. Loève, M. (1977) "Probability Theory", Graduate Texts in Mathematics, Volume 45, 4th
edition, Springer-Verlag, p. 12.
7. Bienaymé, I.-J. (1853) "Considérations à l'appui de la découverte de Laplace sur la loi de
probabilité dans la méthode des moindres carrés", Comptes rendus de l'Académie des
sciences Paris, 37, p. 309–317; digital copy available [1] (http://visualiseur.bnf.fr/CadresFen
etre?O=NUMM-2994&I=313)
8. Bienaymé, I.-J. (1867) "Considérations à l'appui de la découverte de Laplace sur la loi de
probabilité dans la méthode des moindres carrés", Journal de Mathématiques Pures et
Appliquées, Série 2, Tome 12, p. 158–167; digital copy available [2] (http://gallica.bnf.fr/ark:/
12148/bpt6k16411c/f166.image.n19)[3] (http://sites.mathdoc.fr/JMPA/PDF/JMPA_1867_2_1
2_A10_0.pdf)
9. Cornell, J R, and Benjamin, C A, Probability, Statistics, and Decisions for Civil Engineers,
McGraw-Hill, NY, 1970, pp.178-9.
10. Goodman, Leo A. (December 1960). "On the Exact Variance of Products". Journal of the
American Statistical Association. 55 (292): 708–713. doi:10.2307/2281592 (https://doi.org/1
0.2307%2F2281592). JSTOR 2281592 (https://www.jstor.org/stable/2281592).
11. Navidi, William (2006) Statistics for Engineers and Scientists, McGraw-Hill, p. 14.
12. Montgomery, D. C. and Runger, G. C. (1994) Applied statistics and probability for engineers,
page 201. John Wiley & Sons New York
13. Knight K. (2000), Mathematical Statistics, Chapman and Hall, New York. (proposition 2.11)
14. Casella and Berger (2002) Statistical Inference, Example 7.3.3, p. 331
15. Mood, A. M., Graybill, F. A., and Boes, D.C. (1974) Introduction to the Theory of Statistics, 3rd
Edition, McGraw-Hill, New York, p. 229
16. Kenney, John F.; Keeping, E.S. (1951) Mathematics of Statistics. Part Two. 2nd ed. D. Van
Nostrand Company, Inc. Princeton: New Jersey.
http://krishikosh.egranth.ac.in/bitstream/1/2025521/1/G2257.pdf
17. Rose, Colin; Smith, Murray D. (2002) Mathematical Statistics with Mathematica. Springer-
Verlag, New York.
http://www.mathstatica.com/book/Mathematical_Statistics_with_Mathematica.pdf
18. Weisstein, Eric W. (n.d.) Sample Variance Distribution. MathWorld—A Wolfram Web
Resource. http://mathworld.wolfram.com/SampleVarianceDistribution.html
19. Samuelson, Paul (1968). "How Deviant Can You Be?". Journal of the American Statistical
Association. 63 (324): 1522–1525. doi:10.1080/01621459.1968.10480944 (https://doi.org/1
0.1080%2F01621459.1968.10480944). JSTOR 2285901 (https://www.jstor.org/stable/22859
01).
20. 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" (https://doi.org/10.1006%2Fjmaa.1999.6688). J. Math. Anal.
Appl. 243 (1): 163–173. doi:10.1006/jmaa.1999.6688 (https://doi.org/10.1006%2Fjmaa.1999.
6688).
21. Sharma, R. (2008). "Some more inequalities for arithmetic mean, harmonic mean and
variance". Journal of Mathematical Inequalities. 2 (1): 109–114. CiteSeerX 10.1.1.551.9397
(https://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.551.9397). doi:10.7153/jmi-02-11
(https://doi.org/10.7153%2Fjmi-02-11).
22. Fama, Eugene F.; French, Kenneth R. (2010-04-21). "Q&A: Semi-Variance: A Better Risk
Measure?" (https://famafrench.dimensional.com/questions-answers/qa-semi-variance-a-bett
er-risk-measure.aspx). Fama/French Forum.
23. Kocherlakota, S.; Kocherlakota, K. (2004). "Generalized Variance". Encyclopedia of
Statistical Sciences. Wiley Online Library. doi:10.1002/0471667196.ess0869 (https://doi.org/
10.1002%2F0471667196.ess0869). ISBN 0471667196.