data bissau;
  filename bisfile url "http://www.biostat.ku.dk/~pka/epidata/bissau.txt";
  infile bisfile firstobs=2;
  input id fuptime dead bcg dtp age agem;
run;











/* 2.1.3 - 14 */

proc freq data=bissau;
tables bcg*dead2/relrisk nocol nopercent;
run;

data bissau; set bissau;
dead2=dead;
if dead=0 then dead2=2;
if dtp>0 then dtp2=1;
if dtp=0 then dtp2=2;
nybcg=2-bcg;
dtpany=0; if dtp>0 then dtpany=1;
years=fuptime/364.24;
logyears=log(years);
run;

proc freq data=bissau;
tables nybcg*dead/relrisk nocol nopercent;
run;

/* 2.1.3 - 15 */

proc freq data=bissau;
tables dtp2*dead2/relrisk nocol nopercent;
run;

/* 2.1.3 - 16 */

proc freq data=bissau;
tables dtp2*bcg/chisq nocol nopercent;
run;

/* 2.1.3 - 17 */

proc freq data=bissau;
where bcg=1;
tables dtp2*dead2/relrisk nocol nopercent;
run;

proc freq data=bissau;
where bcg=2;
tables dtp2*dead2/relrisk nocol nopercent;
run;

/* 2.1.3 - 18 */

proc means data=bissau sum;
class bcg;
var dead years;
run;

proc means data=bissau sum;
class dtpany;
var dead years;
run;

proc means data=bissau sum;
where bcg=1;
class dtpany;
var dead years;
run;


proc means data=bissau sum;
where bcg=2;
class dtpany;
var dead years;
run;

data regn; 
Rbcg1=(33+92)/(539.0896112+984.4360861);
Rbcg2=(95+2)/(878.0392049+14.9379530);
RRbcg=Rbcg1/Rbcg2;

Rdtp0=(33+95)/(539.0896112+878.0392049);
Rdtp1=(92+2)/(984.4360861+14.9379530);
RRdtp=Rdtp1/Rdtp0;

Rdtp0bcg1=(33)/(539.0896112);
Rdtp1bcg1=(92)/(984.4360861);
RRdtpbcg1=Rdtp1bcg1/Rdtp0bcg1;

Rdtp0bcg2=(95)/(878.0392049);
Rdtp1bcg2=(2)/(14.9379530);
RRdtpbcg2=Rdtp1bcg2/Rdtp0bcg2;

proc print;
run;

proc genmod data=bissau;
class bcg;
model dead=bcg / dist=poisson offset=logyears;
estimate 'rate ratio' bcg 1 -1/exp;
run; 


proc genmod data=bissau;
class dtpany;
model dead=dtpany / dist=poisson offset=logyears;
estimate 'rate ratio' dtpany -1 1/exp;
run; 


proc genmod data=bissau;
where bcg=1;
class dtpany;
model dead=dtpany / dist=poisson offset=logyears;
estimate 'rate ratio' dtpany -1 1/exp;
run; 


proc genmod data=bissau;
where bcg=2;
class dtpany;
model dead=dtpany / dist=poisson offset=logyears;
estimate 'rate ratio' dtpany -1 1/exp;
run;