options ls=72;
options pagesize=50;
data a;
infile 'a:peru.dat';
input age years weight height chin forearm calf systolic;
proc reg;
model systolic=age years weight height chin forearm calf/selection=rsquare best=2;
proc reg;
model systolic=age years weight height chin forearm calf/selection=cp;
proc reg;
model systolic=age years weight height chin forearm calf/selection=stepwise
 slentry=.15 slstay=.15;
proc reg;
model systolic=years weight;
output out=m2 press=d2i;
proc reg;model systolic= years weight chin;
output out=m3 press=d3i;
data b; merge m2 m3;
proc univariate;
var d2i d3i;
run;