new;
#include utils.gau;
#include utils1.gau;
#include utils3.gau;

proc(3)=getres(h1,h2,h3,h4,who,startt,maxt,thmin1,thmax1,dth1);
local res_70_80, res_70_90, res_70_00,botim;

botim=hsec;
output on; "comparing "$+who$+" in 1970 to "$+who$+" in 1980"; output off;
res_70_80=grid_seq_uni(h1,h2,startt,maxt,thmin1,thmax1,dth1);
output on; "done in ";; (hsec-botim)/100;; "seconds"; " "; output off;

botim=hsec;
output on; "comparing "$+who$+" in 1970 to "$+who$+" in 1990"; output off;
res_70_90=grid_seq_uni(h1,h3,startt,maxt,thmin1,thmax1,dth1); 
output on; "done in ";; (hsec-botim)/100;; "seconds"; " ";output off;

botim=hsec;
output on; "comparing "$+who$+" in 1970 to "$+who$+" in 2000"; output off;
res_70_00=grid_seq_uni(h1,h4,startt,maxt,thmin1,thmax1,dth1); 
output on; "done in ";; (hsec-botim)/100;; "seconds"; " ";output off;

retp(res_70_80,res_70_90,res_70_00); 

endp;

botim=hsec;

startt=45; maxt=35;
thmin1=0.9; thmax1=1.4; dth1=0.05;


output file=ba1univ.out, reset;



output on; "the risk considered here is CVD"; " "; output off;


/* WHITE MALES */

load h1=haz_70wm;
load h2=haz_80wm;
load h3=haz_90wm;
load h4=haz_00wm;
h1=h1[.,1]; h2=h2[.,1]; h3=h3[.,1]; h4=h4[.,1];
{res_70_80wm_cvd,res_70_90wm_cvd,res_70_00wm_cvd}=getres(h1,h2,h3,h4,"white males",startt,maxt,thmin1,thmax1,dth1);
save res_70_80wm_cvd,res_70_90wm_cvd,res_70_00wm_cvd;

/* WHITE FEMALES */

load h1=haz_70wf;
load h2=haz_80wf;
load h3=haz_90wf;
load h4=haz_00wf;
h1=h1[.,1]; h2=h2[.,1]; h3=h3[.,1]; h4=h4[.,1];
{res_70_80wf_cvd,res_70_90wf_cvd,res_70_00wf_cvd}=getres(h1,h2,h3,h4,"white females",startt,maxt,thmin1,thmax1,dth1);
save res_70_80wf_cvd,res_70_90wf_cvd,res_70_00wf_cvd;

/* BLACK MALES */

load h1=haz_70bm;
load h2=haz_80bm;
load h3=haz_90bm;
load h4=haz_00bm;
h1=h1[.,1]; h2=h2[.,1]; h3=h3[.,1]; h4=h4[.,1];
{res_70_80bm_cvd,res_70_90bm_cvd,res_70_00bm_cvd}=getres(h1,h2,h3,h4,"black males",startt,maxt,thmin1,thmax1,dth1);
save res_70_80bm_cvd,res_70_90bm_cvd,res_70_00bm_cvd;

/* BLACK FEMALES */

load h1=haz_70bf;
load h2=haz_80bf;
load h3=haz_90bf;
load h4=haz_00bf;
h1=h1[.,1]; h2=h2[.,1]; h3=h3[.,1]; h4=h4[.,1];
{res_70_80bf_cvd,res_70_90bf_cvd,res_70_00bf_cvd}=getres(h1,h2,h3,h4,"black females",startt,maxt,thmin1,thmax1,dth1);
save res_70_80bf_cvd,res_70_90bf_cvd,res_70_00bf_cvd;


/* *********************************************************** */



output on; "the risk considered here is Cancer"; " "; output off;


/* WHITE MALES */

load h1=haz_70wm;
load h2=haz_80wm;
load h3=haz_90wm;
load h4=haz_00wm;
h1=h1[.,2]; h2=h2[.,2]; h3=h3[.,2]; h4=h4[.,2];
{res_70_80wm_can,res_70_90wm_can,res_70_00wm_can}=getres(h1,h2,h3,h4,"white males",startt,maxt,thmin1,thmax1,dth1);
save res_70_80wm_can,res_70_90wm_can,res_70_00wm_can;

/* WHITE FEMALES */

load h1=haz_70wf;
load h2=haz_80wf;
load h3=haz_90wf;
load h4=haz_00wf;
h1=h1[.,2]; h2=h2[.,2]; h3=h3[.,2]; h4=h4[.,2];
{res_70_80wf_can,res_70_90wf_can,res_70_00wf_can}=getres(h1,h2,h3,h4,"white females",startt,maxt,thmin1,thmax1,dth1);
save res_70_80wf_can,res_70_90wf_can,res_70_00wf_can;

/* BLACK MALES */

load h1=haz_70bm;
load h2=haz_80bm;
load h3=haz_90bm;
load h4=haz_00bm;
h1=h1[.,2]; h2=h2[.,2]; h3=h3[.,2]; h4=h4[.,2];
{res_70_80bm_can,res_70_90bm_can,res_70_00bm_can}=getres(h1,h2,h3,h4,"black males",startt,maxt,thmin1,thmax1,dth1);
save res_70_80bm_can,res_70_90bm_can,res_70_00bm_can;

/* BLACK FEMALES */

