瞭解時間序列對象R（頻譜分析）的頻率

Question

我想知道我的數據中是否有6個月的季節性。我想用光譜分析來做到這一點，我會在特定的頻率下觀察頻譜的Ftest。

我尋找有助於瞭解什麼x軸在譜分析裝置（即什麼頻率的意思），如果我計算它基於與頻率12.

這裏一些數據的時間序列的對象上產生。每個月我的數據被收集超過10年：

data = ts(1:120, frequency = 12, start = c(1959, 1)) # data that is collected each month over 10 years. 
library(multitaper) 
resSpec <- spec.mtm(data, nFFT = "default", 
       Ftest = TRUE, 
       jackknife = FALSE, maxAdaptiveIterations = 100, 
       plot = TRUE, na.action = na.fail) 

Output from multitaper

如果我在6個月的尋找光譜功率（每年2次）循環，我在尋找頻率1/6或6/12？

Answer 1

頻率輸出是每年週期，而不是一個月。您可以檢查此爲這樣：

curve <- function(i) { 
    sin(2*pi/6 * i) 
} 

data = ts(sapply(1:120, curve), frequency = 12, start = c(1959, 1)) 

resSpec <- spec.mtm(data, nFFT = "default", 
       Ftest = TRUE, 
       jackknife = FALSE, maxAdaptiveIterations = 100, 
       plot = TRUE, na.action = na.fail, dtUnits="year", deltat=1/12) 

所以，你要尋找的頻率爲2（次/年）