Durbin–Watson Statistic
   HOME

TheInfoList



OR:

In
statistics Statistics (from German language, German: ', "description of a State (polity), state, a country") is the discipline that concerns the collection, organization, analysis, interpretation, and presentation of data. In applying statistics to a s ...
, the Durbin–Watson statistic is a
test statistic Test statistic is a quantity derived from the sample for statistical hypothesis testing.Berger, R. L.; Casella, G. (2001). ''Statistical Inference'', Duxbury Press, Second Edition (p.374) A hypothesis test is typically specified in terms of a tes ...
used to detect the presence of
autocorrelation Autocorrelation, sometimes known as serial correlation in the discrete time case, measures the correlation of a signal with a delayed copy of itself. Essentially, it quantifies the similarity between observations of a random variable at differe ...
at lag 1 in the residuals (prediction errors) from a regression analysis. It is named after James Durbin and Geoffrey Watson. The small sample distribution of this ratio was derived by
John von Neumann John von Neumann ( ; ; December 28, 1903 – February 8, 1957) was a Hungarian and American mathematician, physicist, computer scientist and engineer. Von Neumann had perhaps the widest coverage of any mathematician of his time, in ...
(von Neumann, 1941). Durbin and Watson (1950, 1951) applied this statistic to the residuals from
least squares The method of least squares is a mathematical optimization technique that aims to determine the best fit function by minimizing the sum of the squares of the differences between the observed values and the predicted values of the model. The me ...
regressions, and developed bounds tests for the
null hypothesis The null hypothesis (often denoted ''H''0) is the claim in scientific research that the effect being studied does not exist. The null hypothesis can also be described as the hypothesis in which no relationship exists between two sets of data o ...
that the errors are serially uncorrelated against the alternative that they follow a first order
autoregressive In statistics, econometrics, and signal processing, an autoregressive (AR) model is a representation of a type of random process; as such, it can be used to describe certain time-varying processes in nature, economics, behavior, etc. The autoregre ...
process. Note that the distribution of this test statistic does not depend on the estimated regression coefficients and the variance of the errors. A similar assessment can be also carried out with the Breusch–Godfrey test and the Ljung–Box test.


Computing and interpreting the Durbin–Watson statistic

If e_t is the residual given by e_t = \rho e_+ \nu_t , the Durbin-Watson
test statistic Test statistic is a quantity derived from the sample for statistical hypothesis testing.Berger, R. L.; Casella, G. (2001). ''Statistical Inference'', Duxbury Press, Second Edition (p.374) A hypothesis test is typically specified in terms of a tes ...
is : d = , where T is the number of observations. For large T, ''d'' is approximately equal to 2(1- \hat), where \hat \rho is the sample autocorrelation of the residuals at lag 1.Gujarati (2003) p. 469 ''d = 2'' therefore indicates no autocorrelation. The value of ''d'' always lies between 0 and 4. If the Durbin–Watson statistic is substantially less than 2, there is evidence of positive serial correlation. As a rough rule of thumb, if Durbin–Watson is less than 1.0, there may be cause for alarm. Small values of ''d'' indicate successive error terms are positively correlated. If ''d>2'', successive error terms are negatively correlated. In regressions, this can imply an underestimation of the level of
statistical significance In statistical hypothesis testing, a result has statistical significance when a result at least as "extreme" would be very infrequent if the null hypothesis were true. More precisely, a study's defined significance level, denoted by \alpha, is the ...
. To test for positive autocorrelation at significance \alpha, the test statistic ''d'' is compared to lower and upper critical values (''d_'' and d_): :*If ''d < d_'', there is statistical evidence that the error terms are positively autocorrelated. :*If ''d > d_'', there is no statistical evidence that the error terms are positively autocorrelated. :*If d_ < d < d_ , the test is inconclusive. Positive serial correlation is serial correlation in which a positive error for one observation increases the chances of a positive error for another observation. To test for negative autocorrelation at significance ''\alpha'', the test statistic (4 - d) is compared to lower and upper critical values (''d_'' and d_): :*If ''(4 - d) < d_'', there is statistical evidence that the error terms are negatively autocorrelated. :*If ''(4 - d) > d_'', there is no statistical evidence that the error terms are negatively autocorrelated. :*If d_ < (4 - d) < d_ , the test is inconclusive. Negative serial correlation implies that a positive error for one observation increases the chance of a negative error for another observation and a negative error for one observation increases the chances of a positive error for another. The critical values, ''d_'' and ''d_'', vary by level of significance (''\alpha'') and the degrees of freedom in the regression equation. Their derivation is complex—statisticians typically obtain them from the appendices of statistical texts. If the
design matrix In statistics and in particular in regression analysis, a design matrix, also known as model matrix or regressor matrix and often denoted by X, is a matrix of values of explanatory variables of a set of objects. Each row represents an individual o ...
\mathbf of the regression is known, exact critical values for the distribution of d under the null hypothesis of no serial correlation can be calculated. Under the null hypothesis d is distributed as : \frac , where n is the number of observations and k is number of regression variables; the \xi_i are independent standard normal random variables; and the \nu_i are the nonzero eigenvalues of ( \mathbf - \mathbf ( \mathbf^T \mathbf ) ^ \mathbf^T ) \mathbf, where \mathbf is the matrix that transforms the residuals into the d statistic, i.e. d = \mathbf^T\mathbf\mathbf. A number of computational algorithms for finding percentiles of this distribution are available. Although serial correlation does not affect the consistency of the estimated regression coefficients, it does affect our ability to conduct valid statistical tests. First, the F-statistic to test for overall significance of the regression may be inflated under positive serial correlation because the mean squared error (MSE) will tend to underestimate the population error variance. Second, positive serial correlation typically causes the ordinary least squares (OLS) standard errors for the regression coefficients to underestimate the true standard errors. As a consequence, if positive serial correlation is present in the regression, standard linear regression analysis will typically lead us to compute artificially small standard errors for the regression coefficient. These small standard errors will cause the estimated t-statistic to be inflated, suggesting significance where perhaps there is none. The inflated t-statistic, may in turn, lead us to incorrectly reject null hypotheses, about population values of the parameters of the regression model more often than we would if the standard errors were correctly estimated. If the Durbin–Watson statistic indicates the presence of serial correlation of the residuals, this can be remedied by using the Cochrane–Orcutt procedure. The Durbin–Watson statistic, while displayed by many regression analysis programs, is not applicable in certain situations. For instance, when lagged dependent variables are included in the explanatory variables, then it is inappropriate to use this test. Durbin's h-test (see below) or likelihood ratio tests, that are valid in large samples, should be used.


Durbin h-statistic

The Durbin–Watson statistic is biased for
autoregressive moving average model In the statistical analysis of time series, autoregressive–moving-average (ARMA) models are a way to describe a (weakly) stationary stochastic process using autoregression (AR) and a moving average (MA), each with a polynomial. They are a too ...
s, so that autocorrelation is underestimated. But for large samples one can easily compute the unbiased
normally distributed In probability theory and statistics, a normal distribution or Gaussian distribution is a type of continuous probability distribution for a real number, real-valued random variable. The general form of its probability density function is f(x ...
h-statistic: :h = \left( 1 - \frac d \right) \sqrt, using the Durbin–Watson statistic ''d'' and the estimated variance : \widehat (\widehat\beta_1) of the regression coefficient of the lagged dependent variable, provided : T \cdot \widehat(\widehat\beta_1)<1. \,


Implementations in statistics packages

# R: the dwtest function in the lmtest package, durbinWatsonTest (or dwt for short) function in the car package, and pdwtest and pbnftest for panel models in the plm package. #
MATLAB MATLAB (an abbreviation of "MATrix LABoratory") is a proprietary multi-paradigm programming language and numeric computing environment developed by MathWorks. MATLAB allows matrix manipulations, plotting of functions and data, implementat ...
: the dwtest function in the Statistics Toolbox. #
Mathematica Wolfram (previously known as Mathematica and Wolfram Mathematica) is a software system with built-in libraries for several areas of technical computing that allows machine learning, statistics, symbolic computation, data manipulation, network ...
: the Durbin–Watson (''d'') statistic is included as an option in the LinearModelFit function. # SAS: Is a standard output when using proc model and is an option (dw) when using proc reg. # EViews: Automatically calculated when using OLS regression #
gretl gretl is an open-source statistical package, mainly for econometrics. The name is an acronym for ''G''nu ''R''egression, ''E''conometrics and ''T''ime-series ''L''ibrary. It has both a graphical user interface (GUI) and a command-line interf ...
: Automatically calculated when using OLS regression #
Stata Stata (, , alternatively , occasionally stylized as STATA) is a general-purpose Statistics, statistical software package developed by StataCorp for data manipulation, visualization, statistics, and automated reporting. It is used by researchers ...
: the command estat dwatson, following regress in time series data. Engle's LM test for autoregressive conditional heteroskedasticity (ARCH), a test for time-dependent volatility, the Breusch–Godfrey test, and Durbin's alternative test for serial correlation are also available. All (except -dwatson-) tests separately for higher-order serial correlations. The Breusch–Godfrey test and Durbin's alternative test also allow regressors that are not strictly exogenous. # Excel: although Microsoft Excel 2007 does not have a specific Durbin–Watson function, the ''d''-statistic may be calculated using =SUMXMY2(x_array,y_array)/SUMSQ(array) # Minitab: the option to report the statistic in the Session window can be found under the "Options" box under Regression and via the "Results" box under General Regression. # Python: a durbin_watson function is included in the statsmodels package (statsmodels.stats.stattools.durbin_watson), but statistical tables for critical values are not available there. #
SPSS SPSS Statistics is a statistical software suite developed by IBM for data management, advanced analytics, multivariate analysis, business intelligence, and criminal investigation. Long produced by SPSS Inc., it was acquired by IBM in 2009. Versi ...
: Included as an option in the Regression function. # Julia: the ''DurbinWatsonTest'' function is available in the ''HypothesisTests'' package.


See also

* Time-series regression * ACF / PACF *
Correlation dimension In chaos theory, the correlation dimension (denoted by ''ν'') is a measure of the dimensionality of the space occupied by a set of random points, often referred to as a type of fractal dimension. For example, if we have a set of random points on t ...
* Breusch–Godfrey test * Ljung–Box test


References


Further reading

* * * * * *


External links


Table for high ''n'' and ''k''
* by Mark Thoma {{DEFAULTSORT:Durbin-Watson statistic Time series statistical tests Autocorrelation