seqplot
with confidence bands for the forecast region.
seqplotForecast(xts, forecast, multiplier = 1.96, series = deparse(substitute(observed)), ylim, CI.percent=round((1-2*(1-pnorm(multiplier)))*100,2), main = paste( series, " with forecast + ", CI.percent, "% CI", sep=""), xlab=NULL, ylab=NULL, ...) ## x.at, xlim
multiplier=1.96
).xyplot
.seqplot