/* Import the Excel file treeave.xls as Soil */ proc import out=soil datafile="W:\courses\stat706\treeave.xlsx" dbms=xlsx replace; run; data a; set soil; if soil in ('grass', 'grnhous', 'vermic', 'promix'); run; proc sgplot; vbox score/group=soil; run; proc glm; class soil; model score=soil/e; run;