/* Use the following DATA step to obtain Outputs 10.5 and 10.6                 */

data O_Ring;
 input flt temp td;
datalines;
 1 66 0
 2 70 1
 3 69 0
 4 68 0
 5 67 0
 6 72 0
 7 73 0
 8 70 0
 9 57 1
10 63 1
11 70 1
12 78 0
13 67 0
14 53 1
15 67 0
16 75 0
17 70 0
18 81 0
19 76 0
20 79 0
21 75 1
22 76 0
23 58 1
;

proc freq;
 tables temp*td;
run;

proc genmod data=O_Ring;
 model td=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 not shown in text, but similar to 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;