Example of Forecast with Best ARIMA Model for a non-seasonal model

A marketing analyst wants to use an ARIMA model to generate short-term forecasts for sales of a shampoo product. The analyst collects sales data from the previous three years. The analyst previously examined a time series plot and the autocorrelation function (ACF) plot for the series. Both plots suggest 1 as the starting point for the order of non-seasonal differencing. The data do not exhibit a seasonal pattern on a time series plot, so the analyst chooses to begin with a nonseasonal model. The analyst requests forecasts for the next 3 months.

  1. Open the sample data ShampooSales.MTW.
  2. Choose Stat > Time Series > Forecast with Best ARIMA Model.
  3. In Series, enter Sales.
  4. In Differencing order d, select 1.
  5. Deselect Include the constant term in models.
  6. In Number of forecasts, enter 3.
  7. Select OK.

Interpret the results

The model selection table ranks the models from the search in order by AICc. The ARIMA (0, 1, 2) model has the least AICc. The ARIMA results that follow are for the ARIMA (0, 1, 2) model.

The p-values in the parameters table show that the moving average terms are significant at the 0.05 level. The analyst concludes that the coefficients belong in the model. The p-values for the Modified Box-Pierce (Ljung-Box) statistics are all insignificant at the 0.05 level. The ACF of the residuals and the PACF of the residuals are all within the 0.05 limits on their respective plots. The analyst concludes that the model meets the assumption that the residuals are independent. The analyst concludes that examination of the forecasts is reasonable.

* WARNING * Inestimable ARIMA(p, d, q) models that do not include a constant term:
(2, 1, 2)

Method

Criterion for best modelMinimum AICc
Rows used36
Rows unused0

Model Selection

Model (d = 1)LogLikelihoodAICcAICBIC
p = 0, q = 2*-197.052400.878400.103404.769
p = 1, q = 2-196.989403.311401.978408.199
p = 1, q = 0-201.327407.029406.654409.765
p = 2, q = 0-200.239407.251406.477411.143
p = 1, q = 1-200.440407.655406.880411.546
p = 2, q = 1-201.776412.884411.551417.773
p = 0, q = 1-204.584413.542413.167416.278
p = 0, q = 0-213.614429.350429.229430.784
* Best model with minimum AICc.  Output for the best model follows.

Final Estimates of Parameters

TypeCoefSE CoefT-ValueP-Value
MA   11.2570.1329.520.000
MA   2-0.8820.133-6.620.000
Differencing: 1 Regular
Number of observations after differencing: 35

Model Summary

DFSSMSMSDAICcAICBIC
331310173970.213743.34400.878400.103404.769
MS = variance of the white noise series

Modified Box-Pierce (Ljung-Box) Chi-Square Statistic

Lag12243648
Chi-Square15.9027.15**
DF1022**
P-Value0.1030.206**
* WARNING * Inestimable ARIMA(p, d, q) models that do not include a constant term:
(2, 1, 2)

Forecasts from Time Period 36




95% Limits
Time PeriodForecastSE ForecastLowerUpperActual
37563.19363.0096439.669686.717 
38594.91265.0499467.388722.435 
39594.91276.0553445.813744.010