第七章方差分析和协方差分析（精品）

《第七章方差分析和协方差分析（精品）》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、第七章方差分析和协方差分析（医学统计之星）上次更新日期：方差分析和协方差分析在SAS系统屮由SAS/STAT模块来完成，其屮我们常用的有ANOVA过程和GLM过程。而者运算速度较快，但功能较为有限；后者运算速度较慢，但功能强大，我们做协方差分析时就要用到GUV1过程。本章将首先介绍方差分析所用数据集的建立技巧，然后重点介绍这两个程序步。Q3其实，这里的速度快慢只是相对而言，SAS的处理速度是首屈一指的。举个例子，这个暑假我做了一个有6600条记录的，7因素的，交叉设计的方差分析（是不是已经有人喊头痛了？），我先是用SPSSF

2、ORWIN957.5来做，运行了大约10分钟才出结果。我又换用SASFORWIN956.12来做结果用了一一2.47秒!§7.1方差分析数据集的建立技巧7.1.1方差分析的数据集格式统计分析所用的数据格式和我们在分析整理资料时所用的格式是不同的。般来说，数据集中应至少有一个结果变量，用于记录不同处理因素水平下观察值的大小；至少有一个处理因素变量，用于记录处理因素的类型及其水平数。以单因素方差分析为例，就应有一个结果变量和一个处理因索变量；而两因索的方差分析应有一个结果变量和两个处理因素变量。例7.1某职业病防治院对31名石棉

3、矿工中的石棉肺患者、可疑患者及非患者进行了用力肺活量测定，请给出数据集的结构（卫统P44例5・1）。解：数据集中应有两个变量，x和groupox记录肺活量的大小；group取值为1、2或3,分别代表石棉肺患者、可疑患者及非患者。例7.2某厂医务室测定了10名氟作业工人工前、工中及工后4小时的尿氟浓度,请给出数据集的结构（卫统p46例5.2）。解：数据集中应有三个变量，x、group和*70厂1<8厂。x记录尿氟浓度；group取值为1、2或3,分别代表工前、工中及工后；worker取值为1到10,分别代表10名人。7.1.

4、2方差分析数据集的建立技巧可见方弄分析的数据集其变量取值有一定的规律，因此可以利用循环语句和判断语句来简化输入。例7.3请建立例6.1的数据集。libnamea'c:user';dataa.wtli5_l;group=l;if_n_>llthenif_n_>20thengroup=3;elsegroup=2;解：我们可以一个一个的输入变量group的值，但这里给出一种更酷的建立方法。此处的_n_为SAS系统变量，表明读入的数据为第几条记录，由于石棉肺患者只有X人，因此第12条记录只能是可疑患者的记录，同理，从第20条记录开

5、始的只能是非患者的记录。这里利用判断语句来完成从_n_到group的转换。inputvaluecards;1.81.41.52.13.33.5procprint;run;例7.4请建立例7.2的数据集。解：此例中数据较有规律，各组的例数均相等，这可正是循环语句大显身手的时候。libnameaBc:user1;dataa.wtli5_2;dogroup=lto3;这里的SAS语句非常类似于Foxbase语句.两个循环分别控doworker=lto10;制了group和worker的取值，这在输入数据时是常见的一种技巧0inp

6、utvalueoutput;两个变量只是循环变量，因此要用output语句写入数据集中.end;end;cards;90.5388.4347.37・・・请注意X变童和处理因素变量的对应关系，不要改变数据块的次序，・・・105.2758.95■否则建立的数据集是错的.procprint;run;§7.2引例例7・5请回答在例7.1中三组石棉矿工的用力肺活量有无差别？解：由于是三组做比较，因此要用方差分析法。在建立数据集后可用ANOVAi±程或GLM过程进行分析，同时给出选项进行各组均数的两两比较（此处用SNK法），两过程的输出

7、格式基本一致。对于单因素的方差分析，在比较之前一般要考虑方差齐性的问题。这里可以用GLM过程中的HOVTEST选项实现，实际上这也是在SAS中做方差齐性检验的唯一途径。①设定数据库环境:LIBNAMEA‘C:USEFT;②数据步，建立数据集（略）：③ANOVA或GUV1过程，进行方差分析并同时做两两比较（SNK法）：PROCANOVADATA=A.WTLI5_1;PROCGLMDATA=A.WTLI5_1;CLASSGROUP;CLASSGROUP;MODELVALUE=GROUP;MODELVALUE=GROUP/SS1

8、;MEANGROUP/SNK;MEANGROUP/SNKHOVTEST：RUN;RUN;3这里的GLM过程中有两个选项，SS1要求只输出方差分析SS1模型的结果（默认输出SS1和SS3两种模型的结果），HOVTEST则要求进行方差齐性检验（该选项只能在SAS6.12以上的版本中实现）。§7