options ls=75;

data one;
   do drug=1 to 5;
      do rep=1 to 4;
         input y @;
         output;
      end;
   end;
datalines;
29.6   24.3   28.5   32  
27.3   32.6   30.8   34.8 
5.8   6.2   11   8.3 
21.6   17.4   18.3   19 
29.2   32.8   25   24.2 
;
run;

proc print;run;

proc means mean var;
   class drug;
   var y;
run;
proc glmmodel data=one;
   class drug;
   model  y=drug;
run;

proc glm data=one;
   title "Bonferroni correction for 4 contrasts";
   class drug;
   model  y=drug/clparm alpha=.0125;
   estimate "theta1" drug -1 1;
   estimate "theta2" drug -1 0 1;
   estimate "theta3" drug -1 0 0 1;
   estimate "theta4" drug -1 0 0 0 1;
run;