R-时间从0.1秒到半小时(30分钟)的时间步长，UTC小数小时

Question

我的数据间隔为0.1秒或10行，持续一秒钟因此，基于24*60*60*10。

我想通过将其从0.1第二次步骤汇总到半小时，在我的数据中找到列的平均值（风速和其他变量）。因此，数据将从864000行汇总到48行（一天）

输入：

           tms  Hr  Min Sec Wind speed
7/13/2014 0:00  0   0   0   3.45
7/13/2014 0:00  0   0   0.1 52.34
7/13/2014 0:00  0   0   0.2 1.23
7/13/2014 0:00  0   0   0.3 4.3
7/13/2014 0:00  0   0   0.4 1.34
7/13/2014 0:00  0   0   0.5 3.6

输出我想看到：

Year    Month   Day Hr  Wind speed
7/13/2014   7   13  0   21.92
7/13/2014   7   13  0.5 29.38
7/13/2014   7   13  1   24.18
7/13/2014   7   13  1.5 1.70
7/13/2014   7   13  2   1.80

我的代码以每小时的含义，我想更改为半小时（不是一个小时）的汇总数据。其中dat是没有tms列的数据：因此，我添加了一个日期列。

library(data.table)
library(xts)
dat <- data.table(dat)
tms <- as.POSIXct(seq(0,24*(60*60*10)-1,by=1),origin="2014-07-13",tz="UTC")
xts.ts <- data.frame(xts(dat,tms))

现在，我在我的数据

中添加了tms列

Aut <- data.frame(tms,xts.ts, check.names=FALSE, row.names=NULL) 
mean2 <- aggregate(Aut, 
                   list(hour=cut(as.POSIXct(Aut$tms), "hour")),
                   mean)

但即使每小时也不正确。我想要半小时的数据。有任何建议吗？

Answer 1

正如我在评论中提到的那样，您可以轻松地使用xts::period.apply：

library(xts)
options(digits.secs = 1)  # display fractional seconds
# create 1 day of timestamps that are 0.1 seconds apart
tms <- as.POSIXct(seq(0, 86400-1, by=0.1), origin="2014-07-13", tz="UTC")
# create an xts object with some random data and the times created above
set.seed(21)
xts.ts <- xts(runif(length(tms), 0, 50), tms)
# use period.apply and endpoints to calculate the 30-minute means
mean30min <- period.apply(xts.ts, endpoints(xts.ts, "mins", 30), mean)
# round up to next 30-minute period
mean30min <- align.time(mean30min, 60*30)

如果您希望结果为data.table.table或data.frame，则添加了其他列，则可以在汇总后轻松执行此操作。

library(data.table)
dt.mean30 <- as.data.table(mean30min)
dt.mean30[, Month := .indexmon(mean30min) + 1]
dt.mean30[, Day := .indexmday(mean30min)]
dt.mean30[, Hr := .indexhour(mean30min) + .indexmin(mean30min)/60]