data bcg;
input bcg alder cases controls;
total=cases+controls;
cards;
1 7 1 7593
0 7 1 11719
1 6 11 7143
0 6 14 10184
1 5 28 5611
0 5 22 7561
1 4 16 2208
0 4 28 8117
1 3 20 2438
0 3 19 5588
1 2 36 4356
0 2 11 1625
1 1 47 5245
0 1 6 1234
;
run;

proc genmod data=bcg;
class alder bcg; 
model cases/total = alder bcg/ dist=bin type3;
run;

data bcglille;
input bcg alder cases con1000 conmatch;
total1000=cases+con1000; 
totalmatch=cases+conmatch;
cards;
1 7 1 101 3
0 7 1 137 5
1 6 11 91 48
0 6 14 115 52
1 5 28 82 67
0 5 22 101 133
1 4 16 28 46
0 4 28 87 130
1 3 20 25 50
0 3 19 69 106
1 2 36 63 126
0 2 11 21 62
1 1 47 56 174
0 1 6 24 38
;
run;

proc genmod data=bcglille;
class alder bcg; 
model cases/total1000 = alder bcg/ dist=bin type3;
run;

proc genmod data=bcglille;
class alder bcg; 
model cases/totalmatch = alder bcg/ dist=bin type3;
run;