data sales2;
    set sales;
    Time = _n_;
run;

proc arima data=sales2;
identify var=sales(1,12) crosscorr=(Time);
estimate q=(1)(12)(12);
run;