R中每四分之一小时的标准偏差

Question

我有电力系统频率的原始数据。 86 400个数字。

 frequency=a$Ist_Frq

 plot.ts(frequency, main="System frequency [Hz]", xlab="Time [s]")

参见示例：

Raw data

现在，我必须确定四分之一小时的时间间隔。

 frequency=ts(a$Ist_Frq, start=1, frequency=900)

[quarter-hour time interval][2]

我的问题是：

有没有办法如何确定每个季度小时的标准偏差？

感谢您的回答。

Answer 1

这个问题可能有几个解决方案：这里有一个

#some data
x <- rnorm(10000)

#identify quarter hour segments 
y <- rep(1:ceiling(length(x)/(15 * 60)), each = 15 * 60)[1:length(x)]

#use tapply to find sd of x for every value of y
tapply(x, y, sd)

nb最后一个值可能基于少于900个值