DATA brand;
INPUT GA Male dead nmbrT9 nmbrT10;
GA_24 = GA-24;              /* GA_24:  more "centered" GA */
LINES;
24  0   0    1    1
24  0   1    8    8
24  1   0    2    3
24  1   1   16   15
25  0   0    5    5
25  0   1   20   20
25  1   0    9   11
25  1   1   36   34
26  0   0   14   14
26  0   1   28   28
26  1   0   20   24
26  1   1   40   36
27  0   0   30   30
27  0   1   30   30
27  1   0   20   23
27  1   1   20   17
28  0   0   26   26
28  0   1   13   13
28  1   0   18   20
28  1   1    9    7
29  0   0   20   20
29  0   1    5    5
29  1   0    8    8
29  1   1    2    2
;
RUN;

proc means data = brand mean;
 class male;
 var GA dead;
 FREQ nmbrT9;
RUN;

PROC LOGISTIC DATA = brand DESCENDING ;
  MODEL dead = male  / rl; FREQ nmbrT9;
RUN;

PROC LOGISTIC DATA = brand DESCENDING ;
  MODEL dead = ga_24 / rl; FREQ nmbrT9;
RUN;

PROC LOGISTIC DATA = brand DESCENDING ;
  MODEL dead = ga_24 male / rl ; FREQ nmbrT9;
RUN;



PROC LOGISTIC DATA = brand DESCENDING ;
  MODEL dead = male  / rl; FREQ nmbrT10;
RUN;

PROC LOGISTIC DATA = brand DESCENDING ;
  MODEL dead = ga_24 / rl; FREQ nmbrT10;
RUN;

PROC LOGISTIC DATA = brand DESCENDING ;
  MODEL dead = ga_24 male / rl ; FREQ nmbrT10;
RUN;

proc sort data=brand; by GA Male male;
data b; retain n_alive n_dead n;
 set brand; by GA Male ;
 if first.male then n_alive = nmbrT9;
 if last.male  then do;
   n_dead  = nmbrT9; n = n_dead + n_alive;
   logit2 = log2(n_dead/n_alive);
   logit  =  log(n_dead/n_alive);
   var    =  1/n_dead + 1/n_alive ;
   w      =  1/var ;
   output;
 end;
RUN;
PROC PRINT data=b;RUN;
PROC reg DATA = b  ;
  MODEL logit = ga_24 male  ; WEIGHT w ;
RUN;