library(dygraphs) hw <- HoltWinters(ldeaths) p <- predict(hw, n.ahead = 72, prediction.interval = TRUE) p[,'upr'] <- 3000 p[,'lwr'] <- 1000 dygraph(p, main = "Predicted Lung Deaths (UK)") %>% dySeries(c("lwr", "fit", "upr"), label = "Deaths")