data ethanol;
input dose y;
  x=log10(dose);
label dose='dose (mg/dL)'
      x='log10(dose)'
      y='Ethanol (mL/h)';
lines;
2.5 0.08
2.5 0.12
2.5 0.15
2.5 0.2
2.5 0.24
5.0 0.28
5.0 0.34
5.0 0.36
5.0 0.38
5.0 0.46
10.0 0.52
10.0 0.56
10.0 0.62
10.0 0.65
10.0 0.72
;       


proc glm;
 model y=x/ p clm clparm;

 proc glm;
 model y=x/ p clm clparm alpha=.1;

proc gplot data=ethanol;
 plot y*x;
 symbol i=r1clm v=star h=4 w=2;
 title 'Ethanol Data with Regression Line';
run;

data newethanol;
  set ethanol;
  lackfit=x;
proc glm data=newethanol;
 class lackfit;
 model y=x lackfit / ss1;
title1 'Data of Table 9.1 (p. 308) in Woolson and Clarke';
title2 'Lack-of-Fit Test for the Linear Model';

 
run;