/* read data */
data dat;
   infile 'feed.dat' firstobs=2;
   input trt bwt bwg fi fe cla la;

/* adjusted treatment */
proc sort data=dat; by trt;
proc means; by trt;
   var bwg fi;
proc glm;
   class trt;
   model bwg = trt fi;
   lsmeans trt / stderr pdiff;
   output out=chinfit p=pbwg r=rbwg;
proc plot data=chinfit;
   plot bwg*fi=trt pbwg*fi='*' / overlay;