了解时间序列对象R中的频率（频谱分析）

Question

我试图找出我的数据中是否有六个月的季节性。 我想通过频谱分析来做到这一点，在这里我将研究特定频率下的频谱Ftest。

如果我基于频率为12的时间序列对象进行计算，我正在寻求帮助以了解光谱分析中的x轴含义（即频率的含义）。

这里生成一些数据。 我的数据在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) 

多锥输出

如果我正在寻找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（周期/年）