World Library  
Flag as Inappropriate
Email this Article

Gauss–Markov theorem

Article Id: WHEBN0000170353
Reproduction Date:

Title: Gauss–Markov theorem  
Author: World Heritage Encyclopedia
Language: English
Subject: Regression analysis, Polynomial regression, Errors and residuals, Generalized least squares, Linear least squares (mathematics)
Collection: Statistical Theorems
Publisher: World Heritage Encyclopedia

Gauss–Markov theorem

In statistics, the Gauss–Markov theorem, named after Carl Friedrich Gauss and Andrey Markov, states that in a linear regression model in which the errors have expectation zero and are uncorrelated and have equal variances, the best linear unbiased estimator (BLUE) of the coefficients is given by the ordinary least squares (OLS) estimator. Here "best" means giving the lowest variance of the estimate, as compared to other unbiased, linear estimators. The errors do not need to be normal, nor do they need to be independent and identically distributed (only uncorrelated with mean zero and homoscedastic with finite variance). The requirement that the estimator be unbiased cannot be dropped, since biased estimators exist with lower variance. See, for example, the James–Stein estimator (which also drops linearity) or ridge regression.


  • Statement 1
  • Proof 2
  • Remarks on the proof 3
  • Generalized least squares estimator 4
  • Gauss–Markov theorem as stated in Econometrics 5
    • Linearity 5.1
    • Spherical errors 5.2
    • Exogeneity of independent variables 5.3
    • Full rank 5.4
  • See also 6
    • Other unbiased statistics 6.1
  • Notes 7
  • References 8
  • External links 9


Suppose we have in matrix notation,

\underline{y} = X \underline{\beta} + \underline{\varepsilon},\quad (\underline{y},\underline{\varepsilon} \in \mathbb{R}^n, \beta \in \mathbb{R}^K \text{ and } X\in\mathbb{R}^{n\times K})

expanding to,

y_i=\sum_{j=1}^{K}\beta_j X_{ij}+\varepsilon_i \quad \forall i=1,2,\ldots,n

where \beta_j are non-random but unobservable parameters, X_{ij} are non-random and observable (called the "explanatory variables"), \varepsilon_i are random, and so y_i are random. The random variables \varepsilon_i are called the "disturbance", "noise" or simply "error" (will be contrasted with "residual" later in the article; see errors and residuals in statistics). Note that to include a constant in the model above, one can choose to introduce the constant as a variable \beta_{K+1} with a newly introduced last column of X being unity i.e., X_{i(K+1)} = 1 for all i .

The Gauss–Markov assumptions are

  • E(\varepsilon_i)=0,
  • V(\varepsilon_i)= \sigma^2 < \infty,

(i.e., all disturbances have the same variance; that is "homoscedasticity"), and

  • {\rm cov}(\varepsilon_i,\varepsilon_j) = 0, \forall i \neq j

for i\neq j that is, the error terms are uncorrelated. A linear estimator of \beta_j is a linear combination

\widehat\beta_j = c_{1j}y_1+\cdots+c_{nj}y_n

in which the coefficients c_{ij} are not allowed to depend on the underlying coefficients \beta_j , since those are not observable, but are allowed to depend on the values X_{ij} , since these data are observable. (The dependence of the coefficients on each X_{ij} is typically nonlinear; the estimator is linear in each y_i and hence in each random \varepsilon , which is why this is "linear" regression.) The estimator is said to be unbiased if and only if


regardless of the values of X_{ij} . Now, let \sum_{j=1}^K\lambda_j\beta_j be some linear combination of the coefficients. Then the mean squared error of the corresponding estimation is

E \left(\left(\sum_{j=1}^K\lambda_j(\widehat\beta_j-\beta_j)\right)^2\right);

i.e., it is the expectation of the square of the weighted sum (across parameters) of the differences between the estimators and the corresponding parameters to be estimated. (Since we are considering the case in which all the parameter estimates are unbiased, this mean squared error is the same as the variance of the linear combination.) The best linear unbiased estimator (BLUE) of the vector \beta of parameters \beta_j is one with the smallest mean squared error for every vector \lambda of linear combination parameters. This is equivalent to the condition that

V(\tilde\beta)- V(\widehat\beta)

is a positive semi-definite matrix for every other linear unbiased estimator \tilde\beta.

The ordinary least squares estimator (OLS) is the function


