5

嗨，我希望按时间对一些分钟数据进行子集化。我通常使用xts做类似的事情：

subset.string <- 'T10:00/T13:00' 
xts.min.obj[subset.string]

获取每天上午 10 点到下午 1 点（含）之间的所有行，并将输出作为 xts 格式。但是对于我的目的来说有点慢......例如

j <- xts(rnorm(10e6),Sys.time()-(10e6:1))
system.time(j['T10:00/T16:00'])
   user  system elapsed 
  5.704   0.577  17.115

我知道这data.table是 v 快并且在对大型数据集进行子集化时，所以我想知道是否与fasttime包一起处理快速 POSIXct 创建，是否值得创建一个类似的函数

dt.time.subset <- function(xts.min.obj, subset.string){
  require(data.table)
  require(fasttime)
  x.dt <- data.table(ts=format(index(xts.min.obj),"%Y-%m-%d %H:%M:%S %Z"),
                     coredata(xts.min.obj))
  out <- x.dt[,some.subsetting.operation.using."%between%"]
  xts(out,fastPOSIXct(out[,ts])
}

将 xts.min.obj 转换为 data.table 添加某种字符索引，然后使用 data.table 对相关行进行子集化 使用带有 fasttime 的输出行索引重新创建 xts 输出？或者对于已经高度优化并用 C 编写的东西，这是否是过多的多余操作？

4

1 回答 1

6

如果您可以在 中指定范围UTC，则可以执行以下操作：

j[(.index(j) %% 86400) %between% c(10*3600, 16*3600 + 60)]
# +60 because xts includes that minute; you'll need to offset the times
# appropriately to match with xts unless you live in UTC :)

j <- xts(rnorm(10e6),Sys.time()-(10e6:1))
system.time(j[(.index(j) %% 86400) %between% c(10*3600, 16*3600 + 60)])
#  user  system elapsed 
#  1.17    0.08    1.25 
# likely faster on your machine as mine takes minutes to run the OP bench
于 2013-06-27T16:22:40.893 回答