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个值