goptions reset=all;
* filename grafout 'c:grafout\template5.ps';
* goptions targetdevice=pslepsf;

proc greplay  nofs; igout=work.gseg;  delete _all_;
  ** clear out the existing graphics catalog ***;


data anno;
xsys='2'; ysys='2';
function = 'move';  size = 5;
do Y = -3 to 3 by .02;
X = 1-.2*(exp(-.5*Y*Y)); output;
function="draw"; end;

data weiner; array Y(5);

do i=1 to 5; Y(i) = 0; end; t=0; output;
 do t = 1 to 5000;
   do i=1 to 5; Y(i) = Y(i) + normal(1827655);
 end; x1=50; x2=250; x3=x2*5; x4=x3*5; x5=x4*5;
output;
end;

Data next; set weiner;   %let n=50;
array Y y1-y5;
 title "n=&n";
do i=1 to 5; Y(i) = Y(i)/sqrt(&n); end; t=t/&n;
if t<&n+1 then output;
proc gplot annotate=anno;
plot (y1-Y5)*t/overlay vaxis = -3 to 3 by .5 haxis = 0 to 1 by .1
hminor=0 vminor=0;
label Y1 = "Y / sqrt of &n";
label t = "t/&n";
symbol1 v=none i=join r=5 width=3;
run;

Data next; set weiner;   %let n=100;
array Y y1-y5;
 title "n=&n";
do i=1 to 5; Y(i) = Y(i)/sqrt(&n); end; t=t/&n;
if t<&n+1 then output;
proc gplot annotate=anno;
plot (y1-Y5)*t/overlay vaxis = -3 to 3 by .5 haxis = 0 to 1 by .1
hminor=0 vminor=0;
label Y1 = "Y / sqrt of &n";
label t = "t/&n";
run;

Data next; set weiner;   %let n=500;
array Y y1-y5;
 title "n=&n";
do i=1 to 5; Y(i) = Y(i)/sqrt(&n); end; t=t/&n;
if t<&n+1 then output;
proc gplot annotate=anno;
plot (y1-Y5)*t/overlay vaxis = -3 to 3 by .5 haxis = 0 to 1 by .1
hminor=0 vminor=0;
label Y1 = "Y / sqrt of &n";
label t = "t/&n";
run;


Data next; set weiner;   %let n=5000;
array Y y1-y5;
 title "n=&n";
do i=1 to 5; Y(i) = Y(i)/sqrt(&n); end; t=t/&n;
if t<&n+1 then output;
proc gplot annotate=anno;
plot (y1-Y5)*t/overlay vaxis = -3 to 3 by .5 haxis = 0 to 1 by .1
hminor=0 vminor=0;
label Y1 = "Y / sqrt of &n";
label t = "t/&n";
run;


goptions targetdevice=pscolor;

  proc greplay tc=tempcat nofs igout=work.gseg;
  tdef six des='Six panels'

  1/ llx =  0   lly = 54
     ulx =  0   uly = 96
     urx = 44   ury = 96
     lrx = 44   lry = 54

   2/llx = 48   lly = 54
     ulx = 48   uly = 96
     urx = 92   ury = 96
     lrx = 92   lry = 54


  3/ llx =  0   lly =  4
     ulx =  0   uly = 46
     urx = 44   ury = 46
     lrx = 44   lry =  4

   4/llx = 48   lly =  4
     ulx = 48   uly = 46
     urx = 92   ury = 46
     lrx = 92   lry =  4


;

   template = six;

  treplay 1:gplot  2:gplot1 3:gplot2 4:gplot3;

run ;