我在使用时遇到问题period.apply
函数适用于我的高分辨率时间序列分析案例。
我想以 10 分钟为间隔计算数据的统计数据(不同时期的平均值、标准差等)。计算每小时的平均值工作正常,如中所述这个答案 https://stackoverflow.com/a/16019730/9735690.
它创建一个新的 xts 对象,并计算每列的平均值。如何计算每列的最大值?
这个可重现的示例描述了我的数据的结构:
library(xts)
start <- as.POSIXct("2018-05-18 00:00")
tseq <- seq(from = start, length.out = 1440, by = "10 mins")
Measurings <- data.frame(
Time = tseq,
Temp = sample(10:37,1440, replace = TRUE, set.seed(seed = 10)),
Variable1 = sample(1:200,1440, replace = TRUE, set.seed(seed = 187)),
Variable2 = sample(300:800,1440, replace = TRUE, set.seed(seed = 333))
)
Measurings_xts <- xts(Measurings[,-1], Measurings$Time)
HourEnds <- endpoints(Measurings_xts, "hours")
Measurings_mean <- period.apply(Measurings_xts, HourEnds, mean)
我认为将函数参数更改为很容易mean
to max
, 像这样:
Measurings_max <- period.apply(Measurings_xts, HourEnds, max)
它提供输出,但只有一列包含总体最大值。我需要每列每小时的最大值。一个简单的解决方案将不胜感激。