Explained variation

From HandWiki

In statistics, explained variation measures the proportion to which a mathematical model accounts for the variation (dispersion) of a given data set. Often, variation is quantified as variance; then, the more specific term explained variance can be used. The complementary part of the total variation is called unexplained or residual variation.

Definition in terms of information gain

Information gain by better modelling

Following Kent (1983),[1] we use the Fraser information (Fraser 1965)[2]

[math]\displaystyle{ F(\theta) = \int \textrm{d}r\,g(r)\,\ln f(r;\theta) }[/math]

where [math]\displaystyle{ g(r) }[/math] is the probability density of a random variable [math]\displaystyle{ R\, }[/math], and [math]\displaystyle{ f(r;\theta)\, }[/math] with [math]\displaystyle{ \theta\in\Theta_i }[/math] ([math]\displaystyle{ i=0,1\, }[/math]) are two families of parametric models. Model family 0 is the simpler one, with a restricted parameter space [math]\displaystyle{ \Theta_0\subset\Theta_1 }[/math].

Parameters are determined by maximum likelihood estimation,

[math]\displaystyle{ \theta_i = \operatorname{arg max}_{\theta\in\Theta_i} F(\theta). }[/math]

The information gain of model 1 over model 0 is written as

[math]\displaystyle{ \Gamma(\theta_1:\theta_0) = 2 [ F(\theta_1)-F(\theta_0) ]\, }[/math]

where a factor of 2 is included for convenience. Γ is always nonnegative; it measures the extent to which the best model of family 1 is better than the best model of family 0 in explaining g(r).

Information gain by a conditional model

Assume a two-dimensional random variable [math]\displaystyle{ R=(X,Y) }[/math] where X shall be considered as an explanatory variable, and Y as a dependent variable. Models of family 1 "explain" Y in terms of X,

[math]\displaystyle{ f(y\mid x;\theta) }[/math],

whereas in family 0, X and Y are assumed to be independent. We define the randomness of Y by [math]\displaystyle{ D(Y)=\exp[-2F(\theta_0)] }[/math], and the randomness of Y, given X, by [math]\displaystyle{ D(Y\mid X)=\exp[-2F(\theta_1)] }[/math]. Then,

[math]\displaystyle{ \rho_C^2 = 1-D(Y\mid X)/D(Y) }[/math]

can be interpreted as proportion of the data dispersion which is "explained" by X.

Special cases and generalized usage

Linear regression

Main page: Fraction of variance unexplained

The fraction of variance unexplained is an established concept in the context of linear regression. The usual definition of the coefficient of determination is based on the fundamental concept of explained variance.

Correlation coefficient as measure of explained variance

Let X be a random vector, and Y a random variable that is modeled by a normal distribution with centre [math]\displaystyle{ \mu=\Psi^\textrm{T}X }[/math]. In this case, the above-derived proportion of explained variation [math]\displaystyle{ \rho_C^2 }[/math] equals the squared correlation coefficient [math]\displaystyle{ R^2 }[/math].

Note the strong model assumptions: the centre of the Y distribution must be a linear function of X, and for any given x, the Y distribution must be normal. In other situations, it is generally not justified to interpret [math]\displaystyle{ R^2 }[/math] as proportion of explained variance.

In principal component analysis

Explained variance is routinely used in principal component analysis. The relation to the Fraser–Kent information gain remains to be clarified.

Criticism

As the fraction of "explained variance" equals the squared correlation coefficient [math]\displaystyle{ R^2 }[/math], it shares all the disadvantages of the latter: it reflects not only the quality of the regression, but also the distribution of the independent (conditioning) variables.

In the words of one critic: "Thus [math]\displaystyle{ R^2 }[/math] gives the 'percentage of variance explained' by the regression, an expression that, for most social scientists, is of doubtful meaning but great rhetorical value. If this number is large, the regression gives a good fit, and there is little point in searching for additional variables. Other regression equations on different data sets are said to be less satisfactory or less powerful if their [math]\displaystyle{ R^2 }[/math] is lower. Nothing about [math]\displaystyle{ R^2 }[/math] supports these claims".[3]:58 And, after constructing an example where [math]\displaystyle{ R^2 }[/math] is enhanced just by jointly considering data from two different populations: "'Explained variance' explains nothing."[3][page needed][4]:183

See also

References

  1. Kent, J. T. (1983). "Information gain and a general measure of correlation". Biometrika 70 (1): 163–173. doi:10.1093/biomet/70.1.163. 
  2. Fraser, D. A. S. (1965). "On Information in Statistics". Ann. Math. Statist. 36 (3): 890–896. doi:10.1214/aoms/1177700061. 
  3. 3.0 3.1 Achen, C. H. (1982). Interpreting and Using Regression. Beverly Hills: Sage. pp. 58–59. ISBN 0-8039-1915-8. 
  4. Achen, C. H. (1990). "'What Does "Explained Variance" Explain?: Reply". Political Analysis 2 (1): 173–184. doi:10.1093/pan/2.1.173. 

External links