/* Break in later year */
new;
#include utils.gau;
#include utils2a.gau;

proc(3)=getres(h1,h2,h3,h4,who,startt,maxt,breakt,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(h1,h2,startt,maxt,breakt,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(h1,h3,startt,maxt,breakt,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(h1,h4,startt,maxt,breakt,dth1); 
output on; "done in ";; (hsec-botim)/100;; "seconds"; " ";output off;

retp(res_70_80,res_70_90,res_70_00); 

endp;

botim=hsec;
cls;

startt=45; maxt=35;breakt=20;
dth1=zeros(4,1)+0.025;


output file=redo_wk.out, reset;
output on; "the risks considered here are CVD and cancer"; " "; output off;



/* WHITE MALES */

load h1=haz_70wm;
load h2=haz_80wm;
load h3=haz_90wm;
load h4=haz_00wm;
h1=h1[.,1:2]; h2=h2[.,1:2]; h3=h3[.,1:2]; h4=h4[.,1:2]; 
{res1,res2,res3}=getres(h1,h2,h3,h4,"white males",startt,maxt,breakt,dth1);
save res_70_80wm2=res1,res_70_90wm2=res2,res_70_00wm2=res3;

breakt=15;
/* WHITE MALES */

load h1=haz_70wm;
load h2=haz_80wm;
load h3=haz_90wm;
load h4=haz_00wm;
h1=h1[.,1:2]; h2=h2[.,1:2]; h3=h3[.,1:2]; h4=h4[.,1:2]; 
{res_70_80wm1,res_70_90wm1,res_70_00wm1}=getres(h1,h2,h3,h4,"white males",startt,maxt,breakt,dth1);
save res_70_80wm1,res_70_90wm1,res_70_00wm1;