To install click the Add extension button. That's it.

The source code for the WIKI 2 extension is being checked by specialists of the Mozilla Foundation, Google, and Apple. You could also do it yourself at any point in time.

4,5
Kelly Slayton
Congratulations on this excellent venture… what a great idea!
Alexander Grigorievskiy
I use WIKI 2 every day and almost forgot how the original Wikipedia looks like.
Live Statistics
English Articles
Improved in 24 Hours
Languages
Recent
Show all languages
What we do. Every page goes through several hundred of perfecting techniques; in live mode. Quite the same Wikipedia. Just better.
.
Leo
Newton
Brights
Milds

# 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, are uncorrelated and have equal variances, the best linear unbiased estimator (BLUE) of the coefficients is given by the ordinary least squares (OLS) estimator, provided it exists. 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.

• 1/5
Views:
7 704
8 861
2 148
4 819
854
• ✪ The Gauss-Markov theorem
• ✪ Proof Gauss Markov Theorem (Regression - OLS)
• ✪ The Gauss-Markov theorem
• ✪ Gauss Markov Theorem (BSE)
• ✪ ECONOMETRICS I Gauss Markov Assumptions I Part 1

## Statement

Suppose we have in matrix notation,

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

expanding to,

${\displaystyle y_{i}=\sum _{j=1}^{K}\beta _{j}X_{ij}+\varepsilon _{i}\quad \forall i=1,2,\ldots ,n}$

