- data out2;
- set out;
- pi1=lag(x2);
- cs2=lag2(x3);
- fdi1=lag(x4);
- run;
- proc reg data=out2;
- model x1= pi1 cs2 fdi1;
- output out=res r=r;
- run;
- quit;
- symbol1 ci=red l=1 v=none i=join;
- symbol2 ci=black l=13 v=none i=join;
- proc gplot data=res;
- plot r*date=1 /vref=0;
- run;
- quit;
- proc arima data=res;
- identify var=r stationarity=(adf=3) esacf p=(0:7) q=(0:8);
- run;
- quit;
- proc varmax data=out2;
- model x1 pi1 cs2 fdi1 /p=3 cointtest=(johansen=(type=trace));
- run;
- quit;


雷达卡




京公网安备 11010802022788号







