options ls=80;

data dwaine;
input x1 x2 y;
datalines;
  68.5  16.7  174.4
  45.2  16.8  164.4
  91.3  18.2  244.2
  47.8  16.3  154.6
  46.9  17.3  181.6
  66.1  18.2  207.5
  49.5  15.9  152.8
  52.0  17.2  163.2
  48.9  16.6  145.4
  38.4  16.0  137.2
  87.9  18.3  241.9
  72.8  17.1  191.1
  88.4  17.4  232.0
  42.9  15.8  145.3
  52.5  17.8  161.1
  85.7  18.4  209.7
  41.3  16.5  146.4
  51.7  16.3  144.0
  89.6  18.1  232.6
  82.7  19.1  224.1
  52.3  16.0  166.5
;

proc means data=dwaine;
/*
     variable   n      mean           s.d.
      X1	21    62.0190476    18.6203281    
      X2        21    17.1428571     0.9703460 
      Y         21   181.9047619    36.1913039 

*/

data dwaine; 
set dwaine;
X1P=(X1-62.0190476)/(SQRT(20)*18.6203281);
X2P=(X2-17.1428571)/(SQRT(20)*0.9703460);
YP=(Y-181.9047619)/(SQRT(20)*36.1913039); 

PROC PRINT DATA=DWAINE;

PROC REG DATA=DWAINE;
 MODEL YP=X1P X2P/covb xpx i corrb;

proc reg data=dwaine;
 model y=x1 x2/ covb xpx i corrb;

proc standard data=dwaine mean=0 std=1 out=a; var y x1 x2;
data a;set a;y=y/sqrt(20);x1=x1/sqrt(20);x2=x2/sqrt(20);
proc print data=a;



run;