使用rollapply在当前秒之前2秒到当前秒之后3秒的窗口中计算样本方差



这就是我要做的:

时间
00:00 12
00:01 6
00:02 8
00:03 12
00:04 11
00:05 6
00:06 13
00:07 12
00:08 12

我将期望值设为NA, NA, 8.166, 9.466, 7.466, 6.400, NA, NA, NA,而不是像问题中那样。例如,在00:02

var(c(12,6,8,12,11,6))

要获得这些值,只需使用:将现有代码更改为使用一系列正确的整数——此时您正在获得当前行的方差,该行后退两步,该行前进三步,而不是全部介于两者之间。

zoo::rollapply(data$value, list(-2:3), var, fill = NA)

-2:3为您生成整数c(-2, -1, 0, 1, 2, 3)的向量。

最新更新