《多元统计分析》实验五.docVIP

实验五 心理健康状况与社会经济状况的对应分析 （王学民 编写） 一、实验目的 1.掌握如何使用SAS软件来进行对应分析； 2.看懂和理解SAS输出的结果，并学会以此来作出分析； 3.掌握对实际数据如何来进行对应分析。 二、实验内容 将由个人组成的样本按心理健康状况与社会经济状况进行交叉分类，分类结果见表1，对心理健康状况与社会经济状况进行对应分析。 表1 心理健康状况-社会经济状况数据 社会经济状况 A（高） B C D E（低） 心理健康状况 0（好） 121 57 72 36 21 1（轻微症状形成） 188 105 141 97 71 2（中等症状形成） 112 65 77 54 54 3（受损） 86 60 94 78 71 实验1 使用马赛克图（mosaic plot）展示了行轮廓、行密度以及列轮廓、列密度。 实验2 进行对应分析。 三、实验要求 1.用SAS软件完成对应分析的计算及作图； 2.根据SAS输出结果完成对应分析。 四、实验指导 1.使用马赛克图展示行轮廓、行密度以及列轮廓、列密度。 将表1中的数据按每一个人逐条输入，共1660个观测，有一个行变量和一个列变量，均设置为列名型的，保存为数据集sasuser.examp911。 在inshigt环境下打开数据集sasuser.examp911，见图1。选菜单过程如下： 在图1的界面中，分析盒形图/马赛克图（Y）（B）选rowX（见图2）；columnY输出作图3中的选择确定确定 随即出现如图4所示的马赛克图，展示了行轮廓和行密度。在图4中，对心理健康的每一种状况，A、B、C、D、E五个小方块的宽度显示了行轮廓，0、1、2、3四种心理健康状况的小方块高度显示了行密度。 图1 图2 图3 图4 在图1的界面中，分析盒形图/马赛克图（Y）（B）选columnX（见图5）；rowY 输出作图6中的选择确定确定 随即出现如图7所示的马赛克图，展示了列轮廓和列密度。在图7中，对社会经济的每一种状况，0、1、2、3四个小方块的高度显示了列轮廓，A、B、C、D、E五种社会经济状况的小方块宽度显示了列密度。 图5 图6 图7 2.进行对应分析 编制如下的SAS程序： data examp911; input mental $ A B C D E; cards; 0(好) 121 57 72 36 21 1(轻微症状形成) 188 105 141 97 71 2(中等症状形成) 112 65 77 54 54 3(受损) 86 60 94 78 71 ； proc corresp data=examp911 out=results rp cp short; var A B C D E; id mental; run; proc plot data=results; plot dim1*dim2=mental /box vspace=6 hspace=10 vaxis=-.30 to .30 by .15 haxis=-.30 to .30 by .15; run; 程序说明： “data examp911”是用data步产生一个数据集，并命名为examp911，这是一个临时数据集，属于临时库work； “input mental $ A B C D E”是input语句，用于从cards语句后面的数据区中读取数据，通过读取数据为变量命名、确定变量的类型等； “cards”指定本语句后面的是数据行； “proc corresp”是一个对应分析过程； “data= examp911”规定过程分析的是临时库中的examp911数据集； “out=results”指示生成包含行、列、附加的观测和附加的变量坐标的数据集； rp和cp分别指示打印行轮廓矩阵和列轮廓矩阵； short指示不打印除坐标之外的所有点和坐标统计量； “var A B C D E” 是一个var语句，指示列联表中的列是A、B、C、D、E； “id mental” 是一个id语句，指示用mental的值作为输出列联表中的这些行的标签； “proc plot” 是一个图形过程； “data=results” 规定作图的是临时库中的results数据集； “plot dim1*dim2=mental”是一个plot语句，指示dim1为垂直变量，dim2为水平变量，作图字符由变量mental给出； “box vspace=6 hspace=10 haxis=-.30 to .30 by .15 vaxis=-.30 to .30 by .15”是斜杠“/”后面的选项，box