第七章 方差分析和协方差分析.doc

第七章方差分析和协方差分析
(医学统计之星)
上次更新日期:
方差分析和协方差分析在SAS系统中由SAS/STAT模块来完成,其中我们常用的有ANOVA过程和GLM过程。前者运算速度较快,但功能较为有限;后者运算速度较慢,但功能强大,我们做协方差分析时就要用到GLM过程。本章将首先介绍方差分析所用数据集的建立技巧,然后重点介绍这两个程序步。
其实,这里的速度快慢只是相对而言,SAS的处理速度是首屈一指的。举个例子,这个暑假我做了一个有6600条记录的,7因素的,交叉设计的方差分析(是不是已经有人喊头痛了?),我先是用SPSS FOR WIN95 ,运行了大约10分钟才出结果。我又换用SAS FOR WIN95 ,结果用了――!
§ 方差分析数据集的建立技巧
方差分析的数据集格式
统计分析所用的数据格式和我们在分析整理资料时所用的格式是不同的。一般来说,数据集中应至少有一个结果变量,用于记录不同处理因素水平下观察值的大小;至少有一个处理因素变量,用于记录处理因素的类型及其水平数。以单因素方差分析为例,就应有一个结果变量和一个处理因素变量;而两因素的方差分析应有一个结果变量和两个处理因素变量。
某职业病防治院对31名石棉矿工中的石棉肺患者、可疑患者及非患者进行了用力肺活量测定,请给出数据集的结构(卫统p44 )。
解:数据集中应有两个变量,x和group。x记录肺活量的大小;group取值为1、2或3,分别代表石棉肺患者、可疑患者及非患者。
某厂医务室测定了10名氟作业工人工前、工中及工后4小时的尿氟浓度,请给出数据集的结构(卫统p46 )。
解:数据集中应有三个变量,x、group和worker。x记录尿氟浓度;group取值为1、2或3,分别代表工前、工中及工后;worker取值为1到10,分别代表10名工人。
方差分析数据集的建立技巧
可见方差分析的数据集其变量取值有一定的规律,因此可以利用循环语句和判断语句来简化输入。
。
解:我们可以一个一个的输入变量group的值,但这里给出一种更酷的建立方法。
libname a 'c:\user';
data ;
group=1;
if _n_>11 then
此处的_n_为SAS系统变量,表明读入的数据为第几条记录,由于石棉肺患者只有11人,因此第12条记录只能是可疑患者的记录,同理,从第20条记录开始的只能是非患者的记录。这里利用判断语句来完成从_n_到group的转换。
if _n_>20 then group=3;
else group=2;
;
input value @@;
cards;
... ...
;
proc print;
run;
。
解:此例中数据较有规律,各组的例数均相等,这可正是循环语句大显身手的时候。
libname a 'c:\user';
data ;
do group=1 to 3;
这里的SAS语句非常类似于Foxbase语句,两个循环分