/* Sections 7.2.2.3 (p. 210) Blood pressure in dogs example */

data blood; input drug disease @;
do i = 1 to 6;
  input y @;
  output;
end;
cards;
1 1 42 44 36 13 19 22
1 2 33  . 26  . 33 21
1 3 31 -3  . 25 25 24
2 1 28  . 23 34 42 13
2 2  . 34 33 31  . 36
2 3  3 26 28 32  4 16
3 1  .  .  1 29  . 19
3 2  . 11  9  7  1 -6
3 3 21  1  .  9  3  .
4 1 24  .  9 22 -2 15
4 2 27 12 12 -5 16 15
4 3 22  7 25  5 12  .
;
proc glm;
  class drug disease;
  model y = drug disease;
  means drug/dunnett('4');
proc mixed;
  class drug disease;
  model y = drug disease;
  lsmeans drug/pdiff=control('4') adjust=dunnett cl;
  run;