共分散

関: covariant

WordNet

(statistics) the mean value of the product of the deviations of two variates from their respective means
changing so that interrelations with another variable quantity or set of quantities remain unchanged

Wikipedia preview

出典(authority):フリー百科事典『ウィキペディア（Wikipedia）』「2016/01/13 10:37:56」(JST)

wiki en

In probability theory and statistics, covariance is a measure of how much two random variables change together. If the greater values of one variable mainly correspond with the greater values of the other variable, and the same holds for the smaller values, i.e., the variables tend to show similar behavior, the covariance is positive.^[1] In the opposite case, when the greater values of one variable mainly correspond to the smaller values of the other, i.e., the variables tend to show opposite behavior, the covariance is negative. The sign of the covariance therefore shows the tendency in the linear relationship between the variables. The magnitude of the covariance is not easy to interpret. The normalized version of the covariance, the correlation coefficient, however, shows by its magnitude the strength of the linear relation.

A distinction must be made between (1) the covariance of two random variables, which is a population parameter that can be seen as a property of the joint probability distribution, and (2) the sample covariance, which serves as an estimated value of the parameter.

1 Definition
2 Properties
- 2.1 A more general identity for covariance matrices
- 2.2 Uncorrelatedness and independence
- 2.3 Relationship to inner products
3 Calculating the sample covariance
4 Comments
5 Applications
- 5.1 In genetics and molecular biology
- 5.2 In financial economics
- 5.3 In meteorological data assimilation
6 See also
7 References
8 External links

Definition

The covariance between two jointly distributed real-valued random variables X and Y with finite second moments is defined as^[2]

\operatorname{cov}(X,Y) = \operatorname{E}{\big[(X - \operatorname{E}[X])(Y - \operatorname{E}[Y])\big]},

where E[X] is the expected value of X, also known as the mean of X. By using the linearity property of expectations, this can be simplified to

\begin{align} \operatorname{cov}(X,Y) &= \operatorname{E}\left[\left(X - \operatorname{E}\left[X\right]\right) \left(Y - \operatorname{E}\left[Y\right]\right)\right] \\ &= \operatorname{E}\left[X Y - X \operatorname{E}\left[Y\right] - \operatorname{E}\left[X\right] Y + \operatorname{E}\left[X\right] \operatorname{E}\left[Y\right]\right] \\ &= \operatorname{E}\left[X Y\right] - \operatorname{E}\left[X\right] \operatorname{E}\left[Y\right] - \operatorname{E}\left[X\right] \operatorname{E}\left[Y\right] + \operatorname{E}\left[X\right] \operatorname{E}\left[Y\right] \\ &= \operatorname{E}\left[X Y\right] - \operatorname{E}\left[X\right] \operatorname{E}\left[Y\right]. \end{align}

However, when , this last equation is prone to catastrophic cancellation when computed with floating point arithmetic and thus should be avoided in computer programs when the data has not been centered before.^[3] Numerically stable algorithms should be preferred in this case.

For random vectors and , the m×n cross covariance matrix (also known as dispersion matrix or variance–covariance matrix,^[4] or simply called covariance matrix) is equal to

\begin{align} \operatorname{cov}(\mathbf{X},\mathbf{Y}) & = \operatorname{E} \left[(\mathbf{X} - \operatorname{E}[\mathbf{X}]) (\mathbf{Y} - \operatorname{E}[\mathbf{Y}])^\mathrm{T}\right]\\ & = \operatorname{E}\left[\mathbf{X} \mathbf{Y}^\mathrm{T}\right] - \operatorname{E}[\mathbf{X}]\operatorname{E}[\mathbf{Y}]^\mathrm{T}, \end{align}

where m^T is the transpose of the vector (or matrix) m.

The (i,j)-th element of this matrix is equal to the covariance Cov(X_i, Y_j) between the i-th scalar component of X and the j-th scalar component of Y. In particular, Cov(Y, X) is the transpose of Cov(X, Y).

For a vector $\mathbf{X}= \begin{bmatrix}X_1 & X_2 & \dots & X_m\end{bmatrix}^\mathrm{T}$ of m jointly distributed random variables with finite second moments, its covariance matrix is defined as

