options nocenter;
goptions colors=(none);

data table5;
  infile 'I:\.www\datasets502\example5-1.dat';
  input indiv sbp age;

symbol1 v=circle;

proc reg;
 model sbp=age / dw clb clm cli alpha=.05;
 plot sbp*age / conf pred;
 plot r.*nqq.;
 plot r.*p.;
 output out=new p=pred r=res;

 proc print;

 proc gplot;
  plot res*pred / vref=0 frame;

proc univariate normal pctldef=4;
 var res;
 histogram res / normal kernel (L=2);
 qqplot res / normal (L=1 mu=est sigma=est);
run;