data dirt;
infile 'S:\Courses\stat-renaes\Stat404\soil.csv' dsd firstobs=2;
input soil mix;
run;

proc print data=dirt;
run;

* PROC REG autpomatically gives the diagnostic plots for checking assumptions;
* PROC GLM needs a PLOTS= option to do the same;

proc glm data=dirt plots=(diagnostics residuals);
class mix;
model soil=mix;
means mix;
run;