ap_df <- tsdf(AirPassengers)
# periodic if fixed seasonality; doesn't work well:
ggplot(ap_df, aes(x = x, y = y)) +
stat_stl(frequency = 12, s.window = "periodic")
# seasonality varies a bit over time, works better:
ggplot(ap_df, aes(x = x, y = y)) +
stat_stl(frequency = 12, s.window = 7)
ggplot(ldeaths_df, aes(x = YearMon, y = deaths, colour = sex)) +
geom_point() +
facet_wrap(~sex) +
stat_stl(frequency = 12, s.window = 7) +
ggtitle("Seasonally adjusted lung deaths")Run the code above in your browser using DataLab