HaltWintersの結果プロット

aはある時系列データ(tsオブジェクト)とする
f:id:seinzumtode:20220118234621p:plain

h <- HoltWinters(a)
plot(h)

f:id:seinzumtode:20220118234720p:plain

原信号、トレンド成分、周期成分、残渣成分に分けてプロットする

decomp <- ts.union(y=h$x,Level=h$fitted[,"level"]+h$fitted[,"trend"],Season=h$fitted[,"season"],Residuals=residuals(h))
plot(decomp)

f:id:seinzumtode:20220118234820p:plain