data shut1;
 input temp td @@;
 datalines;
 66 0 70 1 69 0 68 0 67 0 72 0 73 0 70 0 57 1 63 1
 70 1 78 0 67 0 53 1 67 0 75 0 70 0 81 0 76 0 79 0
 75 1 76 0 58 1
 ;
data shut2;
 do i=1 to 50;
  temp=i+29; td=.;
 output;
 end;
data shut3;
 set shut1 shut2;
proc logistic descending data=shut3;
 model td = temp;
 output out=shut4 p=p;
proc sort data=shut4;
 by temp;
goptions; 
 symbol1 color=black value=dot interpol=none;
 symbol2 color=black value=none interpol=join;
proc gplot data=shut4;
 plot td*temp p*temp / overlay;
run;