load h1=haz_70bf;
load h2=haz_80bf;
load h3=haz_90bf;
load h4=haz_00bf;
h1=h1[.,2]; h2=h2[.,2]; h3=h3[.,2]; h4=h4[.,2];
{res_70_80bf_can,res_70_90bf_can,res_70_00bf_can}=getres(h1,h2,h3,h4,"black females",startt,maxt,thmin1,thmax1,dth1);
save res_70_80bf_can,res_70_90bf_can,res_70_00bf_can;


/* *********************************************************** */



output on; "the risk considered here is Cancer Excl Lung"; " "; output off;


/* WHITE MALES */

load h1=haz_70wm;
load h2=haz_80wm;
load h3=haz_90wm;
load h4=haz_00wm;
h1=h1[.,2]-h1[.,3]; h2=h2[.,2]-h2[.,3]; h3=h3[.,2]-h3[.,3]; h4=h4[.,2]-h4[.,3];
{res_70_80wm_can_xl,res_70_90wm_can_xl,res_70_00wm_can_xl}=getres(h1,h2,h3,h4,"white males",startt,maxt,thmin1,thmax1,dth1);
save res_70_80wm_can_xl,res_70_90wm_can_xl,res_70_00wm_can_xl;

/* WHITE FEMALES */

load h1=haz_70wf;
load h2=haz_80wf;
load h3=haz_90wf;
load h4=haz_00wf;
h1=h1[.,2]-h1[.,3]; h2=h2[.,2]-h2[.,3]; h3=h3[.,2]-h3[.,3]; h4=h4[.,2]-h4[.,3];
{res_70_80wf_can_xl,res_70_90wf_can_xl,res_70_00wf_can_xl}=getres(h1,h2,h3,h4,"white females",startt,maxt,thmin1,thmax1,dth1);
save res_70_80wf_can_xl,res_70_90wf_can_xl,res_70_00wf_can_xl;

/* BLACK MALES */

load h1=haz_70bm;
load h2=haz_80bm;
load h3=haz_90bm;
load h4=haz_00bm;
h1=h1[.,2]-h1[.,3]; h2=h2[.,2]-h2[.,3]; h3=h3[.,2]-h3[.,3]; h4=h4[.,2]-h4[.,3];
{res_70_80bm_can_xl,res_70_90bm_can_xl,res_70_00bm_can_xl}=getres(h1,h2,h3,h4,"black males",startt,maxt,thmin1,thmax1,dth1);
save res_70_80bm_can_xl,res_70_90bm_can_xl,res_70_00bm_can_xl;

/* BLACK FEMALES */

load h1=haz_70bf;
load h2=haz_80bf;
load h3=haz_90bf;
load h4=haz_00bf;
h1=h1[.,2]-h1[.,3]; h2=h2[.,2]-h2[.,3]; h3=h3[.,2]-h3[.,3]; h4=h4[.,2]-h4[.,3];
{res_70_80bf_can_xl,res_70_90bf_can_xl,res_70_00bf_can_xl}=getres(h1,h2,h3,h4,"black females",startt,maxt,thmin1,thmax1,dth1);
save res_70_80bf_can_xl,res_70_90bf_can_xl,res_70_00bf_can_xl;


/* *********************************************************** */



output on; "the risk considered here is Lung Cancer"; " "; output off;


/* WHITE MALES */

load h1=haz_70wm;
load h2=haz_80wm;
load h3=haz_90wm;
load h4=haz_00wm;
h1=h1[.,3]; h2=h2[.,3]; h3=h3[.,3]; h4=h4[.,3];
{res_70_80wm_lcan,res_70_90wm_lcan,res_70_00wm_lcan}=getres(h1,h2,h3,h4,"white males",startt,maxt,thmin1,thmax1,dth1);
save res_70_80wm_lcan,res_70_90wm_lcan,res_70_00wm_lcan;

/* WHITE FEMALES */

load h1=haz_70wf;
load h2=haz_80wf;
load h3=haz_90wf;
load h4=haz_00wf;
h1=h1[.,3]; h2=h2[.,3]; h3=h3[.,3]; h4=h4[.,3];
{res_70_80wf_lcan,res_70_90wf_lcan,res_70_00wf_lcan}=getres(h1,h2,h3,h4,"white females",startt,maxt,thmin1,thmax1,dth1);
save res_70_80wf_lcan,res_70_90wf_lcan,res_70_00wf_lcan;

/* BLACK MALES */

load h1=haz_70bm;
load h2=haz_80bm;
load h3=haz_90bm;
load h4=haz_00bm;
h1=h1[.,3]; h2=h2[.,3]; h3=h3[.,3]; h4=h4[.,3];
{res_70_80bm_lcan,res_70_90bm_lcan,res_70_00bm_lcan}=getres(h1,h2,h3,h4,"black males",startt,maxt,thmin1,thmax1,dth1);
save res_70_80bm_lcan,res_70_90bm_lcan,res_70_00bm_lcan;

/* BLACK FEMALES */

load h1=haz_70bf;
load h2=haz_80bf;
load h3=haz_90bf;
load h4=haz_00bf;
h1=h1[.,3]; h2=h2[.,3]; h3=h3[.,3]; h4=h4[.,3];
{res_70_80bf_lcan,res_70_90bf_lcan,res_70_00bf_lcan}=getres(h1,h2,h3,h4,"black females",startt,maxt,thmin1,thmax1,dth1);
save res_70_80bf_lcan,res_70_90bf_lcan,res_70_00bf_lcan;


/* *********************************************************** */