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 = '周期成分:時間変動あり')