Random variables whose covariance is zero are called uncorrelated. Similarly, random vectors whose covariance matrix is zero in every entry outside the main diagonal are called uncorrelated.

The units of measurement of the covariance Cov(X, Y) are those of X times those of Y. By contrast, correlation coefficients, which depend on the covariance, are a dimensionless measure of linear dependence. (In fact, correlation coefficients can simply be understood as a normalized version of covariance.)

The covariance of two discrete sets can be equivalently expressed, without directly referring to the mean^[5]

Properties

Variance is a special case of the covariance when the two variables are identical:

If X, Y, W, and V are real-valued random variables and a, b, c, d are constant ("constant" in this context means non-random), then the following facts are a consequence of the definition of covariance:

\begin{align} \sigma(X, a) &= 0 \\ \sigma(X, X) &= \sigma^2(X) \\ \sigma(X, Y) &= \sigma(Y, X) \\ \sigma(aX, bY) &= ab\, \sigma(X, Y) \\ \sigma(X+a, Y+b) &= \sigma(X, Y) \\ \sigma(aX+bY, cW+dV) &= ac\,\sigma(X,W)+ad\,\sigma(X,V)+bc\,\sigma(Y,W)+bd\,\sigma(Y,V) \end{align}

For a sequence X₁, ..., X_n of random variables, and constants a₁, ..., a_n, we have

\sigma^2\left(\sum_{i=1}^n a_iX_i \right) = \sum_{i=1}^n a_i^2\sigma^2(X_i) + 2\sum_{i,j\,:\,i<j} a_ia_j\sigma(X_i,X_j) = \sum_{i,j} {a_ia_j\sigma(X_i,X_j)}

A more general identity for covariance matrices

Let be a random vector with covariance matrix , and let be a matrix that can act on . The covariance matrix of the vector is:

\Sigma(A\mathbf{X}) = A\, \Sigma(\mathbf{X})\, A^\mathrm{T}

.

This is a direct result of the linearity of expectation and is useful when applying a linear transformation, such as a whitening transformation, to a vector.

Uncorrelatedness and independence

If X and Y are independent, then their covariance is zero. This follows because under independence,

The converse, however, is not generally true. For example, let X be uniformly distributed in [-1, 1] and let Y = X². Clearly, X and Y are dependent, but

\begin{align} \sigma(X, Y) &= \sigma(X, X^2) \\ &= \operatorname{E}\!\left[X \cdot X^2\right] - \operatorname{E}[X] \cdot \operatorname{E}\!\left[X^2\right] \\ &= \operatorname{E}\!\left[X^3\right] - \operatorname{E}[X]\operatorname{E}\!\left[X^2\right] \\ &= 0 - 0 \cdot \operatorname{E}\!\left[X^2\right] \\ &= 0. \end{align}

In this case, the relationship between Y and X is non-linear, while correlation and covariance are measures of linear dependence between two variables. This example shows that if two variables are uncorrelated, that does not in general imply that they are independent. However, if two variables are jointly normally distributed (but not if they are merely individually normally distributed), uncorrelatedness does imply independence.

Relationship to inner products

Many of the properties of covariance can be extracted elegantly by observing that it satisfies similar properties to those of an inner product:

bilinear: for constants a and b and random variables X, Y, Z, σ(aX + bY, Z) = a σ(X, Z) + b σ(Y, Z);
symmetric: σ(X, Y) = σ(Y, X);
positive semi-definite: σ²(X) = σ(X, X) ≥ 0 for all random variables X, and σ(X, X) = 0 implies that X is a constant random variable (K).

In fact these properties imply that the covariance defines an inner product over the quotient vector space obtained by taking the subspace of random variables with finite second moment and identifying any two that differ by a constant. (This identification turns the positive semi-definiteness above into positive definiteness.) That quotient vector space is isomorphic to the subspace of random variables with finite second moment and mean zero; on that subspace, the covariance is exactly the L² inner product of real-valued functions on the sample space.

As a result for random variables with finite variance, the inequality

holds via the Cauchy–Schwarz inequality.

Proof: If σ²(Y) = 0, then it holds trivially. Otherwise, let random variable

Then we have

