/* Use proc logistic with a "freq" statement and a single dependent
 * variable. */

options linesize = 80;

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

data oring;
  input temp td notd total;
  distress = 0;
  count    = notd;
  if count > 0 then output;
  distress = 1;
  count    = td;
  if count > 0 then output;
  drop 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
;

proc print data = oring;

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

data oring;
  set oring dummy;

proc logistic data = oring;
  freq count;
  model distress = 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 distress = .;
  symbol1 v = none i = join;
  symbol2 v = none i = join c = 'red';
  plot p * temp = 1 (l u) * temp = 2 / overlay;
run;