Английская Википедия:Cochrane–Orcutt estimation
Cochrane–Orcutt estimation is a procedure in econometrics, which adjusts a linear model for serial correlation in the error term. Developed in the 1940s, it is named after statisticians Donald Cochrane and Guy Orcutt.[1]
Theory
Consider the model
- <math>y_t = \alpha + X_t \beta+\varepsilon_t,\,</math>
where <math>y_{t}</math> is the value of the dependent variable of interest at time t, <math>\beta</math> is a column vector of coefficients to be estimated, <math>X_{t}</math> is a row vector of explanatory variables at time t, and <math>\varepsilon_t</math> is the error term at time t.
If it is found, for instance via the Durbin–Watson statistic, that if the error term is serially correlated over time, then standard statistical inference as normally applied to regressions is invalid because standard errors are estimated with bias. To avoid this problem, the residuals must be modeled. If the process generating the residuals is found to be a stationary first-order autoregressive structure,[2] <math>\varepsilon_t =\rho \varepsilon_{t-1}+e_t,\ |\rho| <1 </math>, with the errors {<math>e_t</math>} being white noise, then the Cochrane–Orcutt procedure can be used to transform the model by taking a quasi-difference:
- <math>y_t - \rho y_{t-1} = \alpha(1-\rho)+(X_t - \rho X_{t-1})\beta + e_t. \,</math>
In this specification the error terms are white noise, so statistical inference is valid. Then the sum of squared residuals (the sum of squared estimates of <math>e_t^2</math>) is minimized with respect to <math>(\alpha,\beta)</math>, conditional on <math>\rho</math>.
Inefficiency
The transformation suggested by Cochrane and Orcutt disregards the first observation of a time series, causing a loss of efficiency that can be substantial in small samples.[3] A superior transformation, which retains the first observation with a weight of <math>\sqrt{(1-\rho^{2})}</math> was first suggested by Prais and Winsten,[4] and later independently by Kadilaya.[5]
Estimating the autoregressive parameter
If <math>\rho</math> is not known, then it is estimated by first regressing the untransformed model and obtaining the residuals {<math>\hat{\varepsilon}_t</math>}, and regressing <math>\hat{\varepsilon}_t</math> on <math>\hat{\varepsilon}_{t-1}</math>, leading to an estimate of <math>\rho</math> and making the transformed regression sketched above feasible. (Note that one data point, the first, is lost in this regression.) This procedure of autoregressing estimated residuals can be done once and the resulting value of <math>\rho</math> can be used in the transformed y regression, or the residuals of the residuals autoregression can themselves be autoregressed in consecutive steps until no substantial change in the estimated value of <math>\rho</math> is observed.
It has to be noted, though, that the iterative Cochrane–Orcutt procedure might converge to a local but not global minimum of the residual sum of squares.[6][7][8] This problem disappears when using the Prais–Winsten transformation instead, which keeps the initial observation.[9]
See also
- Hildreth–Lu estimation
- Newey–West estimator
- Prais–Winsten estimation
- Feasible generalized least squares
References
Further reading
External links
