****************************************************;
data sf36b; set sf36a;

IF 18<=AGE AND MALE>=0 THEN DO;
SFAGE1=0;SFAGE2=0;SFAGE3=0;SFAGE4=0;SFAGE5=0;
SFAGE6=0;SFAGE7=0;SFAGE8=0;SFAGE9=0;
SFAGE10=0;SFAGE11=0;SFAGE12=0;
END;
IF 18<=AGE<=24 AND MALE=1 THEN SFAGE1=1;
IF 24<AGE<=34 AND MALE=1 THEN SFAGE2=1;
IF 34<AGE<=44 AND MALE=1 THEN SFAGE3=1;
IF 44<AGE<=54 AND MALE=1 THEN SFAGE4=1;
IF 54<AGE<=64 AND MALE=1 THEN SFAGE5=1;
IF 64<AGE AND MALE=1 THEN SFAGE6=1;
IF 18<=AGE<=24 AND MALE=0 THEN SFAGE7=1;
IF 24<AGE<=34 AND MALE=0 THEN SFAGE8=1;
IF 34<AGE<=44 AND MALE=0 THEN SFAGE9=1;
IF 44<AGE<=54 AND MALE=0 THEN SFAGE10=1;
IF 54<AGE<=64 AND MALE=0 THEN SFAGE11=1;
IF 64<AGE AND MALE=0 THEN SFAGE12=1;
*********************************************************;
RUN;

