options nocenter ls = 65 ps = 50;

data a;
   infile 'totca.dat' missover firstobs=2;
   input trt rep volrep ppm mlAAS mgAAS ppmEFF mgmlEFF
      mMEFF volEFF totmgEFF pcttot pctadded;
   lmMEFF = -log10(mMEFF);

proc glm; class trt rep volrep;
   model mMEFF lmMEFF = trt volrep rep(trt) volrep(trt);
   test h = trt e = rep(trt);
   test h = volrep e = volrep(trt);
   means trt/ e = rep(trt) lsd lines;
   means volrep/ e = volrep(trt) lsd lines;

proc sort; by volrep;
proc glm; by volrep;
   class trt;
   model mMEFF lmMEFF = trt;
   means trt/ lsd lines;
   lsmeans trt / out=lsm;

proc sort; by _NAME_;
proc plot data=lsm; by _NAME_;
   plot lsmean*volrep=trt;