Chapter Contents 
Previous 
Next 
The ARIMA Procedure 
The autocorrelation plots for this series, as shown in the previous section, suggest an AR(1) model for the change in SALES. You should check the diagnostic statistics to see if the AR(1) model is adequate. Other candidate models include an MA(1) model, and loworder mixed ARMA models. In this example, the AR(1) model is tried first.
estimate p=1; run;
The ESTIMATE statement fits the model to the data and prints parameter estimates and various diagnostic statistics that indicate how well the model fits the data. The first part of the ESTIMATE statement output, the table of parameter estimates, is shown in Figure 7.10.

The table of parameter estimates is titled "Condidtional Least Squares Estimation," which indicates the estimation method used. You can request different estimation methods with the METHOD= option.
The table of parameter estimates lists the parameters in the model; for each parameter, the table shows the estimated value and the standard error and t value for the estimate. The table also indicates the lag at which the parameter appears in the model.
In this case, there are two parameters in the model. The mean term is labeled MU; its estimated value is .90280. The autoregressive parameter is labeled AR1,1; this is the coefficient of the lagged value of the change in SALES, and its estimate is .86847.
The t values provide significance tests for the parameter estimates and indicate whether some terms in the model may be unnecessary. In this case, the t value for the autoregressive parameter is 15.83, so this term is highly significant. The t value for MU indicates that the mean term adds little to the model.
The standard error estimates are based on large sample theory. Thus, the standard errors are labeled as approximate, and the standard errors and t values may not be reliable in small samples.
The next part of the ESTIMATE statement output is a table of goodnessoffit statistics, which aid in comparing this model to other models. This output is shown in Figure 7.11.

The "Constant Estimate" is a function of the mean term MU and the autoregressive parameters. This estimate is computed only for AR or ARMA models, but not for strictly MA models. See the section "General Notation for ARIMA Models" later in this chapter for an explanation of the constant estimate.
The "Variance Estimate" is the variance of the residual series, which estimates the innovation variance. The item labeled "Std Error Estimate" is the square root of the variance estimate. In general, when comparing candidate models, smaller AIC and SBC statistics indicate the better fitting model. The section "Estimation Details" later in this chapter explains the AIC and SBC statistics.
The ESTIMATE statement next prints a table of correlations of the parameter estimates, as shown in Figure 7.12. This table can help you assess the extent to which collinearity may have influenced the results. If two parameter estimates are very highly correlated, you might consider dropping one of them from the model.

The next part of the ESTIMATE statement output is a check of the autocorrelations of the residuals. This output has the same form as the autocorrelation check for white noise that the IDENTIFY statement prints for the response series. The autocorrelation check of residuals is shown in Figure 7.13.
The test statistics for the residuals series indicate whether the residuals are uncorrelated (white noise) or contain additional information that might be utilized by a more complex model. In this case, the test statistics reject the noautocorrelation hypothesis at a high level of significance. (p=0.0019 for the first six lags.) This means that the residuals are not white noise, and so the AR(1) model is not a fully adequate model for this series.
The final part of the ESTIMATE statement output is a listing of the estimated model using the back shift notation. This output is shown in Figure 7.14.

This listing combines the differencing specification given in the IDENTIFY statement with the parameter estimates of the model for the change in sales. Since the AR(1) model is for the change in sales, the final model for sales is an ARIMA(1,1,0) model. Using B, the back shift operator, the mathematical form of the estimated model shown in this output is as follows:
See the section "General Notation for ARIMA Model" later in this chapter for further explanation of this notation.
An ARMA(1,1) model predicts the change in SALES as an average change, plus some fraction of the previous change, plus a random error, plus some fraction of the random error in the preceding period. An ARMA(1,1) model for the change in sales is the same as an ARIMA(1,1,1) model for the level of sales.
To estimate a mixed autoregressive moving average model, you specify the order of the moving average part of the model with the Q= option on an ESTIMATE statement in addition to specifying the order of the autoregressive part with the P= option. The following statements fit an ARMA(1,1) model to the differenced SALES series:
estimate p=1 q=1; run;
The parameter estimates table and goodnessoffit statistics for this model are shown in Figure 7.15.

The moving average parameter estimate, labeled "MA1,1", is 0.58935. Both the moving average and the autoregressive parameters have significant t values. Note that the variance estimate, AIC, and SBC are all smaller than they were for the AR(1) model, indicating that the ARMA(1,1) model fits the data better without overparameterizing.
The check for white noise residuals is shown in Figure 7.16. The tests show that we cannot reject the hypothesis that the residuals are uncorrelated. Thus, you conclude that the ARMA(1,1) model is adequate for the change in sales series, and there is no point in trying more complex models.

The output showing the form of the estimated ARIMA(1,1,1) model for SALES is shown in Figure 7.17.

The estimated model shown in this output is
Since the model diagnostic tests show that all the parameter estimates are significant and the residual series is white noise, the estimation and diagnostic checking stage is complete. You can now proceed to forecasting the SALES series with this ARIMA(1,1,1) model.
Chapter Contents 
Previous 
Next 
Top 
Copyright © 1999 by SAS Institute Inc., Cary, NC, USA. All rights reserved.