***********************************************************;
**Read in dataset containing US general population means***;
**and standard deviations for each SF-36 scale, by each****;
**age/gender group;
***********************************************************;
DATA GENPOP;
input
#1 PF1  1-5 PF2    6-10 PF3   11-15 PF4   16-20 PF5   21-25 PF6   26-30 
PF7   31-35 PF8   36-40 PF9   41-45 PF10  46-50 PF11  51-55 PF12  56-60 
#2 PF1S 1-5 PF2S   6-10 PF3S  11-15 PF4S  16-20 PF5S  21-25 PF6S  26-30 
PF7S  31-35 PF8S  36-40 PF9S  41-45 PF10S 46-50 PF11S 51-55 PF12S 56-60 
#3 RP1  1-5 RP2    6-10 RP3   11-15 RP4   16-20 RP5   21-25 RP6   26-30 
RP7   31-35 RP8   36-40 RP9   41-45 RP10  46-50 RP11  51-55 RP12  56-60 
#4 RP1S 1-5 RP2S   6-10 RP3S  11-15 RP4S  16-20 RP5S  21-25 RP6S  26-30 
RP7S  31-35 RP8S  36-40 RP9S  41-45 RP10S 46-50 RP11S 51-55 RP12S 56-60 
#5 BP1  1-5 BP2    6-10 BP3   11-15 BP4   16-20 BP5   21-25 BP6   26-30 
BP7   31-35 BP8   36-40 BP9   41-45 BP10  46-50 BP11  51-55 BP12  56-60 
#6 BP1S 1-5 BP2S   6-10 BP3S  11-15 BP4S  16-20 BP5S  21-25 BP6S  26-30 
BP7S  31-35 BP8S  36-40 BP9S  41-45 BP10S 46-50 BP11S 51-55 BP12S 56-60 
#7 GEN1 1-5  GEN2   6-10  GEN3   11-15 GEN4   16-20 GEN5   21-25 GEN6   26-30
GEN7   31-35 GEN8   36-40 GEN9   41-45 GEN10  46-50 GEN11  51-55 GEN12  56-60
#8 GEN1S 1-5 GEN2S   6-10 GEN3S  11-15 GEN4S  16-20 GEN5S  21-25 GEN6S  26-30
GEN7S  31-35 GEN8S  36-40 GEN9S  41-45 GEN10S 46-50 GEN11S 51-55 GEN12S 56-60
#9 ENFT1 1-5  ENFT2    6-10  ENFT3   11-15  ENFT4   16-20 ENFT5   21-25 ENFT6   26-30 
ENFT7   31-35  ENFT8  36-40  ENFT9   41-45  ENFT10  46-50 ENFT11  51-55 ENFT12  56-60 
#10 ENFT1S 1-5 ENFT2S   6-10 ENFT3S  11-15 ENFT4S  16-20 ENFT5S  21-25 ENFT6S  26-30 
ENFT7S  31-35 ENFT8S  36-40 ENFT9S  41-45 ENFT10S 46-50 ENFT11S 51-55 ENFT12S 56-60 
#11SF1    1-5 SF2    6-10   SF3   11-15   SF4   16-20   SF5   21-25   SF6   26-30   
SF7   31-35   SF8   36-40   SF9   41-45   SF10  46-50   SF11  51-55   SF12  56-60   
#12 SF1S 1-5  SF2S   6-10   SF3S  11-15   SF4S  16-20   SF5S  21-25   SF6S  26-30   
SF7S  31-35   SF8S  36-40   SF9S  41-45   SF10S 46-50   SF11S 51-55   SF12S 56-60   
#13           RE1     1-5   RE2    6-10   RE3   11-15   RE4   16-20   RE5   21-25   RE6   26-30   
RE7   31-35   RE8   36-40   RE9   41-45   RE10  46-50   RE11  51-55   RE12  56-60   
#14           RE1S    1-5   RE2S   6-10   RE3S  11-15   RE4S  16-20   RE5S  21-25   RE6S  26-30   
RE7S  31-35   RE8S  36-40   RE9S  41-45   RE10S 46-50   RE11S 51-55   RE12S 56-60   
#15 MH1     1-5   MH2    6-10   MH3   11-15   MH4   16-20   MH5   21-25   MH6   26-30   
MH7   31-35   MH8   36-40   MH9   41-45   MH10  46-50   MH11  51-55   MH12  56-60   
#16           MH1S    1-5   MH2S   6-10   MH3S  11-15   MH4S  16-20   MH5S  21-25   MH6S  26-30   
MH7S  31-35   MH8S  36-40   MH9S  41-45   MH10S 46-50   MH11S 51-55   MH12S 56-60   ;
datalines;
94.1494.9091.3986.5079.8865.7990.1889.1288.0682.8673.0961.86
16.3011.4414.6820.3725.4728.3120.0418.7217.7021.7226.7328.95
93.5091.8689.7685.5876.0359.7284.9186.7383.6579.9371.6156.11
21.3921.0424.7530.2436.6642.5130.7327.9932.2135.3839.8442.53
79.6283.1079.3674.1868.5068.7682.0079.6174.8572.1466.6463.44
21.4718.2221.2124.8126.1025.3721.3120.9422.7423.3425.2627.12
76.9579.4077.5573.1666.5858.6276.4874.8074.2570.4862.8761.64
17.8717.1715.9217.8923.2722.0518.6617.2419.4420.5823.3722.08
65.4164.6665.5063.0563.0057.8059.7158.0459.4360.6258.0855.46
19.1118.9918.6720.4421.3822.5520.0920.8519.7221.3223.4223.51
86.0985.6688.5485.5483.5679.6681.7384.0683.0782.7179.4377.00
20.1819.8318.0922.8521.9626.0020.9821.6623.2720.8427.0227.69
87.4982.1585.5285.4281.1476.9478.6382.3280.0881.9279.5173.38
27.5031.7828.1229.2534.0237.4833.8631.3033.8833.3434.6439.66
78.0274.1376.9876.3876.8777.3771.5372.4573.3274.3673.4074.71
16.0517.2416.4217.6218.7117.4219.4418.6216.7918.0819.7419.88
;
RUN;


****************************************************************;
**T-statistics for the difference between SF-36 sample scores***;
**compared to the US general population, for each US age/gender*;
**group***;
****************************************************************;


DATA TEST;
 IF _N_=1 THEN SET GENPOP;
 SET sf36b;
