seqplot
with confidence bands for the forecast region.seqplot.forecast(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
seqplot