* battery2.sas, battery experiment, Section 4.6.2 (page 96);
* includes Tukey's method for unbalanced one-way anova; 
OPTIONS LS=72;
;
DATA BATTERY;
  INPUT TYPEBAT LIFEUC ORDER;
  LINES;
   1 611  1
   2 923  2
   1 537  3
   4 476  4
   1 542  5
   1 593  6
   2 794  7
   3 445  8
   4 569  9 
   2 827 10
   2 898 11
   3 490 12
   4 480 13
   3 384 14
   4 460 15
   3 413 16
;
PROC PRINT;
;
PROC PLOT;
  PLOT LIFEUC*TYPEBAT
    / VPOS=19 HPOS=50;
;
PROC GLM;
  CLASSES TYPEBAT;
  MODEL LIFEUC = TYPEBAT;
  ESTIMATE 'DUTY'    
          TYPEBAT  1  1 -1 -1 / divisor = 2;
  ESTIMATE 'BRAND'   
          TYPEBAT  1 -1  1 -1 / divisor = 2;
  ESTIMATE 'INTERACTION' 
          TYPEBAT  1 -1 -1  1 / divisor = 2;
  MEANS TYPEBAT / TUKEY CLDIFF ALPHA=0.01;
  CONTRAST 'BRAND'
          TYPEBAT  1 -1  1 -1;
;
DATA B14; SET BATTERY;
  IF ORDER > 2;
PROC GLM;
  CLASSES TYPEBAT;
  MODEL LIFEUC = TYPEBAT;
  MEANS TYPEBAT / TUKEY CLDIFF ALPHA=0.01;