px <- seq(as.POSIXct('1949-01-01'),
as.POSIXct('1960-12-01'), by = 'month')
d0 <- data.frame(px, y = AirPassengers)
library(DT)
datatable(d0)
library('forecast')
s <- mstl(x = msts(d0$y, seasonal.period = 12), s.window = 'periodic')
autoplot(s, main = '周期成分:時間変動なし')
s <- mstl(x = msts(d0$y, seasonal.period = 12))
autoplot(s, main = '周期成分:時間変動あり')