Gauss–Markov theorem

Gauss–Markov theorem Not to be confused with Gauss–Markov process. "BLUE" reindirizza qui. For queue management algorithm, see Blue (queue management algorithm). Part of a series on Regression analysis Models Linear regressionSimple regressionPolynomial regressionGeneral linear model Generalized linear modelDiscrete choiceBinomial regressionBinary regressionLogistic regressionMultinomial logistic regressionMixed logitProbitMultinomial probitOrdered logitOrdered probitPoisson Multilevel modelFixed effectsRandom effectsLinear mixed-effects modelNonlinear mixed-effects model Nonlinear regressionNonparametricSemiparametricRobustQuantileIsotonicPrincipal componentsLeast angleLocalSegmented Errors-in-variables Estimation Least squaresLinearNon-linear OrdinaryWeightedGeneralized PartialTotalNon-negativeRidge regressionRegularized Least absolute deviationsIteratively reweightedBayesianBayesian multivariateLeast-squares spectral analysis Background Regression validationMean and predicted responseErrors and residualsGoodness of fitStudentized residualGauss–Markov theorem Mathematics portal vte In statistics, the Gauss–Markov theorem (or simply Gauss theorem for some authors)[1] states that the ordinary least squares (OLS) estimator has the lowest sampling variance within the class of linear unbiased estimators, if the errors in the linear regression model are uncorrelated, have equal variances and expectation value of zero.[2] 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. Vedere, Per esempio, the James–Stein estimator (which also drops linearity), ridge regression, or simply any degenerate estimator.

The theorem was named after Carl Friedrich Gauss and Andrey Markov, although Gauss' work significantly predates Markov's.[3] But while Gauss derived the result under the assumption of independence and normality, Markov reduced the assumptions to the form stated above.[4] A further generalization to non-spherical errors was given by Alexander Aitken.[5] Contenuti 1 Dichiarazione 1.1 Nota 2 Prova 3 Remarks on the proof 4 Generalized least squares estimator 5 Gauss–Markov theorem as stated in econometrics 5.1 Linearity 5.2 Strict exogeneity 5.3 Full rank 5.4 Spherical errors 6 Guarda anche 6.1 Other unbiased statistics 7 Riferimenti 8 Ulteriori letture 9 External links Statement Suppose we have in matrix notation, {stile di visualizzazione {sottolineare {y}}=X{sottolineare {beta }}+{sottolineare {varepsilon }},quad ({sottolineare {y}},{sottolineare {varepsilon }}in matematica bb {R} ^{n},{sottolineare {beta }}in matematica bb {R} ^{K}{testo{ e }}Xin mathbb {R} ^{ntimes K})} expanding to, {stile di visualizzazione y_{io}=somma _{j=1}^{K}beta _{j}X_{ij}+varepsilon _{io}quad forall i=1,2,ldots ,n} dove {displaystyle beta _{j}} are non-random but unobservable parameters, {stile di visualizzazione X_{ij}} are non-random and observable (called the "explanatory variables"), {displaystyle varepsilon _{io}} are random, e così {stile di visualizzazione y_{io}} are random. The random variables {displaystyle varepsilon _{io}} are called the "disturbance", "rumore" 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., {stile di visualizzazione X_{io(K+1)}=1} per tutti {stile di visualizzazione i} . Note that though {stile di visualizzazione y_{io},} as sample responses, are observable, the following statements and arguments including assumptions, proofs and the others assume under the only condition of knowing {stile di visualizzazione X_{ij},} but not {stile di visualizzazione y_{io}.} The Gauss–Markov assumptions concern the set of error random variables, {displaystyle varepsilon _{io}} : They have mean zero: {nome dell'operatore dello stile di visualizzazione {e} [varepsilon _{io}]=0.} They are homoscedastic, that is all have the same finite variance: {nome dell'operatore dello stile di visualizzazione {Var} (varepsilon _{io})=sigma ^{2}0} In terms of vector multiplication, questo significa {stile di visualizzazione {inizio{bmatrice}K_{1}&cdots &k_{p+1}fine{bmatrice}}{inizio{bmatrice}mathbf {v_{1}} \vdots \mathbf {v} _{p+1}fine{bmatrice}}{inizio{bmatrice}mathbf {v_{1}} &cdots &mathbf {v} _{p+1}fine{bmatrice}}{inizio{bmatrice}K_{1}\vdots \k_{p+1}fine{bmatrice}}=mathbf {K} ^{matematica {T}}{matematico {H}}mathbf {K} =lambda mathbf {K} ^{matematica {T}}mathbf {K} >0} dove {displaystyle lambda } is the eigenvalue corresponding to {displaystyle mathbf {K} } . Inoltre, {displaystyle mathbf {K} ^{matematica {T}}mathbf {K} =somma _{io=1}^{p+1}K_{io}^{2}>0implies lambda >0} Infine, as eigenvector {displaystyle mathbf {K} } was arbitrary, it means all eigenvalues of {stile di visualizzazione {matematico {H}}} are positive, dunque {stile di visualizzazione {matematico {H}}} is positive definite. così, {stile di visualizzazione {simbolo audace {beta }}= sinistra(X^{matematica {T}}Esatto)^{-1}X^{matematica {T}}Y} is indeed a global minimum.