\begin{align} 0 \le \sigma^2(Z) & = \sigma\left(X - \frac{\sigma(X,Y)}{\sigma^2(Y)} Y, X - \frac{\sigma(X,Y)}{\sigma^2(Y)} Y \right) \\[12pt] & = \sigma^2(X) - \frac{ (\sigma(X,Y))^2 }{\sigma^2(Y)}. \end{align}

Calculating the sample covariance

The sample covariance of N observations of K variables is the K-by-K matrix with the entries

,

which is an estimate of the covariance between variable $j$ and variable $k$ .

The sample mean and the sample covariance matrix are unbiased estimates of the mean and the covariance matrix of the random vector , a row vector whose jth element (j = 1, ..., K) is one of the random variables. The reason the sample covariance matrix has in the denominator rather than is essentially that the population mean is not known and is replaced by the sample mean . If the population mean is known, the analogous unbiased estimate is given by

Comments

The covariance is sometimes called a measure of "linear dependence" between the two random variables. That does not mean the same thing as in the context of linear algebra (see linear dependence). When the covariance is normalized, one obtains the correlation coefficient. From it, one can obtain the Pearson coefficient, which gives the goodness of the fit for the best possible linear function describing the relation between the variables. In this sense covariance is a linear gauge of dependence.

Applications

In genetics and molecular biology

Covariance is an important measure in biology. Certain sequences of DNA are conserved more than others among species, and thus to study secondary and tertiary structures of proteins, or of RNA structures, sequences are compared in closely related species. If sequence changes are found or no changes at all are found in noncoding RNA (such as microRNA), sequences are found to be necessary for common structural motifs, such as an RNA loop.

In financial economics

Covariances play a key role in financial economics, especially in portfolio theory and in the capital asset pricing model. Covariances among various assets' returns are used to determine, under certain assumptions, the relative amounts of different assets that investors should (in a normative analysis) or are predicted to (in a positive analysis) choose to hold in a context of diversification.

In meteorological data assimilation

The covariance matrix is important in estimating the initial conditions required for running weather forecast models.

References

^ http://mathworld.wolfram.com/Covariance.html
^ Oxford Dictionary of Statistics, Oxford University Press, 2002, p. 104.
^ Donald E. Knuth (1998). The Art of Computer Programming, volume 2: Seminumerical Algorithms, 3rd edn., p. 232. Boston: Addison-Wesley.
^ W. J. Krzanowski, Principles of Multivariate Analysis, Chap. 7.1, Oxford University Press, New York, 1988
^ 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.

External links

Hazewinkel, Michiel, ed. (2001), "Covariance", Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4
MathWorld page on calculating the sample covariance
Covariance Tutorial using R
Covariance and Correlation

UpToDate Contents

全文を閲覧するには購読必要です。 To read the full text you will need to subscribe.

1. 慢性鼻副鼻腔炎のマネージメント management of chronic rhinosinusitis

English Journal

Extensive green roof CO2 exchange and its seasonal variation quantified by eddy covariance measurements.

Heusinger J1, Weber S2.
The Science of the total environment.Sci Total Environ.2017 Dec 31;607-608:623-632. doi: 10.1016/j.scitotenv.2017.07.052. Epub 2017 Jul 27.
PMID 28709096

Per- and poly-fluoroalkyl compounds in freshwater fish from the Rhône River: Influence of fish size, diet, prey contamination and biotransformation.

Babut M1, Labadie P2, Simonnet-Laprade C3, Munoz G3, Roger MC4, Ferrari BJD5, Budzinski H2, Sivade E6.
The Science of the total environment.Sci Total Environ.2017 Dec 15;605-606:38-47. doi: 10.1016/j.scitotenv.2017.06.111. Epub 2017 Jun 24.
PMID 28654807

Bayesian inference for biomarker discovery in proteomics: an analytic solution.

Dridi N1,2, Giremus A1, Giovannelli JF3, Truntzer C4, Hadzagic M1,5, Charrier JP6, Gerfault L7,8, Ducoroy P4, Lacroix B6, Grangeat P7,8, Roy P9,10,11,12.
EURASIP journal on bioinformatics & systems biology.EURASIP J Bioinform Syst Biol.2017 Dec;2017(1):9. doi: 10.1186/s13637-017-0062-4. Epub 2017 Jul 14.
PMID 28710702