/* Use the following DATA step to obtain Outputs 10.1 through 10.4   */

data Challenger;
 input temp td no_td;
 total=td+no_td;
datalines;
53 1 0
57 1 0
58 1 0
63 1 0
66 0 1
67 0 3
68 0 1
69 0 1
70 2 2
72 0 1
73 0 1
75 1 1
76 0 2
78 0 1
79 0 1
81 0 1
;

proc genmod data=Challenger;
 model td/total=temp/dist=binomial link=logit type1;
 estimate 'logit at 50 deg' intercept 1 temp 50;
 estimate 'logit at 60 deg' intercept 1 temp 60;
 estimate 'logit at 64.7 deg' intercept 1 temp 64.7;
 estimate 'logit at 64.8 deg' intercept 1 temp 64.8;
 estimate 'logit at 70 deg' intercept 1 temp 70;
 estimate 'logit at 80 deg' intercept 1 temp 80;
 ods output estimates=logit;
run;

/* following creates Output 10.4  */
data prob_hat;
 set logit;
  phat=exp(estimate)/(1+exp(estimate));
  se_phat=phat*(1-phat)*stderr;
  prb_LcL=exp(LowerCL)/(1+exp(LowerCL));
  prb_UcL=exp(UpperCL)/(1+exp(UpperCL));

proc print data=prob_hat;
run;