where ${\displaystyle \beta _{j}}$ are non-random but unobservable parameters, ${\displaystyle X_{ij}}$ are non-random and observable (called the "explanatory variables"), ${\displaystyle \varepsilon _{i}}$ are random, and so ${\displaystyle y_{i}}$ are random. The random variables ${\displaystyle \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 ${\displaystyle \beta _{K+1}}$ with a newly introduced last column of X being unity i.e., ${\displaystyle X_{i(K+1)}=1}$ for all ${\displaystyle i}$. Note that though ${\displaystyle y_{i},}$as sample responses, are observable, the following statements and arguments including assumptions, proofs and the others assume under the only condition of knowing ${\displaystyle X_{ij},}$but not ${\displaystyle y_{i}.}$

The Gauss–Markov assumptions concern the set of error random variables, ${\displaystyle \varepsilon _{i}}$:

• They have mean zero: ${\displaystyle \operatorname {E} [\varepsilon _{i}]=0.}$
• They are homoscedastic, that is all have the same finite variance: ${\displaystyle \operatorname {Var} (\varepsilon _{i})=\sigma ^{2}<\infty ,}$ and
• Distinct error terms are uncorrelated: ${\displaystyle {\text{Cov}}(\varepsilon _{i},\varepsilon _{j})=0,\forall i\neq j.}$

A linear estimator of ${\displaystyle \beta _{j}}$ is a linear combination

${\displaystyle {\widehat {\beta }}_{j}=c_{1j}y_{1}+\cdots +c_{nj}y_{n}}$

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

${\displaystyle \operatorname {E} \left[{\widehat {\beta }}_{j}\right]=\beta _{j}}$

regardless of the values of ${\displaystyle X_{ij}}$. Now, let ${\displaystyle \sum \nolimits _{j=1}^{K}\lambda _{j}\beta _{j}}$ be some linear combination of the coefficients. Then the mean squared error of the corresponding estimation is

${\displaystyle \operatorname {E} \left[\left(\sum _{j=1}^{K}\lambda _{j}\left({\widehat {\beta }}_{j}-\beta _{j}\right)\right)^{2}\right],}$

in other words 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 ${\displaystyle \beta }$ of parameters ${\displaystyle \beta _{j}}$ is one with the smallest mean squared error for every vector ${\displaystyle \lambda }$ of linear combination parameters. This is equivalent to the condition that

${\displaystyle \operatorname {Var} \left({\widetilde {\beta }}\right)-\operatorname {Var} \left({\widehat {\beta }}\right)}$

is a positive semi-definite matrix for every other linear unbiased estimator ${\displaystyle {\widetilde {\beta }}}$.

The ordinary least squares estimator (OLS) is the function

${\displaystyle {\widehat {\beta }}=(X'X)^{-1}X'y}$

of ${\displaystyle y}$ and ${\displaystyle X}$ (where ${\displaystyle X'}$ denotes the transpose of ${\displaystyle X}$) that minimizes the sum of squares of residuals (misprediction amounts):

${\displaystyle \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 ${\displaystyle a_{1}y_{1}+\cdots +a_{n}y_{n}}$ whose coefficients do not depend upon the unobservable ${\displaystyle \beta }$ but whose expected value is always zero.

## Proof

Let ${\displaystyle {\tilde {\beta }}=Cy}$ be another linear estimator of ${\displaystyle \beta }$ with ${\displaystyle C=(X'X)^{-1}X'+D}$ where ${\displaystyle D}$ is a ${\displaystyle K\times n}$ non-zero 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 ${\displaystyle {\widehat {\beta }},}$ the OLS estimator. We calculate:

{\displaystyle {\begin{aligned}\operatorname {E} \left[{\tilde {\beta }}\right]&=\operatorname {E} [Cy]\\&=\operatorname {E} \left[\left((X'X)^{-1}X'+D\right)(X\beta +\varepsilon )\right]\\&=\left((X'X)^{-1}X'+D\right)X\beta +\left((X'X)^{-1}X'+D\right)\operatorname {E} [\varepsilon ]\\&=\left((X'X)^{-1}X'+D\right)X\beta &&\operatorname {E} [\varepsilon ]=0\\&=(X'X)^{-1}X'X\beta +DX\beta \\&=(I_{K}+DX)\beta .\\\end{aligned}}}

Therefore, since ${\displaystyle beta}$ is **un**observable, ${\displaystyle {\tilde {\beta }}}$ is unbiased if and only if ${\displaystyle DX=0}$. Then:

{\displaystyle {\begin{aligned}\operatorname {Var} \left({\tilde {\beta }}\right)&=\operatorname {Var} (Cy)\\&=C{\text{ Var}}(y)C'\\&=\sigma ^{2}CC'\\&=\sigma ^{2}\left((X'X)^{-1}X'+D\right)\left(X(X'X)^{-1}+D'\right)\\&=\sigma ^{2}\left((X'X)^{-1}X'X(X'X)^{-1}+(X'X)^{-1}X'D'+DX(X'X)^{-1}+DD'\right)\\&=\sigma ^{2}(X'X)^{-1}+\sigma ^{2}(X'X)^{-1}(DX)'+\sigma ^{2}DX(X'X)^{-1}+\sigma ^{2}DD'\\&=\sigma ^{2}(X'X)^{-1}+\sigma ^{2}DD'&&DX=0\\&=\operatorname {Var} \left({\widehat {\beta }}\right)+\sigma ^{2}DD'&&\sigma ^{2}(X'X)^{-1}=\operatorname {Var} \left({\widehat {\beta }}\right)\end{aligned}}}

Since DD' is a positive semidefinite matrix, ${\displaystyle \operatorname {Var} \left({\tilde {\beta }}\right)}$ exceeds ${\displaystyle \operatorname {Var} \left({\widehat {\beta }}\right)}$ by a positive semidefinite matrix.

## Remarks on the proof

As it has been stated before, the condition of ${\displaystyle \operatorname {Var} \left({\tilde {\beta }}\right)-\operatorname {Var} \left({\widehat {\beta }}\right)}$ is equivalent to the property that the best linear unbiased estimator of ${\displaystyle \ell ^{t}\beta }$ is ${\displaystyle \ell ^{t}{\widehat {\beta }}}$ (best in the sense that it has minimum variance). To see this, let ${\displaystyle \ell ^{t}{\tilde {\beta }}}$ another linear unbiased estimator of ${\displaystyle \ell ^{t}\beta }$.

{\displaystyle {\begin{aligned}\operatorname {Var} \left(\ell ^{t}{\tilde {\beta }}\right)&=\ell ^{t}\operatorname {Var} \left({\tilde {\beta }}\right)\ell \\&=\sigma ^{2}\ell ^{t}(X'X)^{-1}\ell +\ell ^{t}DD^{t}\ell \\&=\operatorname {Var} \left(\ell ^{t}{\widehat {\beta }}\right)+(D^{t}\ell )^{t}(D^{t}\ell )&&\sigma ^{2}\ell ^{t}(X'X)^{-1}\ell =\operatorname {Var} \left(\ell ^{t}{\widehat {\beta }}\right)\\&=\operatorname {Var} \left(\ell ^{t}{\widehat {\beta }}\right)+\|D^{t}\ell \|\\&\geqslant \operatorname {Var} \left(\ell ^{t}{\widehat {\beta }}\right)\end{aligned}}}

Moreover, equality holds if and only if ${\displaystyle D^{t}\ell =0}$. We calculate

{\displaystyle {\begin{aligned}\ell ^{t}{\tilde {\beta }}&=\ell ^{t}\left(((X'X)^{-1}X'+D)Y\right)&&{\text{ from above}}\\&=\ell ^{t}(X'X)^{-1}X'Y+\ell ^{t}DY\\&=\ell ^{t}{\widehat {\beta }}+(D^{t}\ell )^{t}Y\\&=\ell ^{t}{\widehat {\beta }}&&D^{t}\ell =0\end{aligned}}}

This proves that the equality holds if and only if ${\displaystyle \ell ^{t}{\tilde {\beta }}=\ell ^{t}{\widehat {\beta }}}$ which gives the uniqueness of the OLS estimator as a BLUE.

## Generalized least squares estimator

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

## Gauss–Markov theorem as stated in econometrics

In most treatments of OLS, the regressors (parameters of interest) in the design matrix ${\displaystyle \mathbf {X} }$ are assumed to be fixed in repeated samples. This assumption is considered inappropriate for a predominantly nonexperimental science like econometrics.[3] Instead, the assumptions of the Gauss–Markov theorem are stated conditional on ${\displaystyle \mathbf {X} }$.

### Linearity

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 ${\displaystyle y=\beta _{0}+\beta _{1}x^{2},}$ qualifies as linear while ${\displaystyle y=\beta _{0}+\beta _{1}^{2}x}$ can be transformed to be linear by replacing ${\displaystyle \beta _{1}^{2}}$ by another parameter, say ${\displaystyle \gamma }$. An equation with a parameter dependent on an independent variable does not qualify as linear, for example ${\displaystyle y=\beta _{0}+\beta _{1}(x)\cdot x}$, where ${\displaystyle \beta _{1}(x)}$ is a function of ${\displaystyle x}$.

Data transformations are often used to convert an equation into a linear form. For example, the Cobb–Douglas function—often used in economics—is nonlinear:

${\displaystyle Y=AL^{\alpha }K^{1-\alpha }e^{\varepsilon }}$

But it can be expressed in linear form by taking the natural logarithm of both sides:[4]

${\displaystyle \ln Y=\ln A+\alpha \ln L+(1-\alpha )\ln K+\varepsilon =\beta _{0}+\beta _{1}\ln L+\beta _{2}\ln K+\varepsilon }$

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

One should be aware, however, that the parameters that minimize the residuals of the transformed equation not necessarily minimize the residuals of the original equation.

### Strict exogeneity

For all ${\displaystyle n}$ observations, the expectation—conditional on the regressors—of the error term is zero:[5]

${\displaystyle \operatorname {E} [\,\varepsilon _{i}\mid \mathbf {X} ]=\operatorname {E} [\,\varepsilon _{i}\mid \mathbf {x_{1}} ,\dots ,\mathbf {x_{n}} ]=0.}$

where ${\displaystyle \mathbf {x} _{i}={\begin{bmatrix}x_{i1}&x_{i2}&\dots &x_{ik}\end{bmatrix}}^{\mathsf {T}}}$ is the data vector of regressors for the ith observation, and consequently ${\displaystyle \mathbf {X} ={\begin{bmatrix}\mathbf {x_{1}^{\mathsf {T}}} &\mathbf {x_{2}^{\mathsf {T}}} &\dots &\mathbf {x_{n}^{\mathsf {T}}} \end{bmatrix}}^{\mathsf {T}}}$ is the data matrix or design matrix.

Geometrically, this assumption implies that ${\displaystyle \mathbf {x} _{i}}$ and ${\displaystyle \varepsilon _{i}}$ are orthogonal to each other, so that their inner product (i.e., their cross moment) is zero.

${\displaystyle \operatorname {E} [\,\mathbf {x} _{j}\cdot \varepsilon _{i}\,]={\begin{bmatrix}\operatorname {E} [\,{x}_{j1}\cdot \varepsilon _{i}\,]\\\operatorname {E} [\,{x}_{j2}\cdot \varepsilon _{i}\,]\\\vdots \\\operatorname {E} [\,{x}_{jk}\cdot \varepsilon _{i}\,]\end{bmatrix}}=\mathbf {0} \quad {\text{for all }}i,j\in n}$

This assumption is violated if the explanatory variables are stochastic, for instance when they are measured with error, or are endogenous.[6] 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 ${\displaystyle \mathbf {X} }$ must have full column rank.

${\displaystyle \operatorname {rank} (\mathbf {X} )=k}$

Otherwise ${\displaystyle \mathbf {X'X} }$ is not invertible and the OLS estimator cannot be computed.

A violation of this assumption is perfect multicollinearity, i.e. some explanatory variables are linearly dependent. One scenario in which this will occur is called "dummy variable trap," when a base dummy variable is not omitted resulting in perfect correlation between the dummy variables and the constant term.[7]

Multicollinearity (as long as it is not "perfect") can be present resulting in a less efficient, but still unbiased estimate. The estimates will be less precise and highly sensitive to particular sets of data.[8] Multicollinearity can be detected from condition number or the variance inflation factor, among other tests.

### Spherical errors

The outer product of the error vector must be spherical.

${\displaystyle \operatorname {E} [\,{\boldsymbol {\varepsilon }}{\boldsymbol {\varepsilon ^{\mathsf {T}}}}\mid \mathbf {X} ]=\operatorname {Var} [\,{\boldsymbol {\varepsilon }}\mid \mathbf {X} ]={\begin{bmatrix}\sigma ^{2}&0&\dots &0\\0&\sigma ^{2}&\dots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\dots &\sigma ^{2}\end{bmatrix}}=\sigma ^{2}\mathbf {I} \quad {\text{with }}\sigma ^{2}>0}$

This implies the error term has uniform variance (homoscedasticity) and no serial dependence.[9] If this assumption is violated, OLS is still unbiased, but inefficient. The term "spherical errors" will describe the multivariate normal distribution: if ${\displaystyle \operatorname {Var} [\,{\boldsymbol {\varepsilon }}\mid \mathbf {X} ]=\sigma ^{2}\mathbf {I} }$ in the multivariate normal density, then the equation ${\displaystyle f(x)=c}$ is the formula for a ball centered at μ with radius σ in n-dimensional space.[10]

Heteroskedasticity 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. Heteroskedastic 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."[11] 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 spherical errors, the generalized least squares estimator can be shown to be BLUE.[12]

## Notes

1. ^ Aitken, A. C. (1935). "On Least Squares and Linear Combinations of Observations". Proceedings of the Royal Society of Edinburgh. 55: 42–48.
2. ^ Huang, David S. (1970). Regression and Econometric Methods. New York: John Wiley & Sons. pp. 127–147. ISBN 0-471-41754-8.
3. ^ Hayashi, Fumio (2000). Econometrics. Princeton University Press. p. 13. ISBN 0-691-01018-8.
4. ^ Kennedy 2003, p. 110.
5. ^ Hayashi, Fumio (2000). Econometrics. Princeton University Press. p. 7. ISBN 0-691-01018-8.
6. ^ Johnston, John (1972). Econometric Methods (Second ed.). New York: McGraw-Hill. pp. 267–291. ISBN 0-07-032679-7.
7. ^ Wooldridge, Jeffrey (2012). Introductory Econometrics (Fifth international ed.). South-Western. p. 220. ISBN 978-1-111-53439-4.
8. ^ Johnston, John (1972). Econometric Methods (Second ed.). New York: McGraw-Hill. pp. 159–168. ISBN 0-07-032679-7.
9. ^ Hayashi, Fumio (2000). Econometrics. Princeton University Press. p. 10. ISBN 0-691-01018-8.
10. ^ Greene 2012, p. 23-note.
11. ^ Greene 2010, p. 22.
12. ^ Kennedy 2003, p. 135.

## References

Use of BLUE in physics

Basis of this page is in Wikipedia. Text is available under the CC BY-SA 3.0 Unported License. Non-text media are available under their specified licenses. Wikipedia® is a registered trademark of the Wikimedia Foundation, Inc. WIKI 2 is an independent company and has no affiliation with Wikimedia Foundation.