proc arima data=aa;
      identify var=x(0,12);
      estimate p=2 q=(2)(12) method=ml;
      forecast lead=8;
   run;

推荐https://bbs.pinggu.org/thread-1038292-1-1.html