PFISFM=(SFAGE1*PF1)+ (SFAGE2*PF2)+ (SFAGE3*PF3)+ 
       (SFAGE4*PF4)+ (SFAGE5*PF5)+ (SFAGE6*PF6)+ 
       (SFAGE7*PF7)+ (SFAGE8*PF8)+ (SFAGE9*PF9)+ 
       (SFAGE10*PF10)+ (SFAGE11*PF11)+ (SFAGE12*PF12); 
PFIS=(SFAGE1*PF1S*PF1S)+ (SFAGE2*PF2S*PF2S)+
      (SFAGE3*PF3S*PF3S)+ (SFAGE4*PF4S*PF4S)+
      (SFAGE5*PF5S*PF5S)+ (SFAGE6*PF6S*PF6S)+
      (SFAGE7*PF7S*PF7S)+ (SFAGE8*PF8S*PF8S)+
      (SFAGE9*PF9S*PF9S)+ (SFAGE10*PF10S*PF10S)+
      (SFAGE11*PF11S*PF11S)+ (SFAGE12*PF12S*PF12S);
IF SFAGE1=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF1S;
  ELSE IF SFAGE2=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF2S;
  ELSE IF SFAGE3=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF3S;
  ELSE IF SFAGE4=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF4S;
  ELSE IF SFAGE5=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF5S;
  ELSE IF SFAGE6=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF6S;
  ELSE IF SFAGE7=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF7S;
  ELSE IF SFAGE8=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF8S;
  ELSE IF SFAGE9=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF9S;
  ELSE IF SFAGE10=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF10S;
  ELSE IF SFAGE11=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF11S;
  ELSE IF SFAGE12=1 THEN ZPHY10=(PHYFUN10-PFISFM)/PF12S;
RPSFM=(SFAGE1*RP1)+(SFAGE2*RP2)+(SFAGE3*RP3)+          
       (SFAGE4*RP4)+ (SFAGE5*RP5)+ (SFAGE6*RP6)+       
       (SFAGE7*RP7)+ (SFAGE8*RP8)+ (SFAGE9*RP9)+       
       (SFAGE10*RP10)+ (SFAGE11*RP11)+ (SFAGE12*RP12); 
RPS= (SFAGE1*RP1S*RP1S)+ (SFAGE2*RP2S*RP2S)+
      (SFAGE3*RP3S*RP3S)+ (SFAGE4*RP4S*RP4S)+
      (SFAGE5*RP5S*RP5S)+ (SFAGE6*RP6S*RP6S)+
      (SFAGE7*RP7S*RP7S)+ (SFAGE8*RP8S*RP8S)+
      (SFAGE9*RP9S*RP9S)+ (SFAGE10*RP10S*RP10S)+
      (SFAGE11*RP11S*RP11S)+ (SFAGE12*RP12S*RP12S);
IF SFAGE1=1 THEN ZRP=(ROLEP4-RPSFM)/RP1S;
  ELSE IF SFAGE2=1 THEN ZRP=(ROLEP4-RPSFM)/RP2S;
  ELSE IF SFAGE3=1 THEN ZRP=(ROLEP4-RPSFM)/RP3S;
  ELSE IF SFAGE4=1 THEN ZRP=(ROLEP4-RPSFM)/RP4S;
  ELSE IF SFAGE5=1 THEN ZRP=(ROLEP4-RPSFM)/RP5S;
  ELSE IF SFAGE6=1 THEN ZRP=(ROLEP4-RPSFM)/RP6S;
  ELSE IF SFAGE7=1 THEN ZRP=(ROLEP4-RPSFM)/RP7S;
  ELSE IF SFAGE8=1 THEN ZRP=(ROLEP4-RPSFM)/RP8S;
  ELSE IF SFAGE9=1 THEN ZRP=(ROLEP4-RPSFM)/RP9S;
  ELSE IF SFAGE10=1 THEN ZRP=(ROLEP4-RPSFM)/RP10S;
  ELSE IF SFAGE11=1 THEN ZRP=(ROLEP4-RPSFM)/RP11S;
  ELSE IF SFAGE12=1 THEN ZRP=(ROLEP4-RPSFM)/RP12S;