Proof Let {stile di visualizzazione {tilde {beta }}=Cy} be another linear estimator of {displaystyle beta } insieme a {stile di visualizzazione C=(X'X)^{-1}X'+D} dove {stile di visualizzazione D} è un {displaystyle Ktimes 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 {stile di visualizzazione {widehat {beta }},} the OLS estimator. We calculate: {stile di visualizzazione {inizio{allineato}nome operatore {e} sinistra[{tilde {beta }}Giusto]&=operatorname {e} [Cy]\&=operatorname {e} sinistra[sinistra((X'X)^{-1}X'+Dright)(Xbeta +varepsilon )Giusto]\&=left((X'X)^{-1}X'+Dright)Xbeta +left((X'X)^{-1}X'+Dright)nome operatore {e} [varepsilon ]\&=left((X'X)^{-1}X'+Dright)Xbeta &&operatorname {e} [varepsilon ]=0\&=(X'X)^{-1}X'Xbeta +DXbeta \&=(IO_{K}+DX)beta .\end{allineato}}} Perciò, da {displaystyle beta } is unobservable, {stile di visualizzazione {tilde {beta }}} is unbiased if and only if {displaystyle DX=0} . Quindi: {stile di visualizzazione {inizio{allineato}nome operatore {Var} sinistra({tilde {beta }}Giusto)&=operatorname {Var} (Cy)\&=C{testo{ Var}}(y)C'\&=sigma ^{2}CC'\&=sigma ^{2}sinistra((X'X)^{-1}X'+Dright)sinistra(X(X'X)^{-1}+D'right)\&=sigma ^{2}sinistra((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} sinistra({widehat {beta }}Giusto)+sigma ^{2}DD'&&sigma ^{2}(X'X)^{-1}=nome operatore {Var} sinistra({widehat {beta }}Giusto)fine{allineato}}} Since DD' is a positive semidefinite matrix, {nome dell'operatore dello stile di visualizzazione {Var} sinistra({tilde {beta }}Giusto)} exceeds {nome dell'operatore dello stile di visualizzazione {Var} sinistra({widehat {beta }}Giusto)} by a positive semidefinite matrix.

Remarks on the proof As it has been stated before, the condition of {nome dell'operatore dello stile di visualizzazione {Var} sinistra({tilde {beta }}Giusto)-nome operatore {Var} sinistra({widehat {beta }}Giusto)} is a positive semidefinite matrix is equivalent to the property that the best linear unbiased estimator of {displaystyle ell ^{t}beta } è {displaystyle ell ^{t}{widehat {beta }}} (best in the sense that it has minimum variance). Per vedere questo, permettere {displaystyle ell ^{t}{tilde {beta }}} another linear unbiased estimator of {displaystyle ell ^{t}beta } .

{stile di visualizzazione {inizio{allineato}nome operatore {Var} sinistra(ell ^{t}{tilde {beta }}Giusto)&=ell ^{t}nome operatore {Var} sinistra({tilde {beta }}Giusto)ell \&=sigma ^{2}ell ^{t}(X'X)^{-1}ell +ell ^{t}DD^{t}ell \&=operatorname {Var} sinistra(ell ^{t}{widehat {beta }}Giusto)+(D^{t}ell )^{t}(D^{t}ell )&&sigma ^{2}ell ^{t}(X'X)^{-1}ell =operatorname {Var} sinistra(ell ^{t}{widehat {beta }}Giusto)\&=operatorname {Var} sinistra(ell ^{t}{widehat {beta }}Giusto)+|D^{t}ell |\&geq operatorname {Var} sinistra(ell ^{t}{widehat {beta }}Giusto)fine{allineato}}} Inoltre, equality holds if and only if {stile di visualizzazione D^{t}tutto =0} . We calculate {stile di visualizzazione {inizio{allineato}ell ^{t}{tilde {beta }}&=ell ^{t}sinistra(((X'X)^{-1}X'+D)Yright)&&{testo{ 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 =0end{allineato}}} 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,[5] extends the Gauss–Markov theorem to the case where the error vector has a non-scalar covariance matrix.[6] 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.[7] Invece, 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, dire {gamma di stili di visualizzazione } . An equation with a parameter dependent on an independent variable does not qualify as linear, Per esempio {displaystyle y=beta _{0}+beta _{1}(X)cdot x} , dove {displaystyle beta _{1}(X)} is a function of {stile di visualizzazione x} .

Data transformations are often used to convert an equation into a linear form. Per esempio, the Cobb–Douglas function—often used in economics—is nonlinear: {displaystyle Y=AL^{alfa }K^{1-alfa }e^{varepsilon }} But it can be expressed in linear form by taking the natural logarithm of both sides:[8] {displaystyle ln Y=ln A+alpha ln L+(1-alfa )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, però, that the parameters that minimize the residuals of the transformed equation do not necessarily minimize the residuals of the original equation.

Strict exogeneity For all {stile di visualizzazione n} observations, the expectation—conditional on the regressors—of the error term is zero:[9] {nome dell'operatore dello stile di visualizzazione {e} [,varepsilon _{io}mid mathbf {X} ]=nome operatore {e} [,varepsilon _{io}mid mathbf {X} _{1},punti ,mathbf {X} _{n}]=0.} dove {displaystyle mathbf {X} _{io}={inizio{bmatrice}X_{i1}&x_{i2}&cdots &x_{ik}fine{bmatrice}}^{matematica {T}}} is the data vector of regressors for the ith observation, and consequently {displaystyle mathbf {X} ={inizio{bmatrice}mathbf {X} _{1}^{matematica {T}}&mathbf {X} _{2}^{matematica {T}}&cdots &mathbf {X} _{n}^{matematica {T}}fine{bmatrice}}^{matematica {T}}} is the data matrix or design matrix.

Geometricamente, this assumption implies that {displaystyle mathbf {X} _{io}} e {displaystyle varepsilon _{io}} are orthogonal to each other, so that their inner product (cioè., their cross moment) è zero.

{nome dell'operatore dello stile di visualizzazione {e} [,mathbf {X} _{j}cdot varepsilon _{io},]={inizio{bmatrice}nome operatore {e} [,{X}_{j1}cdot varepsilon _{io},]\nome operatore {e} [,{X}_{j2}cdot varepsilon _{io},]\vdots \operatorname {e} [,{X}_{jk}cdot varepsilon _{io},]fine{bmatrice}}=mathbf {0} quad {testo{per tutti }}io,jin n} This assumption is violated if the explanatory variables are stochastic, for instance when they are measured with error, or are endogenous.[10] 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.

{nome dell'operatore dello stile di visualizzazione {rango} (mathbf {X} )= k} Altrimenti {displaystyle mathbf {X} 'mathbf {X} } is not invertible and the OLS estimator cannot be computed.

A violation of this assumption is perfect multicollinearity, cioè. 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.[11] 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.[12] 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.

{nome dell'operatore dello stile di visualizzazione {e} [,{simbolo audace {varepsilon }}{simbolo audace {varepsilon ^{matematica {T}}}}mid mathbf {X} ]=nome operatore {Var} [,{simbolo audace {varepsilon }}mid mathbf {X} ]={inizio{bmatrice}sigma ^{2}&0&cdots &0\0&sigma ^{2}&cdots &0\vdots &vdots &ddots &vdots \0&0&cdots &sigma ^{2}fine{bmatrice}}=sigma ^{2}mathbf {io} quad {testo{insieme a }}sigma ^{2}>0} This implies the error term has uniform variance (homoscedasticity) and no serial dependence.[13] If this assumption is violated, OLS is still unbiased, but inefficient. The term "spherical errors" will describe the multivariate normal distribution: Se {nome dell'operatore dello stile di visualizzazione {Var} [,{simbolo audace {varepsilon }}mid mathbf {X} ]=sigma ^{2}mathbf {io} } in the multivariate normal density, poi l'equazione {stile di visualizzazione f(varepsilon )=c} is the formula for a ball centered at μ with radius σ in n-dimensional space.[14] Heteroskedasticity occurs when the amount of error is correlated with an independent variable. Per esempio, 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. Per esempio, 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." 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.[6] See also Independent and identically distributed random variables Linear regression Measurement uncertainty Other unbiased statistics Best linear unbiased prediction (BLUP) Minimum-variance unbiased estimator (MVUE) References ^ See chapter 7 of Johnson, RA.; Wichern, D.W. (2002). Applied multivariate statistical analysis. vol. 5. Prentice hall. ^ Theil, Henri (1971). "Best Linear Unbiased Estimation and Prediction". Principles of Econometrics. New York: John Wiley & Sons. pp. 119–124. ISBN 0-471-85845-5. ^ Plackett, R. l. (1949). "A Historical Note on the Method of Least Squares". Biometrika. 36 (3/4): 458–460. doi:10.2307/2332682. ^ David, F. N.; Neyman, J. (1938). "Extension of the Markoff theorem on least squares". Statistical Research Memoirs. 2: 105–116. OCLC 4025782. ^ Salta su: a b Aitken, UN. C. (1935). "On Least Squares and Linear Combinations of Observations". Proceedings of the Royal Society of Edinburgh. 55: 42–48. doi:10.1017/S0370164600014346. ^ Salta su: a b Huang, David S. (1970). Regression and Econometric Methods. New York: John Wiley & Sons. pp. 127–147. ISBN 0-471-41754-8. ^ Hayashi, Fumio (2000). Econometrics. Stampa dell'Università di Princeton. p. 13. ISBN 0-691-01018-8. ^ Walters, UN. UN. (1970). An Introduction to Econometrics. New York: w. w. Norton. p. 275. ISBN 0-393-09931-8. ^ Hayashi, Fumio (2000). Econometrics. Stampa dell'Università di Princeton. p. 7. ISBN 0-691-01018-8. ^ Johnston, John (1972). Econometric Methods (Seconda ed.). New York: McGraw-Hill. pp. 267–291. ISBN 0-07-032679-7. ^ Wooldridge, Jeffrey (2012). Introductory Econometrics (Fifth international ed.). South-Western. p. 220. ISBN 978-1-111-53439-4. ^ Johnston, John (1972). Econometric Methods (Seconda ed.). New York: McGraw-Hill. pp. 159–168. ISBN 0-07-032679-7. ^ Hayashi, Fumio (2000). Econometrics. Stampa dell'Università di Princeton. p. 10. ISBN 0-691-01018-8. ^ Ramanathan, Ramu (1993). "Nonspherical Disturbances". Statistical Methods in Econometrics. Stampa accademica. pp. 330–351. ISBN 0-12-576830-3. Further reading Davidson, Giacomo (2000). "Statistical Analysis of the Regression Model". Econometric Theory. Oxford: Blackwell. pp. 17–36. ISBN 0-631-17837-6. Goldberger, Artù (1991). "Classical Regression". A Course in Econometrics. Cambridge: Stampa dell'Università di Harvard. pp. 160–169. ISBN 0-674-17544-1. Theil, Henri (1971). "Least Squares and the Standard Linear Model". Principles of Econometrics. New York: John Wiley & Sons. pp. 101–162. ISBN 0-471-85845-5. 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 hide vte Least squares and regression analysis Computational statistics Least squaresLinear least squaresNon-linear least squaresIteratively reweighted least squares Correlation and dependence Pearson product-moment correlationRank correlation (Spearman's rhoKendall's tau)Partial correlationConfounding variable Regression analysis Ordinary least squaresPartial least squaresTotal least squaresRidge regression Regression as a statistical model Linear regression Simple linear regressionOrdinary least squaresGeneralized least squaresWeighted least squaresGeneral linear model Predictor structure Polynomial regressionGrowth curve (statistics)Segmented regressionLocal regression Non-standard Nonlinear regressionNonparametricSemiparametricRobustQuantileIsotonic Non-normal errors Generalized linear modelBinomialPoissonLogistic Decomposition of variance Analysis of varianceAnalysis of covarianceMultivariate AOV Model exploration Stepwise regressionModel selection Mallows's CpAICBICModel specificationRegression validation Background Mean and predicted responseGauss–Markov theoremErrors and residualsGoodness of fitStudentized residualMinimum mean-square errorFrisch–Waugh–Lovell theorem Design of experiments Response surface methodologyOptimal designBayesian design Numerical approximation Numerical analysisApproximation theoryNumerical integrationGaussian quadratureOrthogonal polynomialsChebyshev polynomialsChebyshev nodes Applications Curve fittingCalibration curveNumerical smoothing and differentiationSystem identificationMoving least squares Regression analysis categoryStatistics category Mathematics portalStatistics outlineStatistics topics Categories: Teoremi in statistica

Se vuoi conoscere altri articoli simili a Gauss–Markov theorem puoi visitare la categoria Teoremi in statistica.

lascia un commento

L'indirizzo email non verrà pubblicato.

Vai su

Utilizziamo cookie propri e di terze parti per migliorare l'esperienza dell'utente Maggiori informazioni