Applying a statistical test on sliding window

Question

Applying functions on rolling windows of zoo objects is normally quite straightfoward, e.g. a moving average:

z <- zoo(1:10, as.Date(31:40))
rollapply(z, 4, mean, align="right")

Now I want to do the same thing with a statistical test, i.e. apply a Augmented Dickey-Fuller test on each window and get the test statistic like I got the mean in the above example.

So basically I am looking for the equivalent of the following piece of code (which of course doesn't work!):

rollapply(z, 4, ADF.test, align="right")

You are apparently using the uroot package, which is [no longer on CRAN](http://cran.r-project.org/web/packages/uroot/). — Vincent Zoonekynd, Feb 13 '12 at 08:36

score 6 · Accepted Answer · answered Feb 13 '12 at 08:29

6

The following works for me.

library(zoo)
library(tseries)
z <- zoo(rnorm(100), as.Date(1:100))
rollapplyr(z, 20, adf.test)

In case you just want the p-value:

rollapplyr(z, 20, function(u) adf.test(u)$p.value)

answered Feb 13 '12 at 08:29

Vincent Zoonekynd

31,893
5
69
78

Applying a statistical test on sliding window

1 Answers1