您根本没有提供太多示例数据来真正让我们了解您想要实现的目标,但这里有一些可能会让您入门。最终,建议归结为将您的数据转换为“xts”,并用于period.apply
在“任意”时间范围内聚合和绘制数据。
这是一些示例数据,首先作为data.frame
从 2012-10-22 12:10:21 EDT 开始每 11 分钟收集一次的随机值,然后转换为xts
:
set.seed(1)
temp <- data.frame(date = seq(as.POSIXct("2012-10-22 12:10:21 EDT", tz="EDT"),
length.out = 100, by = "11 min"),
values = rnorm(100))
library(xts)
x.temp <- xts(temp[-1], temp[, 1])
head(x.temp)
# values
# 2012-10-22 12:10:21 -0.6264538
# 2012-10-22 12:21:21 0.1836433
# 2012-10-22 12:32:21 -0.8356286
# 2012-10-22 12:43:21 1.5952808
# 2012-10-22 12:54:21 0.3295078
# 2012-10-22 13:05:21 -0.8204684
在这里,我们将任意汇总过去每 93 分钟的总和。
my93minAgg <- period.apply(x.temp,
INDEX = endpoints(x.temp, on = "mins", k = 93),
FUN = sum)
my93minAgg
# values
# 2012-10-22 13:16:21 0.3133101
# 2012-10-22 14:55:21 1.1543989
# 2012-10-22 16:23:21 2.1290954
# 2012-10-22 18:02:21 0.5205041
# 2012-10-22 19:30:21 -0.6007850
# 2012-10-22 21:09:21 1.5058904
# 2012-10-22 22:37:21 0.9999036
# 2012-10-23 00:16:21 1.1553448
# 2012-10-23 01:44:21 1.9799783
# 2012-10-23 03:23:21 -1.7152471
# 2012-10-23 04:51:21 2.9882778
# 2012-10-23 06:19:21 0.4580653
这是您绘制它的方式:
plot(my93minAgg, main = "My 93 Minute Aggregation")