PROC RANK DATA=diag NORMAL=BLOM OUT=normset;
     VAR resid;
     RANKS rankres;
PROC PLOT DATA=normset VPERCENT=50;
     PLOT rankres*resid / VREF=0;
TITLE 'NORMAL PROBABILITY PLOT OF RESIDUALS';