load q[57,11]=poidat.asc;
        lbl = { AYRES,  LCOTTON,  LYNBLK,  LBLK52,  BKEVENT,  
                PRPCGNP,  LAGLYN,  LAGL52,  LAGBEVT, year,year2 };

q=q[2:49,.];
   create f0 =lynch 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);

library dstat;
#include dstat.ext;
dstatset;
        var = { AYRES,  LCOTTON,  LYNBLK,  LBLK52,  BKEVENT,  
                PRPCGNP,  LAGLYN,  LAGL52,  LAGBEVT, year,year2 };
__miss=1;
call corr("lynch",var);