options ls=75 nodate;

data one;
    input diet location rep y;
    cards;
1  1  1  2.10
1  2  1  2.02
1  3  1  2.16
1  4  1  1.98
1  1  2  2.32
1  2  2  2.04
1  3  2  2.18
1  4  2  1.86
2  1  1  2.24
2  2  1  2.30
2  3  1  2.22
2  4  1  1.64
2  1  2  2.22
2  2  2  2.12
2  3  2  2.18
2  4  2  1.73
3  1  1  2.28
3  2  1  2.14
3  3  1  2.26
3  4  1  1.83
3  1  2  2.24
3  2  2  2.17
3  3  2  2.21
3  4  2  1.89
;
run;

symbol1 i=join value=dot;
*symbol2 i=join value=circle;
*symbol3 i=join value=plus;
*symbol4 i=join value=diamond;
                                                                                
proc means nway noprint;
   class diet location;
   var y;
   output out=two mean=ymean;
run;

proc gplot data=two; 
   plot ymean*location=diet;
   plot ymean*diet=location;
run;

proc mixed data=one method=type3;
   class location diet;
   model y=diet/s alpha=0.1;
   random location location*diet;
   lsmeans diet/diff adj=bon;
run;