時系列オブジェクトの頻度を理解するR（スペクトル解析）

Question

私のデータに6ヶ月の季節性があるかどうかを調べようとしています。スペクトル分析でそれをやりたいのですが、特定の周波数でスペクトルのFtestを見ます。私は助けを探しています

を理解するためにどのようなスペクトル解析手段におけるx軸（すなわち周波数が何を意味するか）私はここで、周波数12

いくつかのデータと時系列オブジェクトに基づいて、それを計算した場合生成される。私は6ヶ月（年2回）サイクルにおけるスペクトルパワーを探していた場合は、

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

私は周波数を探しています1/6：私のデータは、10年間にわたって毎月収集されます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（サイクル/年）である