TimeSeriesAnalysis
Forecast
forecast a time series using an exponential smoothing model
Calling Sequence
Parameters
Description
Options
Examples
References
Compatibility
Forecast(model, n, extraparameters, opts)
Forecast(model, ts, n, extraparameters, opts)
model
-
Exponential smoothing model
ts
Time series consisting of a single data set
n
nonnegative integer indicating how many data points should be forecast
extraparameters
(optional) table of parameter values
opts
(optional) equation(s) of the form optionname = value, where optionname is one of output, runs, and evolve
The Forecast command generates a forecast for a time series using an exponential smoothing model.
If a time series is not specified in the calling sequence, then the time series that model is created from will be used.
The number n is the number of data points that will be forecast.
The value of any of the parameters can be overridden by supplying the extraparameters option. It is a table in the format returned by Initialize: its indices are global names corresponding to parameters, and the values are their values as explained in the Exponential smoothing model help page.
output= point, runs, percentiles(p1, p2, ..., pn), quantiles(q1, q2, ..., qn), or confidenceintervals(c1, c2, ..., cn)
By default, Maple will generate a point forecast by assuming the future residuals will all be 0. Other forms of output can be selected using the output option.
output = percentiles(p1, p2, ..., pn) returns n time series that represent the p1th, p2th, ..., pnth percentile of the distribution of the values generated by the model. This is obtained via simulation.
output = quantiles(q1, q2, ..., qn) returns n time series that represent the q1th, q2th, ..., qnth quantile of the distribution of the values generated by the model. The only difference with percentiles is that, e.g., p1 = 5 will be interpreted in the same way as q1 = 0.05.
output = confidenceintervals(c1, c2, ..., cn) returns 2⁢n+1 time series: the 50 - ci/2 and 50 + ci/2 percentiles for i from 1 to n, and the 50th percentile. For example, confidenceintervals(80, 95) will return the 2.5, 10, 50, 90, and 97.5 percentiles.
output = runs returns all simulated future time series. This requires that the runs = k option is also included. (See below.)
output = point selects the default behavior of point forecasts obtained by setting all future residuals to 0.
runs= nonnegative integer
For any setting of the output option other than point, Maple will run multiple forecasts. The runs option determines the number of forecasts used to arrive at this output. If output is of the form percentiles(p1, p2, ..., pn) or quantiles(q1, q2, ..., qn) or confidenceintervals(c1, c2, ..., cn), then by default, Maple chooses the number of runs between 50 and 2000 depending on the particular values chosen. For example, for percentiles(10, 50, 70), the default number of runs is 225. The runs option overrides this default. It is a required option if output = runs is selected: there is no default value. If output = point, then the runs option is ignored.
evolve= true, false, or a nonnegative integer
Normally, the values of the model parameters represent the model at the beginning of the time series, and the values we would like to forecast are after the end of the time series. In this case, the state of the model needs to evolve through the time series in order to obtain the state at the end. This is the default behavior. If you would like to use the model as is, you can supply the option evolve = false. If you want the model to evolve only a limited number of steps, say k, you can supply the option evolve = k. All of this applies only to evolution of the model before the actual forecast is started.
with⁡TimeSeriesAnalysis:
Consider the following time series. It represents international tourist visitor nights in Australia.
ts≔TimeSeries⁡41.7,24.0,32.3,37.3,46.2,29.3,36.5,43.0,48.9,31.2,37.7,40.4,51.2,31.9,41.0,43.8,55.6,33.9,42.1,45.6,59.8,35.2,44.3,47.9,startdate=2005,frequency=quarterly,header=Visitor nights
ts≔Time seriesVisitor nights24 rows of data:2005-Jan-01 - 2010-Oct-01
Fit an exponential smoothing model to it.
esm≔ExponentialSmoothingModel⁡ts
esm≔< an ETS(M,A,M) model >
Now get the point forecast of this model for two years.
pfc≔Forecast⁡esm,8
pfc≔Time seriesVisitor nights (forecast)8 rows of data:2010-Dec-31 - 2012-Sep-30
TimeSeriesPlot⁡ts,pfc,color=Red
Or specify the time series explicitly.
pfc2≔Forecast⁡esm,ts,8
pfc2≔Time seriesVisitor nights (forecast)8 rows of data:2010-Dec-31 - 2012-Sep-30
TimeSeriesPlot⁡ts,pfc2,color=Red
Add confidence intervals.
cis≔Forecast⁡esm,8,output=confidenceintervals⁡80,95
cis≔Time seriesVisitor nights (forecast - 2 percentile), ..., Visitor nights (forecast - 98 percentile)8 rows of data:2010-Dec-31 - 2012-Sep-30
TimeSeriesPlot⁡ts,cis,color=Red..Green
If you use a larger number of runs, the confidence intervals will be slightly more accurate.
cis2≔Forecast⁡esm,8,output=confidenceintervals⁡80,95,runs=5000
cis2≔Time seriesVisitor nights (forecast - 2 percentile), ..., Visitor nights (forecast - 98 percentile)8 rows of data:2010-Dec-31 - 2012-Sep-30
TimeSeriesPlot⁡ts,cis2,color=Red..Green
You can also see how the model would forecast the data of 2009-2010. That means you need to evolve the model for only four years, or 16 time steps, before starting to forecast.
cis3≔Forecast⁡esm,8,evolve=16,output=confidenceintervals⁡80,95
cis3≔Time seriesVisitor nights (forecast - 2 percentile), ..., Visitor nights (forecast - 98 percentile)8 rows of data:2009-Jan-01 - 2010-Oct-01
TimeSeriesPlot⁡ts,thickness=3,cis3,color=Red..Green,transparency=0.25
Hyndman, R.J. and Athanasopoulos, G. (2013) Forecasting: principles and practice. http://otexts.org/fpp/. Accessed on 2013-10-09.
Hyndman, R.J., Koehler, A.B., Ord, J.K., and Snyder, R.D. (2008) Forecasting with Exponential Smoothing: The State Space Approach. Springer Series in Statistics. Springer-Verlag Berlin Heidelberg.
The TimeSeriesAnalysis[Forecast] command was introduced in Maple 18.
For more information on Maple 18 changes, see Updates in Maple 18.
The TimeSeriesAnalysis[Forecast] command was updated in Maple 2015.
See Also
TimeSeriesAnalysis[OneStepForecasts]
Download Help Document