data a;
set fly;
proc glm;
class chemical day group;
model response=chemical day group;
output out=outa r=resid p=yhat;
proc plot data=outa;
plot resid*(yhat chemical day group);
run;