我有电力系统频率的原始数据。86 400个数字
frequency=a$Ist_Frq
plot.ts(frequency, main="System frequency [Hz]", xlab="Time [s]")
见示例:
原始数据
现在,我得确定四分之一小时的时间间隔。
frequency=ts(a$Ist_Frq, start=1, frequency=900)
[quarter-hour time interval][2]
我的问题是:
有没有办法确定标准偏差在每四分之一小时?
谢谢你的回答。
发布于 2016-10-03 07:20:07
对于这个问题,可能有几种解决方案:
#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)
注:最后一个值可能基于少于900个值。
https://stackoverflow.com/questions/39825575
复制相似问题