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, xlimmultiplier=1.96).xyplot.seqplot