06-08-2016 10:35 AM
Hello, I am new to using ARIMA statistics in SAS. I am using SAS 9.4. I also have access to SAS Enterprise Guide 6.1 I need help reproducing this output. I am not sure what code was used exactly. I am analysing ARIMA models to find the best model for this particular dataset.
The data file (lets call it "a") is a list of 70 years and numbers for each year. Ex:
Here is an example of the part of the code that I have so far:
proc arima data = a;
estimate p=2 q=1;
forecast lead=0 printall;
This produces an ARMA(2,1) model and various statistics.
However, here is the output I would like to replicate:
What is the exact code or required steps I would need to replicate these statistics of fit variables exactly as shown here for each particular ARIMA model? I have tried outstat= and outest= but they do not provide all these variables.
Thanks in advance for your help.
07-22-2016 06:33 AM
You are right. Some of these statistics are not provided by PROC ARIMA out of the box.
You will need to combine the available outputs to calculate those.
Some will be easy: ARSQ can be calculated from k,n,R2 (NPARMS,NUMRESID,SSE - from OUTSTAT=, SST -from the output of IDENTIFY statement)
For some you will need to use the OUT= data set. For example MAXERR needs the RESIDUAL variable in OUT=.
And probably the most interesting one: RWRSQ - you will need to fit a RandomWalk model, and use its results in the calculation.
You could check other SAS/ETS procedures, maybe some can fit the same (or similar) model, but provide the statistics you need.
PROC UCM, PROC VARMAX.