R的每個小時的標準偏差

我有電力系統頻率的原始數據。 86 400號碼。R的每個小時的標準偏差

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]

我的問題是：

有什麼辦法如何確定每季度小時非標準偏差？

感謝您的回答。

來源

2016-10-03 Joseph

對不起，這是我的計算器 – Joseph

沒有問題的第一個問題。這就是我們在這裏評論發佈的原因。 – buhtz

可能有幾種解決方案，這一問題：這裏是一個

#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個值

來源

2016-10-03 07:20:07

R的每個小時的標準偏差

回答

相關問題