/*Industry sales for printing paper in thousands of francs*/

ods html;

ods graphics on;

data paper1;

set paper;

time=_n_;

run;

/*proc print data=paper1;

run;*/

proc arima data=paper1;

identify var=y(1,12) esacf scan minic stationarity=(adf=(0));

estimate q=(1)(12) noconstant printall plot;

estimate p=(1,2,3) q=(12) noconstant printall plot;

run;

ods graphics off;

ods html close;