options ls=80;

data manpower;
input y x1 x2;
cards;
1015.09 1160.33 9390.83
1105.18 1047.67 14942.33
1598.09 4435.67 14189.08
1204.65 5797.72 4998.58
2037.35 15409.60 12134.67
2060.42 7419.25 20267.75
2400.30 38561.67 11715.00
3183.13 36047.17 18358.29
3217.26 40000.00 20000.00
2776.20 35000.00 15000.00
;

proc reg data=manpower;
 model y=x1 x2/partial;
output out=a r=resid1;


proc capability data=a noprint;
 probplot resid1;

data a;set a;
x1a=sqrt(x1);
x1b=x1**(1/3);
x1c=log(x1);
x1d=x1**(1/4);
proc reg data=a;
 model y=x1a x2/partial; 
output out=b r=resid2;

proc capability data=b noprint;
 probplot resid2; 


proc reg data=a; 
 model y=x1b x2/partial;
output out=c r=resid3;

proc capability data=c noprint;
 probplot resid3;
                     
proc reg data=a;
 model y=x1c x2/partial;
output out=d r=resid4;

proc capability data=d noprint;
 probplot resid4;  


proc reg data=a;
 model y=x1d x2/partial;
output out=e r=resid5;

proc capability data=e noprint;
 probplot resid5; 



proc reg data=a;
 model y=x1b x2/influance r vif;


run;