DATA MORLOK;                                                                    
INPUT V1 V2 V3 V4 V5 V6 V7 V8 V9 SAVPERC;                                       
v6=-1/sqrt(v6);v22=v2**2;                                                        
CARDS;                                                                          
  55.0  0.10 17.95  379.45  15.0 1327.02  4.16  0.00  10.68 19.8                
  67.0  1.00 11.94  767.52  20.0 16120.02 0.12  1.00   7.21 29.0                
  20.0  1.00 17.77  92.40   15.0 3060.24  15.58 1.00  11.56 50.0                
  66.0  0.89 20.18  810.98   0.0 1861.85   3.33 1.00   8.02 48.0                
  35.0  1.00  7.57  220.00   0.0  806.69   4.38 1.00   6.11 30.0                
  69.0  0.45 19.57  468.22  16.0 1014.00   3.52 0.00   7.28 30.0                
  35.0  0.43 11.00  184.30   0.0   35.78   0.76 1.00  6.14   0.1                
  48.0  1.00 13.94  206.44   0.0   24.15   1.20 1.00   7.37 27.0                
  45.0  0.83 14.21  255.00   0.0  38.99    0.47 0.88   7.34  3.1                
  45.0  0.80 13.32  145.00   0.0  27.38    0.47 0.67   7.19  0.1                
  26.0  1.00 12.63  173.00   0.0  25.08    0.47 1.00   6.09 30.2                
  20.0  1.00 19.03  405.10   0.0 1606.77   9.57 0.00   9.28 45.9                
  45.0  0.00 29.93   74.41  15.0 1327.02   4.16 0.00   8.34 49.0                
  45.0  0.00 33.29   38.80  25.0 1327.02   4.16 1.00   6.38 50.0                
  60.0  0.75 10.13  297.20  12.0 7869.54  16.68 1.00   7.80 43.0                
  45.0  0.71 16.83  360.27  25.0 5179.78  12.36 1.00   8.84 22.9                
   9.0  1.00 13.90   27.93   0.0  514.62   2.00 0.00   6.94 10.0                
;                                                                               
OPTIONS LS=70 NODATE NONUMBER;                                              
PROC REG OUTEST=B;                                                              
MODEL SAVPERC=V1 V2 V3 V4 V5 V6 V7 V8 V9;                                   
OUTPUT OUT=C R=R P=P;                                                           
DATA D; SET B;                                                                  
DO I=1 TO 17;                                                                   
     B1=V1; B2=V2; B3=V3;                                                       
     B4=V4; B5=V5; B6=V6;                                                       
     B7=V7; B8=V8; B9=V9;                                                       
KEEP B1 B2 B3 B4 B5 B6 B7 B8 B9;                                                
OUTPUT; END;                                                                    
DATA E; MERGE C D;                                                              
C1=B1*V1+R; C2=B2*V2+R; C3=B3*V3+R;                                             
C4=B4*V4+R; C5=B5*V5+R; C6=B6*V6+R;                                             
C7=B7*V7+R; C8=B8*V8+R; C9=B9*V9+R;                                                                                                        
PROC PLOT;*PLOT C1*V1 C2*V2 C3*V3 C4*V4 C5*V5 C6*V6 C7*V7 C8*V8 C9*V9 R*P;      
plot c2*v2/vpos=17;  title 'Component Plus Rsidual Plot Before Transformation';
run;

PROC REG OUTEST=E;                                                          
MODEL SAVPERC=V1 V2 V3 V4 V5 V6 V7 V8 V9 v22/ tol influence;                    
OUTPUT OUT=F R=R P=P;                                                         
DATA G; SET E;                                                                  
DO I=1 TO 17;                                                                   
     B1=V1; B2=V2; B3=V3;                                                       
     B4=V4; B5=V5; B6=V6;                                                       
     B7=V7; B8=V8; B9=V9;                                                       
KEEP B1 B2 B3 B4 B5 B6 B7 B8 B9;                                                
OUTPUT; END;                                                                    
DATA H; MERGE F G;                                                              
C1=B1*V1+R; C2=B2*V2+R; C3=B3*V3+R;                                             
C4=B4*V4+R; C5=B5*V5+R; C6=B6*V6+R;                                             
C7=B7*V7+R; C8=B8*V8+R; C9=B9*V9+R;                                                                                                                  
PROC PLOT;*PLOT C1*V1 C2*V2 C3*V3 C4*V4 C5*V5 C6*V6 C7*V7 C8*V8 C9*V9 R*P;      
plot c2*v2;  title 'Component Plus Residual Plot after including V2 square';
run;