第9章统计实验统计实验(对应分析).docVIP

实验五 对应分析 1．实验目的： 本实验讨论利用对应分析从众多变量和样品信息中找出变量间、样品间、变量与样品间的本质联系。通过该实验，能够起到如下的效果：(1) 理解对应分析的作用、思想、数学基础、方法和步骤；(2) 熟悉如何利用对应分析，提出问题、分析问题、解决问题、得出结论； （3）会调用SAS软件实现对应分析的各个步骤，根据计算的结果进行分析，得出正确的结论，解决实际的问题。 2．知识准备： 对应分析是从众多变量和样品信息中找出变量间、样品间、变量与样品间的本质联系。其思想是：对于某份数据(n份样品、p维数据)，其变量点（n维空间的点，坐标为该变量在各个样品处的值）的协差阵和样品点（p维空间的点，坐标为该样品在各个变量处的值）的协差阵有本质的联系，而且有相同的特征值，特征向量也具有某种联系。利用该联系进行适当的、尽量保留较多信息的降维，就会既反映变量间、样品间的本质联系，又反映变量与样品间的本质联系。对应分析的步骤大体分为：首先把指标进行正向化；然后计算过渡矩阵，消除原始数据量纲的影响，使样品和变量具有某种意义下的对等性，以便可以在同一坐标轴中进行描述；然后对数据进行R型因子分析，根据过渡矩阵的相关阵的特征根和累计贡献率选取适当的公因子，计算出R型（变量点对应的）因子载荷和Q型（样品点对应的）因子载荷；然后把样品点和变量点根据它们变换后的坐标（R型因子得分或Q型因子的得分），描述到同一坐标轴中；最后根据样品点和变量点间的距离进行分析，得出结论。 3．实验内容： 下面表１的数据是2009年广东省城镇居民家庭平均每人全年消费性支出构成的基本数据，其中的单位是百分比，数据来源于《广东省2009统计年鉴》： 表1 广东省城镇居民家庭平均每人全年消费性支出构成的基本数据 居　民　经 济　成　份 食品 衣着 居住 家庭设备用品 医疗保健 交通和通讯 教育文化娱乐服务 其他消费 最低收入户 52.84 3.64 14.97 4.15 5.53 10.36 6.24 2.27 困 难 户 54.52 3.43 15.24 4.65 5.55 8.63 5.70 2.28 低收 入户 52.77 4.78 13.03 4.56 4.67 10.29 7.57 2.33 中等偏下户 46.92 5.46 11.20 5.69 4.91 13.49 9.68 2.65 中等收入户 39.23 6.37 10.26 5.83 5.97 16.72 11.92 3.70 中等偏上户 35.55 6.59 9.65 6.36 5.56 18.77 13.46 4.06 高收 入户 29.56 6.58 12.08 7.04 5.08 20.35 15.39 3.92 最高收入户 28.70 7.17 9.86 7.05 5.91 21.10 15.71 4.50 利用对应分析对该数据进行处理，给出R型、Q型因子载荷，并结合该数据，给出适当的结论。 4．实验步骤： SAS程序: 读入数据： Data consumption; input type X1-X8; cards; 1 52.84 3.64 14.97 4.15 5.53 10.36 6.24 2.27 2 54.52 3.43 15.24 4.65 5.55 8.63 5.70 2.28 3 52.77 4.78 13.03 4.56 4.67 10.29 7.57 2.33 4 46.92 5.46 11.20 5.69 4.91 13.49 9.68 2.65 5 39.23 6.37 10.26 5.83 5.97 16.72 11.92 3.70 6 35.55 6.59 9.65 6.36 5.56 18.77 13.46 4.06 7 29.56 6.58 12.08 7.04 5.08 20.35 15.39 3.92 8 28.70 7.17 9.86 7.05 5.91 21.10 15.71 4.50 ; run; 进行对应分析，并画出散点图： Proc corresp data=consumption out=result; var X1-X8; id type; Proc plot data=result; plot dim1*dim2=*$type/ haxis=-0.06 to 0.1 by 0.02 vaxis=-0.35 to 0.35 by 0.1 vspace=3 hspace=10 HREF=0 VREF=0; run; 语句解释： “Proc corresp”指调用对应分析程序； “var X1-X8;”指变量是“X1-X8”； “id type;”指样品名是变量“type”；