1、第七章_方差分析和协方差分析第七章 方差分析和协方差分析(医学统计之星)上次更新日期: 方差分析和协方差分析在 SAS 系统中由 SAS/STAT 模块来完成,其中我们常用的有 ANOVA 过程和 GLM 过程。前者运算速度较快,但功能较为有限;后者运算速度较慢,但功能强大,我们做协方差分析时就要用到 GLM 过程。本章将首先介绍方差分析所用数据集的建立技巧,然后重点介绍这两个程序步。其实,这里的速度快慢只是相对而言,SAS 的处理速度是首屈一指的。举个例子,这个暑假我做了一个有 6600 条记录的,97.5 来做,运行了大约 10 分钟才出结果。我又换用 SAS FOR WIN95 6.12
2、 来做,结果用了2.47.1 方差分析数据集的建立技巧7.1 某职业病防治院对 31 名石棉矿工中的石棉肺患者、可疑患者及非患者进行了用力肺活量测定,请给出数据集的结构(卫统 p44 例 5.1) 。解:数据集中应有两个变量,x 和 group。x 记录肺活量的大小;group 取值为 1、2 或 3,分别代表石棉肺患者、可疑患者及非患者。例7.1.2 方差分析数据集的建立技巧可见方差分析的数据集其变量取值有一定的规律,因此可以利用循环语句和判断语句来简化输入。例7.4 请建立例7.37 58.95请注意 x 变量和处理因素变量的对应关系,不要改变数据块的次序,否则建立的数据集是错的。;pro
3、c print;run;7.5 请回答在例7.3 ANOVA 过程ANOVA 过程是 SAS 系统中用于方差分析的几个程序步之一,它适用于平衡的数据。这意味着如果按所分析的变量将观测值分类,则每个格子内的观测值数目应相同。如果数据不满足这一要求,则 ANOVA 过程的分析结果并不可靠,分析者应改用 GLM 过程来处理数据。7.3.2 语法说明Means 语句的选项主要用来指定两两比较的方法和检验水准,主要有:两两比较方法 BON(Bonferroni 检验) 、DUNCAN(新复极差法) 、DUNNETT、GABRIEL、LSD(最小显著差法) 、SCHEFFE、SIDAK 、SNK(q 检验
4、) 、TUKEY 、WALLER 其它 ALPHA= p 值 更改水准7.2 的数据为例,ANOVA 过程的输出结果分两大部分,第一部分是分析变量的列表描述,第二部分是统计结果,先是所用方差分析模型的检验,然后是各个处理因素和交互作用的检验结果,最后则是两两比较的结果。Analysis of Variance Procedure 方差分析过程Class Level Information 处理因素取值情况Class Levels Values处理因素变量名 取值数 具体值GROUP 3 1 2 3WORKER 10 1 2 3 4 5 6 7895.8728667066242 4.51 0.0024Error 18 175338111Corrected Total 29 65261.43814667R-Square C.V. Root MSE 应变量名 Mean所用模型的决定系数 变异系数 剩余标准差 应变量均数