SAS主成分分析的实现

壹壹九

/*通过MSA做KMO检验， 在分析结果中看MSA值，大于０.６则适合作主成分分析*/

procfactor data=grade  rotate=varimax msa ;

varg1-g5;

run;

procprincomp data=grade out=graderesult outstat=ssss ; /*主成分分析过程*/

varg1-g5 ; /*定义进行主成分分析的变量*/

run;

/*从主成份的outstat指定的数据集中提取特征根pc1 pc2*/

datazhuchenfen;

setssss;

if_TYPE_="EIGENVAL" then do;

        call symput("pc1",g1);

        call symput("pc2",g2);

end;

run;

/*计算综合得分*/

datazonghe;

setgraderesult;

zonghe= Prin1 * &pc1 + Prin2 * &pc2;

run;

/*以综合得分排名*/

procsort data=zonghe;

byzonghe;

run;