libname g350 "c:\sasdat\g350";                                                                                                          
options linesize=100;                                                                                                                   
options pagesize=200;    

data g350.van01(compress=yes);                                                                                                          
        set g350.temp01;                                                                                                                
                                                                                                                                        
        p_04=(v7+v26)/v5;                       label p_04="population under age 5";                                                    
        p_15=(v8+v9+v27+v28)/v5;                label p_15="population share 5-14";                                                     
        p_20=(v10+v29)/v5;                      label p_20="population share 15-19";                                                    
        p_75=(v20+v21+v39+v40)/v5;              label p_75="population share 65-74";                                                    
        p_85=(v22+v23+v24+v41+v42+v43)/v5;      label p_85="population share over 75";                                                  
                                                                                                                                        
        pchange=v3;             label pchange="percentage population change";                                                           
        single=v45/v44;         label single="single population share";                                                                 
        divorce=v48/v44;        label divorce="divorced population share";                                                              
        separ=v47/v44;          label separ="separated population share";                                                               
        mckids=v57/v53;         label mckids="married couples with children";                                                           
        flone=v68/v53;          label flone="female lone parent households";                                                            
                                                                                                                                        
        hown=v99/v96;           label hown="homeownership rate";                                                                        
                                                                                                                                        
        vm_ch=v704/v702;        label vm_ch="vismin Chinese";                                                                           
        vm_sa=v705/v702;        label vm_sa="vismin South Asian";                                                                       
        vm_bl=v706/v702;        label vm_bl="vismin Black";                                                                             
        vm_fl=v707/v702;        label vm_fl="vismin Filipino";                                                                          
        vm_se=v709/v702;        label vm_se="vismin Southeast Asian";                                                                   
                                                                                                                                        
        lang=v217/v213;         label lang="share without knowledge of eng or french";                                                  
        aborig=v718/v717;       label aborig="total aboriginal identity population";                                                    
                                                                                                                                        
        lico=v1620;             label lico="total incidence of low income";                                                             
        mhhinc=(v1634/49940);   label mhhinc="median household income ratio";                                                           
        elite=v1632/v1621;      label elite="share of households over 100k";                                                            
        avval=(v1674/294847);   label avval="average dwelling value ratio";                                                             
        unemp=v742;             label unemp="unemployment rate";                                                                        
                                                                                                                                        
                                                                                                                                        
        mig_p=v387/v381;        label mig_p="intraprovincial movers 1 year";                                                            
        mig_c=v388/v381;        label mig_c="interprovincial movers 1 year";                                                            
        mig_i=v389/v381;        label mig_i="international movers 1 year";                                                              
                                                                                                                                        
        immig=v507/v402;        label immig="recent immigrants, 1996-2001";                                                             
                                                                                                                                        
        drivers=(v1255+v1264)/v1253;            label drivers="share of workers driving";                                               
        transit=(v1257+v1266)/v1253;            label transit="share of workers transit";                                               
        wbike=(v1258+v1259+v1267+v1268)/v1253;  label wbike="share of workers walking/biking";                                          
                                                                                                                                        
        rcost=v1668/v1666;      label rcost="renters paying >30% for rent";                                                             
        ocost=v1672/v1670;      label ocost="owners paying >30% for rent";                                                              
        b46=v105/v96;           label b46="share units built before 1946";                                                              
        b60=v106/v96;           label b60="share units built 1946-1960";                                                                
        b70=v107/v96;           label b70="share units built 1961-1970";                                                                
        b80=v108/v96;           label b80="share units built 1971-1980";                                                                
        b95=v110/v96;           label b95="share units built 1991-1995";                                                                
        b01=v111/v96;           label b01="share units built 1996-2001";                                                                
                                                                                                                                        
        h_sf=v113/v112;                 label h_sf="single family housing";                                                             
        h_ro=(v114+v115+v116)/v112;     label h_ro="semi-detached, row, duplex housing";                                                
        h_hi=(v117)/v112;               label h_hi="buildings 5+ stories";                                                              
                                                                                                                                        
        hs_hrs=(v1298+v1299)/v1293;             label hs_hrs="all persons doing 30+hrs housework";                                      
        hf_hrs=(v1312+v1313)/(v1298+v1299);     label hf_hrs="females as share of 30+hrs housework";                                    
        cc_hrs=(v1319+v1320)/v1314;             label cc_hrs="all persons doing 30+hrs childcare";                                      
        cf_hrs=(v1333+v1334)/(v1319+v1320);     label cf_hrs="females as share of 30+hrs childcare";                                    
        el_hrs=(v1341)/v1335;                   label el_hrs="all persons doing 20+hrs elder care";                                     
                                                                                                                                        
        run;                                                                                                                            
                                                                                                                                        
                               
proc factor data=g350.van01 n=10 score rotate=varimax out=factemp;                                                                      
        where _NAME_ ne 't0' and _NAME_ ne 't251';                                                                                      
        var                                                                                                                             
        pchange p_04 p_15 p_20 p_75 p_85                                                                                                
        lico mhhinc elite unemp                                                                                                         
        avval hown h_sf h_ro h_hi rcost ocost b46 b60 b70 b80 b95 b01                                                                   
        single divorce separ mckids flone                                                                                               
        hs_hrs hf_hrs cc_hrs cf_hrs el_hrs                                                                                              
        mig_p mig_c mig_i immig                                                                                                         
        vm_ch vm_sa vm_bl vm_fl vm_se lang aborig                                                                                       
        drivers transit wbike;                                                                                                          
        title 'New Vancouver Factorial Ecology';                                                                                        
run;