/* Use the following SAS statements in conjuction with the DATA step used      */
/*    to create data set DATA=Challlenger (Output 10.1)                        */
/*    to obtain Outputs 10.7 and 10.8                                          */ 

proc genmod data=Challenger;
 model td/total=temp/dist=binomial link=probit 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.8  */
data prob_hat;
 set logit;
  phat=probnorm(estimate);
  pi=3.14159;
  invsqrt=1/(sqrt(2*pi));
  se_phat=invsqrt*exp(-0.5*(estimate**2))*stderr;
  prb_LcL=probnorm(LowerCL);
  prb_UcL=probnorm(UpperCL);
proc print data=prob_hat;
run;