options nocenter ps=1000 ls=75 formdlim="-";

data one;
   input trt bulkdensity;
   cards;
1  1.83
1  2.01
1  1.94
1  1.79
2  1.74
2  1.68
2  1.85
2  1.72
3  1.53
3  1.60
3  1.56
3  1.62
;

proc means data=one;
   class trt;
   var bulkdensity;
   output out=rao8_4 N=n sum=bdsum mean = bdmean var=bdvar;
run;
proc print data=rao8_4;
   title "summary statistics for problem 8.4";
run;