*NKNW964.sas, one-way random effects anova using data in Table 24.1;

options nocenter;
data a1; infile 'u:\.www\datasets512\CH24TA01.DAT';
   input rating officer;
proc print data=a1; run;

title1 'Plot of the data';
symbol1 v=circle i=none c=black;
proc gplot data=a1; 
   plot rating*officer/frame;
run;

proc means data=a1; 
   output out=a2 mean=avrate;
   var rating;
   by officer;

title1 'Plot of the means';
symbol1 v=circle i=join c=black;
proc gplot data=a2; 
   plot avrate*officer/frame;
run;

proc glm data=a1; 
   class officer;
   model rating=officer;
   random officer;
run;

proc varcomp data=a1; 
   class officer;
   model rating=officer;
run;

proc mixed data=a1 cl; 
   class officer;
   model rating=;
   random officer/vcorr;
run;