* Required macros: 'srs', 'sys', 'est_srs' ;

 data random ;
  do id = 1 to 200 ;
    y = 50 + 4*rannor(0) ;
output ;
  end ;

 data ordered ;
  do id = 1 to 200 ;
    y = id + rannor(0) ;
output ;
  end ;

 data periodic ;
  do id = 1 to 200 ;
    y = sin(id/5) ;
    output ;
  end ;

proc sgplot data = random; scatter x=id y=y ;
proc sgplot data = ordered ; scatter x=id y=y ;
proc sgplot data = periodic; scatter x=id y=y ; run; 

%srs(frame=random,npop=200,n=20,sample=random_srs,seed=5273) ;
%sys(frame=random,npop=200,k=10,sample=random_sys,seed=5273) ;
proc sgplot data = random_srs; scatter x=id y=y ;
proc sgplot data = random_sys ; scatter x=id y=y ; run;
proc means data = random ; var y ;
%est_srs(sample=random_srs,npop=200,response=y,param=mean) ;
%est_srs(sample=random_sys,npop=200,response=y,param=mean) ;

%srs(frame=ordered,npop=200,n=20,sample=ordered_srs,seed=73) ;
%sys(frame=ordered,npop=200,k=10,sample=ordered_sys,seed=73) ;
proc sgplot data = ordered_srs ; scatter x=id y=y ;
proc sgplot data = ordered_sys ; scatter x=id y=y ; run;
proc means data = ordered ; var y ;
%est_srs(sample=ordered_srs,npop=200,response=y,param=mean) ;
%est_srs(sample=ordered_sys,npop=200,response=y,param=mean) ;

%srs(frame=periodic,npop=200,n=6,sample=periodic_srs,seed=590) ;
%sys(frame=periodic,npop=200,k=31,sample=periodic_sys,seed=590) ;
proc sgplot data = periodic_srs ; scatter x=id y=y ;
proc sgplot data = periodic_sys ; scatter x=id y=y ; run;
proc means data = periodic ; var y ;
%est_srs(sample=periodic_srs,npop=200,response=y,param=mean) ;
%est_srs(sample=periodic_sys,npop=200,response=y,param=mean) ;

proc print ; run ;