p = predict(fit3,500) ts.plot(ts(sspot),p$pred,p$pred+2*p$se, p$pred-2*p$se, col=c("black","red","blue","blue"))