实验八主成分分析.docVIP

实验八 主成分分析 一、实验目的和要求 能利用原始数据与相关矩阵、协主差矩阵作主成分分析,并能理解标准化变量主成分与原始数据主成分的联系与区别； 能根据SAS输出结果选出满足要求的几个主成分． 实验要求:编写程序，结果分析． 实验内容： 1.何为主成分？主成分的求法、性质？主成分个数选取的依据？ 主成分定义： 满足：（1）系数向量单位化； （2）各主成分不相关，无重叠信息，； （3）主成分方差由此递减． 步骤： 1）解，求的p个特征值； 2）对应的正交单位化的特征向量分别为； 其中 ， 3）的第个主成分为 ； 为个主成分构成的随机向量，则，其中 为正交矩阵， 且 故主成分向量的协方差阵为 各主成分的总方差 4）主成分的贡献率与累计贡献率 ——第k个主成分的贡献率 它反映了第k个主成分提取全部信息的多少 ——前k个主成分的累积贡献率 它反映了前k个主成分共同提取全部信息的多少． 性质： （1）系数向量单位化； （2）各主成分不相关，无重叠信息，； （3）主成分方差由此递减． 依据： 正交单位化特征向量 第K个主成分可以记为： 2.利用样品求主成分的步骤？ 1）求样本协方差矩阵特征值； 2）相应的正交单位化特征向量，； 3）第个样本主成分 ， 4） —— 第k个样本主成分的贡献率 ——前k个主成分的累积贡献率 5）第个样本主成分的个观测值称为第个样本主成分的得分，可以依据得分对各组样本观测数据进行排序： （） 样本主成分的观测数据（得分向量） 3.书上4.5 4.6 4.5 (1) 程序如下： data examp4_5 ; input id x1-x8; cards; 1 8.35 23.53 7.51 8.62 17.42 10.00 1.04 11.21 2 9.25 23.75 6.61 9.19 17.77 10.48 1.72 10.51 3 8.19 30.50 4.72 9.78 16.28 7.60 2.52 10.32 4 7.73 29.20 5.42 9.43 19.29 8.49 2.52 10.00 5 9.42 27.93 8.20 8.14 16.17 9.42 1.55 9.76 6 9.16 27.98 9.01 9.32 15.99 9.10 1.82 11.35 7 10.06 28.64 10.52 10.05 16.18 8.39 1.96 10.81 8 9.09 28.12 7.40 9.62 17.26 11.12 2.49 12.65 9 9.41 28.20 5.77 10.80 16.36 11.56 1.53 12.17 10 8.70 28.12 7.21 10.53 19.45 13.30 1.66 11.96 11 6.93 29.85 4.54 9.49 16.62 10.65 1.88 13.61 12 8.67 36.05 7.31 7.75 16.67 11.68 2.38 12.88 13 9.98 37.69 7.01 8.94 16.15 11.08 0.83 11.67 14 6.77 38.69 6.01 8.82 14.79 11.44 1.74 13.23 15 8.14 37.75 9.61 8.49 13.15 9.76 1.28 11.28 16 7.67 35.71 8.04 8.31 15.13 7.76 1.41 13.25 17 7.90 39.77 8.49 12.94 19.27 11.05 2.04 13.29 18 7.18 40.91 7.32 8.94 17.60 12.75 1.14 14.80 19 8.82 33.70 7.59 10.98 18.82 14.73 1.78 10.10 20 6.25 35.02 4.72 6.28 10.03 7.15 1.93 10.39 21 10.60 52.41 7.70 9.98 12.53 11.70 2.31 14.69 22 7.27 52.65 3.84 9.16 13.03 15.26 1.98 14.57 23 13.45 55.85 5.50 7.45 9.55 9.52 2.21 16.30 24 10.85 44.68 7.32 14.51 17.13 12.08 1.26 11.57 25 7.21 45.79 7.66 10.36 16.56 12.86 2.25 11.69 26 7.68 50.37