@ program reads in three variables from lynching dataset:  @
@ {WEST4  WEST1  PRPCGNP}                                  @
@ and bootstraps the correlation between the whitened Ayres@
@ index (WEST4) and the whitened log-lynching series WEST1 @
new;
/*
load x[49,3]=resamp.lyn;
q=x~seqa(1,1,rows(x));
        lbl = {
          west4,west1,prpcgnp,yr
               };


   create f0 =samplyn with ^lbl,0,8;
    if f0 == -1;
        errorlog "Can't open output file";
        end;
    endif;
    if writer(f0,q) /= rows(q);
        errorlog "Disk Full";
        end;
    endif;
    f0 = close(f0);
*/
repls=500;
let stackr[1,1]= . ;
i=0;
do while i < repls;
n_samp = exctsmpl("samplyn","bootlyn",79.6);
library dstat;
#include dstat.ext;
dstatset;
var = {
        west4,west1
               };
__miss=1;

{ cor,vc,cases,nms,des } = corr("bootlyn",var);

stackr=stackr|cor[2,1];
print cor[2,1];
i=i+1;
closeall;
endo;

stackr=trimr(stackr,1,0);