続 〜 周期変動部分の分解

Rで周波数解析や成分分解」や「RでstatsパッケージのSTLによる結果から周期変動部分、トレンド部分、残差部分を取り出す」で取り扱ったことを違うデータでやってみる。

データは「2つのデータの相関、ccfなど」で使った、fdataとvi。

まずはviについて。

>vi.ts<-ts(as.numeric(vi.ts),frequency=80)
> plot(stl(vi.ts,s.window="per"))

次に、fdataについて。

> fdata.close.ts<-ts(as.numeric(fdata.close.ts),frequency=80)
> plot(stl(fdata.close.ts,s.window="per"))