sspot = ts(dat[,1],freq=12*11) hw = HoltWinters(sspot) p = predict(hw,200) ts.plot(sspot,p,col=c("black","red"))