滚动均方根而不是r中的均值

Question

假设我有 1 小时的数据，我可以使用以下代码通过滚动平均值将数据聚合为 2 小时和 3 小时的数据。

library(zoo)
h1<-c(2, 4, 6)
h2<-rollmean(h1,2)
h3<-rollmean(h1,3)

这是我从中得到的：

h2
[1] 3 5

h3
[1] 4

但我想在聚合时使用下面的root mean square 函数而不是mean 函数。

sqrt(mean(x^2))

以下是所需的输出：

h2_rms
[1] 3.162278 5.099020
h3_rms
[1] 4.320494

Answer 1

您可以使用rollapply 应用您自己的功能：

library(zoo)

h1<-c(2, 4, 6)
rms <- function(x) sqrt(mean(x^2))

h2 <- rollapply(h1,2, rms)
h3 <- rollapply(h1,3, rms)

h2
#[1] 3.162278 5.099020

h3
#[1] 4.320494