主成分分析聚类分析VIP免费

1 随机抽取管理学院10 名学生，对其4 门课程的考试成绩进行统计，如下表所示，这 4门课程分别为多元统计分析1X ，运筹学2X，经济学3X，管理学4X. 使用主成分分析方法对学生成绩进行分析. 学生编号 多元统计分析 运筹学 经济学 管理学 1 77 82 67 81 2 63 78 80 81 3 75 73 71 81 4 55 72 63 68 5 31 55 60 73 6 67 81 82 67 7 70 81 78 80 8 66 81 73 71 9 70 68 72 63 10 57 73 55 60 >> x1=[77 63 75 55 31 67 70 66 70 57]; >> x2=[82 78 73 72 55 81 81 81 68 73]; >> x3=[67 80 71 63 60 82 78 73 72 55]; >> x4=[81 81 81 68 73 67 80 71 63 60]; >> data=[x1;x2;x3;x4]'; %输入观测值数据矩阵 >> [n,m]=size(data); >> for i=1:m %将数据矩阵中心标准化 sddata(:,i)=(data(:,i)-mean(data(:,i)))./std(data(:,i),1); end >> [P,score,egenvalue,t2]=princomp(sddata) %做主成分分析 P = -0.5511 0.3268 -0.3624 0.6769 -0.5588 0.3358 -0.2089 -0.7289 -0.5110 -0.1193 0.8460 0.0944 -0.3505 -0.8753 -0.3307 -0.0398 score = -1.3489 -0.2567 -1.2840 -0.0315 -1.2458 -0.9740 0.5341 -0.2698 -0.8704 -0.7467 -0.5863 0.7336 1.1642 0.3096 -0.2225 -0.2707 3.3634 -1.5690 0.3950 -0.0572 -1.1054 0.8480 1.1534 -0.2345 -1.5954 -0.5212 0.0958 -0.1867 -0.6992 0.4872 0.1005 -0.4105 0.4738 0.9799 0.5768 1.0306 2 1.8637 1.4430 -0.7628 -0.3031 egenv alu e = 2.7502 0.9334 0.5275 0.2334 t2 = 3.8620 2.4333 3.8297 1.0034 7.0609 3.9724 1.3834 1.1732 6.2912 4.9906 >> for k=1:m gx l(k)=su m(egenv alu e(1:k))/su m(egenv alu e); end >> gx l %输出累计贡献率 gx l = 0.6188 0.8288 0.9475 1.0000 >> plot(score(:,1),score(:,2),'r+') %画出第一第二主成分的散点图 >> gname 第一主成分43211XXXXY0.3505-0.5110-0.5588--0.5511，所有科目考试成绩的系数均为负，且差异不大，故1Y 可解释为学生的综合学习成绩，该主成分得分越小（散点图中的位置越靠左），综合成绩越好. 第二主成分43212XXXXY0.8753-0.1193-0.33580.3268，数学科目考试成绩的系数均为正，专业科目考...