# Seasonal decomposition
fit = stl(myts, s.window="periodic")
plot(fit,col="purple")

x=fit$time.series;

(season_dec=x[1:12,1])
plot(season_dec,pch=20,col="red")

tendency=x[,2]
plot(tendency,pch=20,col="red")

#is the noise white ?
acf(x[,3])

monthplot(myts);
library(forecast)
seasonplot(myts)