options ls=80;

data kenton;
input y design $ store $;
cards;
  11  1  1
  17  1  2
  16  1  3
  14  1  4
  15  1  5
  12  2  1
  10  2  2
  15  2  3
  19  2  4
  11  2  5
  23  3  1
  20  3  2
  18  3  3
  17  3  4
  27  4  1
  33  4  2
  22  4  3
  26  4  4
  28  4  5
;

proc glm data=kenton;
 class design;
  model y=design;

contrast 'Design 1 vs Design 2' design 1 -1 0 0;
contrast 'Design 1 vs Design 3' design 1 0 -1 0;  
contrast 'Design 2 vs Design 4' design 0 1 0 -1; 
contrast 'Design 3 vs Design 4' design 0 0 1 -1;


estimate 'LSM overall' intercept 1;
estimate 'LSM design 1' intercept 1 design 1;
estimate 'LSM design 2' intercept 1 design 0 1;
estimate 'LSM design 3' intercept 1 design 0 0 1;
estimate 'LSM design 4' intercept 1 design 0 0 0 1;
estimate 'LSM design 1-2' design 1 -1 0 0;
estimate 'LSM design 1-3' design 1 0 -1 0;
estimate 'LSM design 2-4' design 0 1 0 -1;
estimate 'LSM design 3-4' design 0 0 1 -1; 
estimate 'LSM  (1+2)-(3+4)' design 1 1 -1 -1;
estimate 'LSM  (1+2)-(3+4)' design 0.5 0.5 -0.5 -0.5;
estimate 'LSM  (1+3)-(2+4)' design 0.5 -0.5 0.5 -0.5;

means design/tukey scheffe bon lsd duncan waller;    

means design/tukey alpha=0.1;


proc glm data=kenton;
 class design;
  model y=design;

estimate 'Test 1:' design 1 -1 0 0;
estimate 'Test 2:' design 1 0 -1 0;
estimate 'Test 3:' design 0 1 0 -1;
estimate 'Test 4:' design 0 0 1 -1;
estimate 'Test 5:' design 0.5 0.5 -0.5 -0.5;
estimate 'Test 6:' design 0.5 -0.5 0.5 -0.5;