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