Rで周波数解析や成分分解 - 日常メモ今回は、まずこんなデータを考えてみる。 期間は2006/1/4〜2015/11/20。 Rのspec.pgramを使用してピリオドグラムを求める。 spec.pgramはデフォルトでは対数グラフを作成するが、今回は対数化は行わないので、log="no"を指定する。 > spec.pgram(y,log="no") このピリオドグラムから周波数0.1のところにピークがあることがわかる。つまり、原系列は周波数0.1（周期10）のデータであることがわかる。 このように、時系列データに隠されている周期性を解析できたわけだが、このような手法がスペクトル分析だ（理系学部を卒業した人には、常識中の常識ですね・・・）。 今回はスペクトル分析を、Rのパッケージstatsのspec.pgram関数を用いて行ったということ。 ちなみに、既にわかっていると思うが、原系列のデータは次のようにして作成した周期0.1の
