options nocenter ls=75;

data potency;
 input sub y @@;
 cards;
 1 29 1 28 1 23 1 26 1 26 
 1 19 1 25 1 29 1 26 1 28
 2 17 2 25 2 24 2 19 2 28
 2 21 2 20 2 25 2 19 2 24 
 3 17 3 16 3 21 3 22 3 23
 3 18 3 20 3 17 3 25 3 21
 4 18 4 20 4 25 4 24 4 16 
 4 20 4 20 4 17 4 19 4 17
 ;

proc glm;
 class sub;
 model y = sub;
 means sub / lsd bon tukey snk scheffe;
run;
quit;

data oxygen;
 input loc y @@;
 if loc=1 then dist=0;  if loc=2 then dist=10;
 if loc=3 then dist=20; if loc=4 then dist=30;
 cards;
 1 4 1 5 1 6 
 2 6 2 6 2 6 
 3 7 3 8 3 9 
 4 8 4 9 4 10
 ;

 proc glm;
  class loc;
  model y = loc;
  contrast 'linear' loc -3 -1 1 3;
  contrast 'quadratic' loc 1 -1 -1 1;
  contrast 'cubic' loc -1 3 -3 1;

proc reg;
  model y = dist;

run;
quit;