多元统计分析案例实验使用SAS软件对心理健康状况与社会经济状况的对应分析.doc

实验五 心理健康状况与社会经济状况的对应分析 一、实验目的 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.examp911mosaic（注：下载的本书数据集中有）。 在inshigt环境下打开数据集sasuser.examp911mosaic，见图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程序： proc corresp data=sasuser.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 haxis=-.30 to .30 by .15 vaxis=-.30 to .30 by .15; run; 程序说明： “proc corresp”是一个对应分析过程；“data= sasuser.examp911”规定过程分析的是sasuser库中的examp911数据集；“out=results”指示生成包含行、列、附加的观测和附加的变量坐标的数据集，并将其命名为results，属work（临时）库；选项“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 vaxis=-.30 to .30 by .15 haxis=-.30 to .30 by .15”是斜杠“/”后面的选项，box要求画出的边框围住整个图形，vspace=6规定垂直轴上刻度标记之间输出的行数目是6，hspace=10规定水平轴上刻度标记之间输出的列数目是10，vaxis=-.30 to .30 by .15定义垂直轴上的刻度从-0.30至0.30，增量为0.15，haxis=-.30 to .30 by .15定义水平轴上的刻度从-0.30至0.30，增量为0.15。 程序运行后输出如下结果： 上述对应分析图显示，心理健康状况1和2具有类似的轮廓（两个行点彼此靠近），0和3的轮廓非常不同（两个行点远离）；社会经济状况B和C的轮廓相似（两个列点靠得近），A和E的轮廓很不相似（