ods html;

ods graphics on;

 

data simu;

y1=0; y11=0;

do i=-10 to 100;

y=.8*y1-.5*y11+rannor(1234);

if i>0 then output;

y11=y1; y1=y;

end;

proc print data=simu;

run;

data simu1;

set simu;

t=_n_;

run;

proc statespace data=simu1 out=out lead=10;

      var y y1 ; id t;    

run;

 

data out1;

set out;

t=_n_;

run;

 

proc print data=out1;

      id t;

      where t > 50;

   run;

proc gplot data=out1;

      plot for1*t=1 y*t=2 /

           overlay href=200.5;

      symbol1 v=circle i=join;

     

      where t > 50;

   run;

 

 

 

 

ods graphics off;

ods html close;