BPSFM=(SFAGE1*BP1)+(SFAGE2*BP2)+(SFAGE3*BP3)+          
       (SFAGE4*BP4)+ (SFAGE5*BP5)+ (SFAGE6*BP6)+       
       (SFAGE7*BP7)+ (SFAGE8*BP8)+ (SFAGE9*BP9)+       
       (SFAGE10*BP10)+ (SFAGE11*BP11)+ (SFAGE12*BP12); 
BPS=(SFAGE1*BP1S*BP1S)+ (SFAGE2*BP2S*BP2S)+
      (SFAGE3*BP3S*BP3S)+ (SFAGE4*BP4S*BP4S)+
      (SFAGE5*BP5S*BP5S)+ (SFAGE6*BP6S*BP6S)+
      (SFAGE7*BP7S*BP7S)+ (SFAGE8*BP8S*BP8S)+
      (SFAGE9*BP9S*BP9S)+ (SFAGE10*BP10S*BP10S)+
      (SFAGE11*BP11S*BP11S)+ (SFAGE12*BP12S*BP12S);

IF SFAGE1=1 THEN ZBP=(SFPAIN-BPSFM)/BP1S;
  ELSE IF SFAGE2=1 THEN ZBP=(SFPAIN-BPSFM)/BP2S;
  ELSE IF SFAGE3=1 THEN ZBP=(SFPAIN-BPSFM)/BP3S;
  ELSE IF SFAGE4=1 THEN ZBP=(SFPAIN-BPSFM)/BP4S;
  ELSE IF SFAGE5=1 THEN ZBP=(SFPAIN-BPSFM)/BP5S;
  ELSE IF SFAGE6=1 THEN ZBP=(SFPAIN-BPSFM)/BP6S;
  ELSE IF SFAGE7=1 THEN ZBP=(SFPAIN-BPSFM)/BP7S;
  ELSE IF SFAGE8=1 THEN ZBP=(SFPAIN-BPSFM)/BP8S;
  ELSE IF SFAGE9=1 THEN ZBP=(SFPAIN-BPSFM)/BP9S;
  ELSE IF SFAGE10=1 THEN ZBP=(SFPAIN-BPSFM)/BP10S;
  ELSE IF SFAGE11=1 THEN ZBP=(SFPAIN-BPSFM)/BP11S;
  ELSE IF SFAGE12=1 THEN ZBP=(SFPAIN-BPSFM)/BP12S;
GENSFM=(SFAGE1*GEN1)+(SFAGE2*GEN2)+(SFAGE3*GEN3)+      
       (SFAGE4*GEN4)+ (SFAGE5*GEN5)+ (SFAGE6*GEN6)+    
       (SFAGE7*GEN7)+ (SFAGE8*GEN8)+ (SFAGE9*GEN9)+    
       (SFAGE10*GEN10)+ (SFAGE11*GEN11)+ (SFAGE12*GEN12);
GENS=(SFAGE1*GEN1S*GEN1S)+ (SFAGE2*GEN2S*GEN2S)+
      (SFAGE3*GEN3S*GEN3S)+ (SFAGE4*GEN4S*GEN4S)+
      (SFAGE5*GEN5S*GEN5S)+ (SFAGE6*GEN6S*GEN6S)+
      (SFAGE7*GEN7S*GEN7S)+ (SFAGE8*GEN8S*GEN8S)+
      (SFAGE9*GEN9S*GEN9S)+ (SFAGE10*GEN10S*GEN10S)+
      (SFAGE11*GEN11S*GEN11S)+ (SFAGE12*GEN12S*GEN12S);
