[R] Durbin-Watson test in packages "car" and "lmtest"
Torsten.Hothorn at rzmail.uni-erlangen.de
Fri Apr 19 10:49:58 CEST 2002
> P-values in Durbin-Watson test obtained through the use of
> functions available in packages "lmtest" and "car" are different. The
> difference is quite significant. function "dwtest" in "lmtest" is much
> faster than "burbinwatson" in "car". Actually, you can take a nap while
> the latter trying to calculated Durbin-Watson test. My question is which
> p-value is better?
The answer is essencially given in ?durbin.watson and ?dwtest. The latter
The p value is computed
using a Fortran version of the Applied Statistics Algorithm AS 153
by Farebrother (1980, 1984). This algorithm is called "pan" or
"gradsol". For large sample sizes the algorithm might fail to
compute the p value; in that case a warning is printed and an
approximate p value will be given; this p value is computed using
a normal approximation with mean and variance of the Durbin-Watson
while ?durbin.watson says
simulate: if `TRUE' p-values will be estimated by bootstrapping.
What is a "quite significant" difference for p-values?
Looking at the example from ?durbin.watson gives:
R> durbin.watson(lm(fconvict ~ tfr + partic + degrees + mconvict,
lag Autocorrelation D-W Statistic p-value
1 0.688345 0.6168636 0
R> dwtest(fconvict ~ tfr + partic + degrees + mconvict, data=Hartnagel)
data: fconvict ~ tfr + partic + degrees + mconvict
DW = 0.6169, p-value = 6.96e-09
which is fairly close, so you might give us more details (that is: a
working example) to see what the "difference" is (that is: a bug in
either function or a difference due to simulation error / bad
> Thank you,
> Ahmad Abu Hammour
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
More information about the R-help