Here is my function ploting the temporal data and the forecast of a given time series:
Function[{dates, trees, keys, diff},
Manipulate[Module[{ts, tsm},
ts = TimeSeries[Transpose@{dates, trees[key]}];
tsm = TimeSeriesModelFit[ts, "SARMA"];
DateListPlot[{tsm["TemporalData"], TimeSeriesForecast[tsm, {90}]},
PlotLabel -> tsm,
PlotLegends -> {"TemporalData", "Forecast"},
ImageSize -> imageSize]]
, {key, keys}]
];
I want to show the 50% and 95% confidence interval as bands on the forecasting part, is there a way to modify my function to show it?
trees
is just the name of my dataset.