IF SFAGE1=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN1S;
  ELSE IF SFAGE2=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN2S;
  ELSE IF SFAGE3=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN3S;
  ELSE IF SFAGE4=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN4S;
  ELSE IF SFAGE5=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN5S;
  ELSE IF SFAGE6=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN6S;
  ELSE IF SFAGE7=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN7S;
  ELSE IF SFAGE8=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN8S;
  ELSE IF SFAGE9=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN9S;
  ELSE IF SFAGE10=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN10S;
  ELSE IF SFAGE11=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN11S;
  ELSE IF SFAGE12=1 THEN ZGENH=(SFGENH5-GENSFM)/GEN12S;
ENFTSFM=(SFAGE1*ENFT1)+(SFAGE2*ENFT2)+(SFAGE3*ENFT3)+    
       (SFAGE4*ENFT4)+ (SFAGE5*ENFT5)+ (SFAGE6*ENFT6)+   
       (SFAGE7*ENFT7)+ (SFAGE8*ENFT8)+ (SFAGE9*ENFT9)+   
       (SFAGE10*ENFT10)+ (SFAGE11*ENFT11)+ (SFAGE12*ENFT12);
ENFTS=(SFAGE1*ENFT1S*ENFT1S)+ (SFAGE2*ENFT2S*ENFT2S)+
      (SFAGE3*ENFT3S*ENFT3S)+ (SFAGE4*ENFT4S*ENFT4S)+
      (SFAGE5*ENFT5S*ENFT5S)+ (SFAGE6*ENFT6S*ENFT6S)+
      (SFAGE7*ENFT7S*ENFT7S)+ (SFAGE8*ENFT8S*ENFT8S)+
      (SFAGE9*ENFT9S*ENFT9S)+ (SFAGE10*ENFT10S*ENFT10S)+
      (SFAGE11*ENFT11S*ENFT11S)+ (SFAGE12*ENFT12S*ENFT12S);
IF SFAGE1=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT1S;
  ELSE IF SFAGE2=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT2S;
  ELSE IF SFAGE3=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT3S;
  ELSE IF SFAGE4=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT4S;
  ELSE IF SFAGE5=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT5S;
  ELSE IF SFAGE6=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT6S;
  ELSE IF SFAGE7=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT7S;
  ELSE IF SFAGE8=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT8S;
  ELSE IF SFAGE9=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT9S;
  ELSE IF SFAGE10=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT10S;
  ELSE IF SFAGE11=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT11S;
  ELSE IF SFAGE12=1 THEN ZENFT=(ENFAT4-ENFTSFM)/ENFT12S;
SFSFM=(SFAGE1*SF1)+(SFAGE2*SF2)+(SFAGE3*SF3)+               
       (SFAGE4*SF4)+ (SFAGE5*SF5)+ (SFAGE6*SF6)+            
       (SFAGE7*SF7)+ (SFAGE8*SF8)+ (SFAGE9*SF9)+            
       (SFAGE10*SF10)+ (SFAGE11*SF11)+ (SFAGE12*SF12);      

SFS=(SFAGE1*SF1S*SF1S)+ (SFAGE2*SF2S*SF2S)+
      (SFAGE3*SF3S*SF3S)+ (SFAGE4*SF4S*SF4S)+
      (SFAGE5*SF5S*SF5S)+ (SFAGE6*SF6S*SF6S)+
      (SFAGE7*SF7S*SF7S)+ (SFAGE8*SF8S*SF8S)+
      (SFAGE9*SF9S*SF9S)+ (SFAGE10*SF10S*SF10S)+
      (SFAGE11*SF11S*SF11S)+ (SFAGE12*SF12S*SF12S);
