options ls=72 pagesize=50;
data a;
infile 'a:\peru.dat';
input age years weight height chin forearm calf systolic;
proc reg outest=b outvif outstb ridge=0 to 4.0 by .1;
model systolic=age years weight height chin forearm calf;
plot / ridgeplot;
data c; set b;
if _type_='RIDGESTB' or _type_='RIDGEVIF';
proc sort data=c;
by _type_;
proc gplot; by _type_;
plot (age years weight height chin forearm calf)*_RIDGE_/overlay;
data d; set b;
if _type_='RIDGE';
proc print data=d;
run;