title: Produce nice linear regression plot (fitted line, confidence / prediction bands, etc)
body: <p>I have this sample 10-year regression in the future.</p>
  
  <pre><code>date&lt;-as.Date(c("2015-12-31", "2014-12-31", "2013-12-31", "2012-12-31"))
  value&lt;-c(16348, 14136, 12733, 10737)
  #fit linear regression
  model&lt;-lm(value~date)
  #build predict dataframe
  dfuture&lt;-data.frame(date=seq(as.Date("2016-12-31"), by="1 year", length.out = 10))
  #predict the futurne
  predict(model, dfuture, interval = "prediction")
  </code></pre>
  
  <p>How can I add confidence bands to this?</p>