IF SFAGE1=1 THEN ZSF=(SOCFUN2-SFSFM)/SF1S;
  ELSE IF SFAGE2=1 THEN ZSF=(SOCFUN2-SFSFM)/SF2S;
  ELSE IF SFAGE3=1 THEN ZSF=(SOCFUN2-SFSFM)/SF3S;
  ELSE IF SFAGE4=1 THEN ZSF=(SOCFUN2-SFSFM)/SF4S;
  ELSE IF SFAGE5=1 THEN ZSF=(SOCFUN2-SFSFM)/SF5S;
  ELSE IF SFAGE6=1 THEN ZSF=(SOCFUN2-SFSFM)/SF6S;
  ELSE IF SFAGE7=1 THEN ZSF=(SOCFUN2-SFSFM)/SF7S;
  ELSE IF SFAGE8=1 THEN ZSF=(SOCFUN2-SFSFM)/SF8S;
  ELSE IF SFAGE9=1 THEN ZSF=(SOCFUN2-SFSFM)/SF9S;
  ELSE IF SFAGE10=1 THEN ZSF=(SOCFUN2-SFSFM)/SF10S;
  ELSE IF SFAGE11=1 THEN ZSF=(SOCFUN2-SFSFM)/SF11S;
  ELSE IF SFAGE12=1 THEN ZSF=(SOCFUN2-SFSFM)/SF12S;
RESFM=(SFAGE1*RE1)+(SFAGE2*RE2)+(SFAGE3*RE3)+               
       (SFAGE4*RE4)+ (SFAGE5*RE5)+ (SFAGE6*RE6)+            
       (SFAGE7*RE7)+ (SFAGE8*RE8)+ (SFAGE9*RE9)+            
       (SFAGE10*RE10)+ (SFAGE11*RE11)+ (SFAGE12*RE12);      
RES=(SFAGE1*RE1S*RE1S)+ (SFAGE2*RE2S*RE2S)+
      (SFAGE3*RE3S*RE3S)+ (SFAGE4*RE4S*RE4S)+
      (SFAGE5*RE5S*RE5S)+ (SFAGE6*RE6S*RE6S)+
      (SFAGE7*RE7S*RE7S)+ (SFAGE8*RE8S*RE8S)+
      (SFAGE9*RE9S*RE9S)+ (SFAGE10*RE10S*RE10S)+
      (SFAGE11*RE11S*RE11S)+ (SFAGE12*RE12S*RE12S);
IF SFAGE1=1 THEN ZRE=(ROLEE3-RESFM)/RE1S;
  ELSE IF SFAGE2=1 THEN ZRE=(ROLEE3-RESFM)/RE2S;
  ELSE IF SFAGE3=1 THEN ZRE=(ROLEE3-RESFM)/RE3S;
  ELSE IF SFAGE4=1 THEN ZRE=(ROLEE3-RESFM)/RE4S;
  ELSE IF SFAGE5=1 THEN ZRE=(ROLEE3-RESFM)/RE5S;
  ELSE IF SFAGE6=1 THEN ZRE=(ROLEE3-RESFM)/RE6S;
  ELSE IF SFAGE7=1 THEN ZRE=(ROLEE3-RESFM)/RE7S;
  ELSE IF SFAGE8=1 THEN ZRE=(ROLEE3-RESFM)/RE8S;
  ELSE IF SFAGE9=1 THEN ZRE=(ROLEE3-RESFM)/RE9S;
  ELSE IF SFAGE10=1 THEN ZRE=(ROLEE3-RESFM)/RE10S;
  ELSE IF SFAGE11=1 THEN ZRE=(ROLEE3-RESFM)/RE11S;
  ELSE IF SFAGE12=1 THEN ZRE=(ROLEE3-RESFM)/RE12S;
MHSFM=(SFAGE1*MH1)+(SFAGE2*MH2)+(SFAGE3*MH3)+               
       (SFAGE4*MH4)+ (SFAGE5*MH5)+ (SFAGE6*MH6)+            
       (SFAGE7*MH7)+ (SFAGE8*MH8)+ (SFAGE9*MH9)+            
       (SFAGE10*MH10)+ (SFAGE11*MH11)+ (SFAGE12*MH12);      
