data tomato;
 input group $ yield;
datalines;
a 29.9
a 11.4
b 26.6
b 23.7
a 25.3
b 28.5
b 14.2
b 17.9
a 16.5
a 21.1
b 24.3
;

run;

proc boxplot;
  plot yield*group;
  title 'High resolution boxplots for each group';
run;
       
proc ttest;
  class group;
  var yield;
  title 'T-test of tomato yield';
run;
     
proc glm;
  class group;
  model yield =  group;
  title 'One-way ANOVA of tomato yield';
run;