I am using intra-day data that starts at 9:50am and would like to convert it into 20 minute time intervals so the first period would be from 09:50 to 10:09:59 and the second tim
Here is a useful trick that should maybe be more prominent in the xts documentation.
Start with an xts
object
R> set.seed(42) ## fix seed
R> X <- xts(cumsum(rnorm(100))+100, order.by=Sys.time()+cumsum(runif(100)))
R> head(X)
[,1]
2012-10-05 06:42:20.299761 101.371
2012-10-05 06:42:20.816872 100.806
2012-10-05 06:42:21.668803 101.169
2012-10-05 06:42:22.111599 101.802
2012-10-05 06:42:22.269479 102.207
2012-10-05 06:42:22.711804 102.100
Given this irregular series, we want to subset at regular intervals we impose. Here, I create a two-second interval. Any other would work if it is in the same type as the index, here POSIXct
.
R> ind <- seq(start(X) - as.numeric(start(X)-round(start(X))) + 1,
+ end(X), by="2 secs")
R> head(ind)
[1] "2012-10-05 06:42:21 CDT" "2012-10-05 06:42:23 CDT"
[3] "2012-10-05 06:42:25 CDT" "2012-10-05 06:42:27 CDT"
[5] "2012-10-05 06:42:29 CDT" "2012-10-05 06:42:31 CDT"
R>
The trick now is to merge the regular series with the irregular one, call na.locf()
on it to call the last good irregular obs onto the new time grid -- and to then subset at the time grid:
R> na.locf(merge(X, xts(,ind)))[ind]
X
2012-10-05 06:42:21 100.8063
2012-10-05 06:42:23 102.1004
2012-10-05 06:42:25 105.4730
2012-10-05 06:42:27 107.2635
2012-10-05 06:42:29 104.9588
2012-10-05 06:42:31 101.7505
2012-10-05 06:42:33 104.6884
2012-10-05 06:42:35 103.6441
2012-10-05 06:42:37 101.6476
2012-10-05 06:42:39 98.6246
2012-10-05 06:42:41 97.9922
2012-10-05 06:42:43 97.7545
2012-10-05 06:42:45 101.0187
2012-10-05 06:42:47 98.0331
2012-10-05 06:42:49 100.7752
2012-10-05 06:42:51 103.0702
2012-10-05 06:42:53 102.6578
2012-10-05 06:42:55 103.1342
2012-10-05 06:42:57 103.4714
2012-10-05 06:42:59 102.3683
2012-10-05 06:43:01 105.0394
2012-10-05 06:43:03 103.9775
R>
Voila.