options ls=80;
   
data one;
   *infile "root.dat" firstobs=2 dlm='09'x; 
   input pretrt cutting variety y;
   cards;
1  1  1  11
1  1  2  26
1  1  3  17
1  1  4  08
1  2  1  05
1  2  2  13
1  2  3  30
1  2  4  05
1  3  1  07
1  3  2  15
1  3  3  21
1  3  4  05
2  1  1  15
2  1  2  20
2  1  3  15
2  1  4  15
2  2  1  17
2  2  2  21
2  2  3  29
2  2  4  12
2  3  1  04
2  3  2  20
2  3  3  28
2  3  4  10
3  1  1  03
3  1  2  05
3  1  3  06
3  1  4  10
3  2  1  01
3  2  2  04
3  2  3  03
3  2  4  10
3  3  1  06
3  3  2  04
3  3  3  04
3  3  4  05
run;

proc glm;
   *class pretrt variety;
   class variety pretrt; 
   *model y=pretrt*variety;  /* one-way */
   model y=pretrt|variety;   /* two-factor analysis */
   *lsmeans pretrt*variety/slice=pretrt;
   lsmeans pretrt*variety/slice=variety;
   lsmeans pretrt*variety/pdiff adj=scheffe; 
   output out=two p=ymean;
   /*estimate "pt3 - pt1 for v=2" pretrt -1 0 1 
      pretrt*variety 0 -1 0 0  0 0 0 0  0 1 0 0; */
   estimate "pt3 - pt1 for v=2" pretrt -1 0 1 
      pretrt*variety 0 0 0  -1 0 1  0 0 0  0 0 0 ;
   estimate "contrast h)" pretrt*variety 1 1 -2 -1 -1 2 -1 -1 2 1 1 -2/divisor=2; 
   contrast "contrast h)" pretrt*variety 1 1 -2 -1 -1 2 -1 -1 2 1 1 -2;
run;
endsas;

/*goptions device=pslepsf colors=black;*/

symbol1 i=join l=1 value=dot;
symbol2 i=join l=1 value=square;
symbol3 i=join l=1 value=diamond;
symbol4 i=join l=1 value=plus;

options ps=20 ls=70;

*proc plot;
proc gplot;
   title "Interaction plot for 13.23";
   plot ymean*pretrt=variety;
run;