data cars93a; set ma2612.cars93a; d2=displace**2; h2=hp**2; r2=rpm**2; run; proc reg data=cars93a; model highmpg= hp displace rpm d2 h2 r2; output out=resid_out r=resid rstudent=studres p=fit; run; * Labeled Plot Studentized residuals versus fitted values; proc sgplot data=resid_out; scatter x=fit y=studres/datalabel=model; run; * Partially Labeled Plot Studentized residuals versus fitted values; data resid_out; set resid_out; label=model; if abs(studres)<2 then label=''; run; proc sgplot data=resid_out; scatter x=fit y=studres/datalabel=label; run;