options ls=79 ps=50 nocenter;

/* table of means */
data dat;
   infile 'growth.dat' firstobs=2;
   input trt BA TDZ code $ advplt;
   ladvplt = log10(1+advplt);
   /*
   if advplt=0 then delete;
   */
   if trt=20 then delete;
proc glm;
   class BA TDZ;
   model ladvplt = BA TDZ BA*TDZ;
   lsmeans BA*TDZ / stderr pdiff;
proc anova;
   class trt;
   model ladvplt = trt;
   means trt / lsd lines;

/* one-factor analysis of variance*/
proc glm;
   class trt;
   model ladvplt = trt / ss1;
   /* trt means and standard deviations */
   means trt;
   /* multiple comparisons based on F statistics */
   means trt / lsd bon scheffe regwf waller lines;
   /* multiple comparisons based on range statistics */
   means trt / tukey snk regwq duncan lines;
   /* Fisher's LSD adjusting for sample size */
   lsmeans trt / stderr pdiff;
run;