[R] adf.test help

Arnaud Battistella arnaudb25 at gmail.com
Wed Feb 17 18:14:06 CET 2010


Hi,

I am trying to test whether a series is return series stationary, but
before proceeding I wanted to make sure I understand correctly how to
use the adf.test function and interpret its output... Could you please
let me know whether I am correct in my interpretations?

ex: I take x such as I know it doesn't have a unit root, and is
therefore stationary

1/
> x <- rnorm(1000)
> adf.test(x, "stationary", k=0)

        Augmented Dickey-Fuller Test

data:  x
Dickey-Fuller = -31.8629, Lag order = 0, p-value = 0.01
alternative hypothesis: stationary

Warning message:
In adf.test(x, "stationary", k = 0) : p-value smaller than printed p-value

If I understand correctly, I am told that the probability of x having
a unit root and therefore being non-stationary is 0.01, so the test
tells me that there is a very high probability that x is stationary.
Then I can conclude that x is mean-reverting. Am I correct?

2/ I would like to see critical values also, so I tried with ur.df

> summary(ur.df(x, "trend", lag=0))

<snip>

Value of test-statistic is: -31.8629 338.4156 507.6231

Critical values for test statistics:
      1pct  5pct 10pct
tau3 -3.96 -3.41 -3.12
phi2  6.09  4.68  4.03
phi3  8.27  6.25  5.34

Here if I understand correctly, as my first critical value is
significantly less than the 1% critical value I reject the null
hypothesis that x has a unit root, so x is stationary and then mean
reverting.

Thanks,

-Arnaud



More information about the R-help mailing list