MHS=(SFAGE1*MH1S*MH1S)+ (SFAGE2*MH2S*MH2S)+
      (SFAGE3*MH3S*MH3S)+ (SFAGE4*MH4S*MH4S)+
      (SFAGE5*MH5S*MH5S)+ (SFAGE6*MH6S*MH6S)+
      (SFAGE7*MH7S*MH7S)+ (SFAGE8*MH8S*MH8S)+
      (SFAGE9*MH9S*MH9S)+ (SFAGE10*MH10S*MH10S)+
      (SFAGE11*MH11S*MH11S)+ (SFAGE12*MH12S*MH12S);
IF SFAGE1=1 THEN ZMHI=(EMOT5-MHSFM)/MH1S;
  ELSE IF SFAGE2=1 THEN ZMHI=(EMOT5-MHSFM)/MH2S;
  ELSE IF SFAGE3=1 THEN ZMHI=(EMOT5-MHSFM)/MH3S;
  ELSE IF SFAGE4=1 THEN ZMHI=(EMOT5-MHSFM)/MH4S;
  ELSE IF SFAGE5=1 THEN ZMHI=(EMOT5-MHSFM)/MH5S;
  ELSE IF SFAGE6=1 THEN ZMHI=(EMOT5-MHSFM)/MH6S;
  ELSE IF SFAGE7=1 THEN ZMHI=(EMOT5-MHSFM)/MH7S;
  ELSE IF SFAGE8=1 THEN ZMHI=(EMOT5-MHSFM)/MH8S;
  ELSE IF SFAGE9=1 THEN ZMHI=(EMOT5-MHSFM)/MH9S;
  ELSE IF SFAGE10=1 THEN ZMHI=(EMOT5-MHSFM)/MH10S;
  ELSE IF SFAGE11=1 THEN ZMHI=(EMOT5-MHSFM)/MH11S;
  ELSE IF SFAGE12=1 THEN ZMHI=(EMOT5-MHSFM)/MH12S;

IF PHYFUN10=. THEN PFISFM=.;  IF PFISFM=.  THEN PHYFUN10=.;
IF ROLEP4=.   THEN RPSFM=.;   IF RPSFM=.   THEN ROLEP4=.;
IF PAIN2=.    THEN BPSFM=.;   IF BPSFM=.   THEN PAIN2=.;
IF BPSFM=.    THEN SFPAIN=.;                        
IF GENH5=.    THEN GENSFM=.;  IF GENSFM=.  THEN GENH5=.;
IF GENSFM=.   THEN SFGENH5=.;                      
IF ENFAT4=.   THEN ENFTSFM=.; IF ENFTSFM=. THEN ENFAT4=.;
IF SOCFUN2=.  THEN SFSFM=.;   IF SFSFM=.   THEN SOCFUN2=.;
IF ROLEE3=.   THEN RESFM=.;   IF RESFM=.   THEN ROLEE3=.;
IF EMOT5=.    THEN MHSFM=.;   IF MHSFM=.   THEN EMOT5=.;

label pfisfm="Physical functioning (general pop.)";
label rpsfm="Role limitation physical (general pop.)";
label bpsfm="Pain (general pop.)";
label gensfm="General health (general pop.)";
label mhsfm="Emotional well being (general pop.)";
label resfm="Role limitation emotional (general pop.)";
label enftsfm="Energy (general pop.)";
label sfsfm="Social functioning (general pop.)";
RUN;
*****************************************************;      

DATA TEST;SET TEST;RUN;
PROC MEANS;                                                 
VAR                                                         
PHYFUN10 ROLEP4 PAIN2 SFPAIN GENH5 SFGENH5
EMOT5 ROLEE3 ENFAT4 SOCFUN2;RUN;
PROC MEANS N MEAN;                                          
VAR                                                         
PFISFM RPSFM BPSFM GENSFM
MHSFM RESFM ENFTSFM SFSFM;RUN;
*****************************************************;      
PROC MEANS T PRT;                                           
 VAR ZPHY10 ZRP ZBP ZGENH ZENFT ZSF ZRE ZMHI;RUN;