/* Lesson 15-41 */
 /*    File Name = les1541.sas   07/24/03   */
/*
  protein.prn に対して
    主成分分析を行う。
    因子分析を行う。
*/

options linesize=72 pagesize=30;

data protein;
  infile 'protein.prn' firstobs=2;
  input country $15. Red White Eggs Milk Fish Cereals Starch Nuts Fr_Veg;
run;

proc print data=protein (obs=5);
run;

title '*** Principal Component Analysis for Protein Data ***';
proc princomp data=protein out=pscore;
  var Red--Fr_Veg;
run;
proc print data=pscore;
run;
proc plot data=pscore;                     
  plot prin1*prin2=country/vref=0.0 href=0.0;  
  plot prin1*prin3=country/vref=0.0 href=0.0;  
  plot prin1*prin4=country/vref=0.0 href=0.0;  
  plot prin2*prin3=country/vref=0.0 href=0.0;  
  plot prin2*prin4=country/vref=0.0 href=0.0;  
  plot prin3*prin4=country/vref=0.0 href=0.0;  
run;                                       

title '*** Factor Analysis for Protein Data ***';
proc factor data=protein;
  var Red--Fr_Veg;
run;

proc factor data=protein nfactors=4 out=fscore;
  var Red--Fr_Veg;
run;
proc print data=fscore;
run;
proc plot data=fscore;                     
  plot factor1*factor2=country/vref=0.0 href=0.0;  
  plot factor1*factor3=country/vref=0.0 href=0.0;  
  plot factor1*factor4=country/vref=0.0 href=0.0;  
  plot factor2*factor3=country/vref=0.0 href=0.0;  
  plot factor2*factor4=country/vref=0.0 href=0.0;  
  plot factor3*factor4=country/vref=0.0 href=0.0;  
run;                                       

title '*** Factor Analysis for Protein Data with Rotation ***';
proc factor data=protein nfactors=4 rotate=varimax out=fscore2;
  var Red--Fr_Veg;
run;
proc print data=fscore2;
run;
proc plot data=fscore2;
  plot factor1*factor2=country/vref=0.0 href=0.0;  
  plot factor1*factor3=country/vref=0.0 href=0.0;  
  plot factor1*factor4=country/vref=0.0 href=0.0;  
  plot factor2*factor3=country/vref=0.0 href=0.0;  
  plot factor2*factor4=country/vref=0.0 href=0.0;  
  plot factor3*factor4=country/vref=0.0 href=0.0;  
run;