of y and X (where X' denotes the transpose of X ) that minimizes the sum of squares of residuals (misprediction amounts):

\sum_{i=1}^n\left(y_i-\widehat{y}_i\right)^2=\sum_{i=1}^n\left(y_i-\sum_{j=1}^K\widehat\beta_j X_{ij}\right)^2.

The theorem now states that the OLS estimator is a BLUE. The main idea of the proof is that the least-squares estimator is uncorrelated with every linear unbiased estimator of zero, i.e., with every linear combination a_1y_1+\cdots+a_ny_n whose coefficients do not depend upon the unobservable \beta but whose expected value is always zero.


Let \tilde\beta = Cy be another linear estimator of \beta and let C be given by (X'X)^{-1}X' + D , where D is a k \times n nonzero matrix. As we're restricting to unbiased estimators, minimum mean squared error implies minimum variance. The goal is therefore to show that such an estimator has a variance no smaller than that of \hat\beta , the OLS estimator.

The expectation of \tilde\beta is:

\begin{align} E(Cy) &= E(((X'X)^{-1}X' + D)(X\beta + \varepsilon)) \\ &= ((X'X)^{-1}X' + D)X\beta + ((X'X)^{-1}X' + D)\underbrace{E(\varepsilon)}_0 \\ &= (X'X)^{-1}X'X\beta + DX\beta \\ &= (I_k + DX)\beta. \\ \end{align}

Therefore, \tilde\beta is unbiased if and only if DX = 0 .

The variance of \tilde\beta is

\begin{align} V(\tilde\beta) &= V(Cy) = CV(y)C' = \sigma^2 CC' \\ &= \sigma^2((X'X)^{-1}X' + D)(X(X'X)^{-1} + D') \\ &= \sigma^2((X'X)^{-1}X'X(X'X)^{-1} + (X'X)^{-1}X'D' + DX(X'X)^{-1} + DD') \\ &= \sigma^2(X'X)^{-1} + \sigma^2(X'X)^{-1} (\underbrace{DX}_{0})' + \sigma^2 \underbrace{DX}_{0} (X'X)^{-1} + \sigma^2DD' \\ &= \underbrace{\sigma^2(X'X)^{-1}}_{V(\hat\beta)} + \sigma^2DD'. \end{align}

Since DD' is a positive semidefinite matrix, V(\tilde\beta) exceeds V(\hat\beta) by a positive semidefinite matrix.

Remarks on the proof

As it has been stated before, the condition of V(\tilde\beta)- V(\widehat\beta) is equivalent to the property that the best linear unbiased estimator of l^t\beta is l^t\widehat\beta (best in the sense that it has minimum variance). To see this, let l^t\tilde\beta another linear unbiased estimator of l^t\beta .

\begin{align} V(l^t\tilde\beta) &= l^t V(\tilde\beta) l=\underbrace{\sigma^2 l^t (X'X)^{-1}l}_{V(l^t\hat\beta)}+l^tDD^t \\ &= {V(l^t\hat\beta)}+(D^tl)(D^tl)={V(l^t\hat\beta)}+||D^tl||\geq {V(l^t\hat\beta)}\\ \end{align}

Therefore, V(l^t\tilde\beta)\geq V(l^t\hat\beta) .

Moreover, suppose that the equality holds ( V(l^t\tilde\beta)= V(l^t\hat\beta) ). It happens if and only if D^tl=0 . Remembering that, from the proof above, we have \tilde\beta= ((X'X)^{-1}X' + D) Y , then:

\begin{align} l^t\tilde\beta= & l^t(X'X)^{-1}X'Y + l^tDY = & l^t\widehat\beta +\underbrace{(D^tl)^t}_{=0}Y=l^t\widehat\beta \end{align}

This proves that the equality holds if and only if l^t\tilde\beta=l^t\widehat\beta which gives the unicity of the OLS estimator as a BLUE.

Generalized least squares estimator

The generalized least squares (GLS) or Aitken estimator extends the Gauss–Markov theorem to the case where the error vector has a non-scalar covariance matrix – the Aitken estimator is also a BLUE.[1]

Gauss–Markov theorem as stated in Econometrics

In most treatments of OLS, the data *X* is assumed to be fixed. This assumption is considered inappropriate for a predominantly nonexperimental science like econometrics.[2] Instead, the assumptions of the Gauss–Markov theorem are stated conditional on *X*


The dependent variable is assumed to be a linear function of the variables specified in the model. The specification must be linear in its parameters. This does not mean that there must be a linear relationship between the independent and dependent variables. The independent variables can take non-linear forms as long as the parameters are linear. The equation y = \alpha + \beta x^2, \, qualifies as linear while y = \alpha + \beta^2 x can be transformed to be linear by replacing (beta)^2 by another parameter, say gamma. An equation with a parameter dependent on an independent variable does not qualify as linear, for example y = alpha + beta(x) * x, where beta(x) is a function of x.

Data transformations are often used to convert an equation into a linear form (see, however, Santos Silva and Tenreyro, 2006). For example, the Cobb–Douglas function—often used in economics—is nonlinear:

Y=AL^{\alpha}K^{\beta}\varepsilon \,

But it can be expressed in linear form by taking the natural logarithm of both sides:[3] ln Y=ln A + \alpha ln L + \beta lnK + ln\varepsilon

This assumption also covers specification issues: assuming that the proper functional form has been selected and there are no omitted variables.

Spherical errors

\operatorname{Var}[\,\varepsilon|X\,] = \sigma^2 I_n,

Error terms are assumed to be spherical otherwise the OLS estimator is inefficient. The OLS estimator remains unbiased, however. Spherical errors occur when errors have both uniform variance (homoscedasticity) and are uncorrelated with each other.[4] The term "spherical errors" will describe the multivariate normal distribution: if \operatorname{Var}[\,\varepsilon|X\,] = \sigma^2 I_n in the multivariate normal density, then the equation f(x)=c is the formula for a “ball” centered at μ with radius σ in n-dimensional space.[5]

Heteroskedacity occurs when the amount of error is correlated with an independent variable. For example, in a regression on food expenditure and income, the error is correlated with income. Low income people generally spend a similar amount on food, while high income people may spend a very large amount or as little as low income people spend. Heteroskedacity can also be caused by changes in measurement practices. For example, as statistical offices improve their data, measurement error decreases, so the error term declines over time.

This assumption is violated when there is autocorrelation. Autocorrelation can be visualized on a data plot when a given observation is more likely to lie above a fitted line if adjacent observations also lie above the fitted regression line. Autocorrelation is common in time series data where a data series may experience "inertia."[6] If a dependent variable takes a while to fully absorb a shock. Spatial autocorrelation can also occur geographic areas are likely to have similar errors. Autocorrelation may be the result of misspecification such as choosing the wrong functional form. In these cases, correcting the specification is one possible way to deal with autocorrelation.

In the presence of non-spherical errors, the generalized least squares estimator can be shown to be BLUE.[7]

Exogeneity of independent variables

\operatorname{E}[\,\varepsilon|X\,] = 0.

This assumption is violated if the variables are endogenous. Endogeneity can be the result of simultaneity, where causality flows back and forth between both the dependent and independent variable. Instrumental variable techniques are commonly used to address this problem.

Full rank

The sample data matrix must have full rank or OLS cannot be estimated. There must be at least one observation for every parameter being estimated and the data cannot have perfect multicollinearity.[8] Perfect multicollinearity will occur in a "dummy variable trap" when a base dummy variable is not omitted resulting in perfect correlation between the dummy variables and the constant term.

Multicollinearity (as long as it is not "perfect") can be present resulting in a less efficient, but still unbiased estimate.

See also

Other unbiased statistics


  1. ^ A. C. Aitken, "On Least Squares and Linear Combinations of Observations", Proceedings of the Royal Society of Edinburgh, 1935, vol. 55, pp. 42–48.
  2. ^ Hayashi, "Econometrics", 2000, p.13
  3. ^ Kennedy 2003, p. 110.
  4. ^ Kennedy 2003, p. 133.
  5. ^ Greene 2012, p. 23-note.
  6. ^ Greene 2010, p. 22.
  7. ^ Kennedy 2003, p. 135.
  8. ^ Kennedy 2003, p. 205.


  • Greene, William H. (2012, 7th ed.) Econometric Analysis, Prentice Hall.

Use of BLUE in physics

  • L. Lyons, D. Gibaut, P. Clifford (1998). "How to combine correlated estimates of a single physical quantity". Nucl. Instr. and Meth. A270: 110. 
  • L. Lyons, A. J. Martin, D. H. Saxon (1990). "On the determination of the b lifetime by combining the results of different experiments". Phys. Rev. D41: 982–985.  

External links

  • Earliest Known Uses of Some of the Words of Mathematics: G (brief history and explanation of the name)
  • Proof of the Gauss Markov theorem for multiple linear regression (makes use of matrix algebra)
  • A Proof of the Gauss Markov theorem using geometry
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a non-profit organization.

Copyright © World Library Foundation. All rights reserved. eBooks from Project Gutenberg are sponsored by the World Library Foundation,
a 501c(4) Member's Support Non-Profit Organization, and is NOT affiliated with any governmental agency or department.