proc sgplot data=iris;
vbar species / response=s_length stat=mean;
run;