********************************************; 
title 'example of medical statics 2003 5';
********************************************;  
data adata;                                    
infile '/st5/eisei4/sasuser/data?.txt';         
  input (no area group total var1-var12) (5.);
********************************************;  
  rate98=var1/var9*100;                      
  rate99=(var1+var2+var3-var4);            
********************************************;  
  area2=1;                                     
  if area>=4 then area2=2;                     
********************************************;  
 run;                                          
********************************************;  
proc means;run;                                
proc plot;plot no*total;run;                   
********************************************;  
proc ttest;class group;var total;run;
********************************************;  
***　以下を実行するときは第１カラムの「*」を除くこと　*************;  
********************************************;  
*proc means;by area;var var1;run;             
********************************************;  
*proc glm;                                      
 *class area;
 *model var1=area/ss1;
 *means area/tukey;run;
********************************************;  
*proc glm; 
 *class group area;
 *model total=group area group*area/ss2;                  
 *means group*area/tukey;
 *run;                       
********************************************;  
*proc reg;
 *model total=var1 rate98 rate99/stb;         
 *run;
********************************************;  
*proc factor method=prin nfact=3 rotate=varimax plot;
 *var var1-var7;                                
 *run; 
********************************************;