options linesize = 80;

ods html file='o-ring.html';

data oring;
  input temp td notd total;
  cards;
53 1 0 1
57 1 0 1
58 1 0 1
63 1 0 1
66 0 1 1
67 0 3 3
68 0 1 1
69 0 1 1
70 2 2 4
72 0 1 1
73 0 1 1
75 1 1 2
76 0 2 2
78 0 1 1
79 0 1 1
81 0 1 1
;

data dummy;
  do temp = 30 to 85;
    output;
  end;

data oring;
  set oring dummy;

proc logistic data = oring;
  model td / total = temp;
  output out = oringout predicted = p l = l u = u;

proc plot data = oringout;
  plot (l p u) * temp / overlay;

proc gplot data = oringout;
  where total = .;
  symbol1 v = none i = join;
  symbol2 v = none i = join c = 'red';
  plot p * temp = 1 (l u) * temp = 2 / overlay;
run;