options nocenter ls = 65 ps = 24;

data a;
   infile 'solca.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;