/* For the purpose of this exercise, we will assume that the
first medicine is a placebo and we want only to make all pairwise
comparisons against the mean for this treatment, using Dunnett's
procedure. 

How many such comparisons are there when there are 8 treatments
plus the one control?

Note that Dunnett's procedure is based upon a distribution which 
is not covered in the text or in the notes. I only mention it
because questions about it were raised several times in lecture.*/

options ls=75 nodate;

data one; 
   array medlabel{9} $7. ("placebo","b","c","d","e","f","g","h","i");
      do i=1 to 9;
         input temp @@;
         medication=medlabel{i};
         output;
      end;
    cards;
       97.6 97.8 97.9 97.5 97.8   98 97.6 97.9   98
       97.4 97.7 98.1 97.4 97.8 97.7 97.5   98 97.6
       97.3 97.6 97.8 97.5 97.7 97.8 97.6 97.9   98
       97.5 97.7 97.8 97.6 97.7 97.9 97.5 97.9 97.9
       97.5 97.7 97.6 97.7 97.8 97.8 97.3 97.8 97.9
run;

*insert PROC GLM code here;

proc glm;
   title "Dunnett's procedure: multiple comparisons with a control";
   class medication;
   model temp=medication;
   means medication/dunnett ("placebo");
run;