例3.7.sas

来自「第一章 时间序列分析简介 第二章 时间序列的预处理 第三章 平稳时间序列分析 第」· SAS 代码 · 共 20 行

data a;
x_1=0;
e_1=0;
do t=-100 to 1000;
e=rannor(12345);
x=0.5*x_1+e-0.8*e_1;
x_1=x;
e_1=e;
if t>0 then output ;
end;
data a;
set a;
keep t x;
proc arima;
identify var=x nlag=20 outcov=out;
proc gplot data=out;
plot corr*lag partcorr*lag ;
symbol c=red i=needle v=none;
run;