w<-read.csv("BOGAMBNS.csv")
x<-ts(w[,2],start=c(1959,1),freq=12) #必须转化为时间序列格式
用stl函数直接做季节分解
bstl<-stl(x,"per")
plot(bstl)
用HoltWinters函数做指数平滑或季节分解(这里趋势和水平结合起来相当于stl中的趋势项)
(b<-HoltWinters(z))
Holt-Winters