求高手帮帮忙

我想通过SAS生成一个AR(1)残差序列图，程序请看附件,中间有许多矛盾的地方，请帮忙修改！

data  userexample5;                                                                                                                     

x1=0.5;                                                                                                                                 

n=-50;                                                                                                                                 

do i=-50 to 1000;                                                                                                                       

e=rannor(123);

xt=0.4*x1+e;

x1=xt;                                                                                                                        

Rt=xt-0.4*x1;

x1=xt;                                                                                                                                                                                                                     

n=n+1;                                                                                                                                 

if i>0 then output;                                                                                                                     

end;                                                                                                                                    

run;                                                                                                                                    

proc print data=userexample5;                                                                                                           

var Rt;                                                                                                                                 

proc gplot data=userexample5;                                                                                                           

symbol i=spline c=blue;                                                                                                                  

plot Rt*n;                                                                                                                              

run;