sas第九章 t检验和方差分析

第九章 t 检验和方差分析

在科研中，我们往往是根据样本之间的差异，去推断其总体之间是否有差异。样本差异可能是由抽样误差所致，也可能是由本质的不同所致。应用统计学方法来处理这类问题，称为“差异的显著性检验”。若已知总体为正态分布，进行差异的显著性检验，称为“参数性检验”，SAS 中MEANS 、TTEST 、ANOVA 、GLM 等均属此类检验；若未知总体分布，进行差异的显著性检验，称为“非参数性检验”，SAS 中采用NPAR1WAY 过程。

第一节 t 检验

9.1.1 简介

t 检验是用于两组数据均值间差异的显著性检验。它常用于以下场合：

1．样本均值与总体(理论)均值差别的显著性检验

检验所测得的一组连续资料是否抽样于均值已知的总体

根据大量调查的结果或以往的经验，可得到某事物的平均数(例如生理生化的正常值)，以此作总体均值看待。

SAS 中采用MEANS 过程，计算出观察与总体均值的差值，再对该差值的均值进行t 检验。

2.同一批对象实验前后差异的显著性检验(自身对照比较)或配对资料差异的显著性检验(配对比较检验)

比如，在医学研究中，我们常常对同一批病人治疗前后的某些生理生化指标(如血压、体温等)进行测量，以观察疗效；或对同一批人群进行预防接种，以观察预防效果；或把实验对象配成对进行测定，比较其实验结果。

SAS 中采用MEANS 过程，计算出两样本观察的差值(如治疗前、后实验数据的差值)，再对该差值的均值进行t 检验。

3．两样本均值差异的显著性检验

作两样本均值差异比较的两组原始资料各自独立，没有成对关系。两组样本所包含的个数可以相等，也可以不相等。每组观测值都是来自正态总体的样本。 设1X 与2X 为两样本的均值，1n 与2n 为两样本数，21s ，22s 为两样本方差，分两种情形，其数学模型为：

(1)方差齐(相等)时：

)/1/1(2122

1n n s x x t +-=

)2/(])1()1[(212222112-+-+-=n n s n s n s

(2)方差不齐时： 2221212

1//n s n s x x t +-=

SAS 中采用TTEST 过程，先作方差齐性检验(F 检验)，然后根据方差齐(EQUAL)和方差不齐(UNEQUAL)输出t 值和P 值以及基本统计量。

在作方差齐性检验时，用F 检验。F 值计算公式为：

)

,(),(22212221S S Min S S Max F =

9.1.2 用MEANS 过程作t 检验

1．过程格式

PROC MEANS MEAN STD STDERR T PRT ；

VAR 变量表；

2．说明

(1)PROC MEANS 语句中，选择了5个统计量：均值、标准差、标准误差、t 值、P 值。

(2)VAR 语句中的变量是分析变量。缺省时，计算所有数值型变量。

3．举例

例1: 样本均值与总体均值差别的显著性检验。

已知某水样中含CaCO 3的真值为20.7mg ／L ，现用某方法重复测定该

水样11次，CaCO 3的含量为：20.99，20.41，20.10，20.00，20.91，22.60，

20.99，20.41，20，23，22。问用该法测CaCO 3的含量所得的均值与真值

有无显著差别 ?

程序： (yp111.sas)

编程说明：

在数据步中，变量x 读取测定值，产生一个差值变量y(y=x-20.7)，在过程步中，计算出Y 的均值、标准差、标准误差、t 值、P 值。

结果说明：

因t ＝1.0636907，0.05

的均值与真值间无显著差异。此法可信。

例2：配对比较的t 检验。

研究食物中维生素E 与肝脏中维生素A 含量的关系。将大白鼠按性别、体重配对。每对随机分配，一个用正常饲料，一个用缺乏维生素E 的饲料。经过一个时期饲养，杀死动物测定肝中维生素A 的含量，结果如下表：

大白鼠肝脏中维生素A含量（IU/g）

配对号 1 2 3 4 5 6 7 8

正常饲料3550 2000 3000 3950 3800 3750 3450 3050 缺乏E饲料2450 2400 1800 3200 3250 2700 2500 1750

程序：（yp112.sas）

编程说明：

数据步中，把每对数据中的一个作为x，另一个作为y，计算出差数d(d=x-y)，在过程步中，计算差数d的均值、标准差、标准误差、t值、P值。

结果说明：

因t=4.21 p=0.0040<0.05，故有非常显著差异，即正常饲料组鼠肝维生素A含量比维生素E缺乏组的含量大。

例3 自身对照比较的t检验。

应用克矽平治疗矽肺患者10名，治疗前后血红蛋白的含量如下表，问该药是否会引起血红蛋白的变化?

治疗前后血红蛋白的含量（mg%）

治疗前11.3 15.0 15.0 13.5 12.8 10.0 11.0 12.0 13.0 12.3 治疗后14.0 13.8 14.0 13.5 13.5 12.0 14.7 11.4 13.8 12.0

程序：（yp113.sas）

编程说明：

数据步中，把每对数据中的一个作为x，另一个作为y，计算出差数d(d=x-y)，在过程步中，计算差数d的均值、标准差、标准误差、t值、P值。

结果说明：

因t=1.1989377，O.05

9.1.3用TTEST过程作t检验

1．过程格式

PROC TTEST [DATA=数据集]；

CLASS变量；

VAR变量表；

2．说明

(1)CLASS语句中的变量是分类变量，其水平值只能有两个，并对应两组观察。是必选语句。

(2)VAR语句中的变量是被分析的变量，如果缺省，则对所有的数值型变量进行分析。

3．举例

观察某药物对大白鼠肉瘤的影响。数据如下表，试作差异性检验。

对照组56 55 54 53 56 52 57 54 52 56

实验组50 48 49 49 50 50 60 55 43 52 56 57 程序： (yp114.sas)

编程说明：

数据步中，用循环控制变量a作分类变量，其水平值为1(对照组)和2(实验组)，变量n表示样本数，分别为10和12，用变量x读取原始数据，在过程步中，用CLASS语句标识分类变量，用VAR语句标识分析变量。

结果说明：

输出的最后：H0：Variances are equal(方差相等)表示给出的是对方差相等假设的结果。

先看方差齐性检验(F检验)，然后根据F检验的结果，选择方差齐(Equal)或不齐(Unequal)的t值和p值，以及两组观察的均值、标准差、标准误差、最大值、最小值。

找Prob>F’=右边的值，即方差相等假设检验的P值，如果此值大于0.01则在方差相等假设下继续进行，否则只能使用近似T检验。因P=0.0074，说明方差差异显著，即方差不齐。

当方差不齐时，T=2.0000，0.05

第二节方差分析

当试验结果受到多个因素的影响，而且也受到每个因素的各水平的影响，为从数量上反映各因素以至各因素诸水平对试验结果的影响时使用方差分析的方法。

方差分析的基本思想是把全部数据关于总均值的离差平方和分解成几个部分，每一部分表示某因素交互作用所产生的效应，将各部分均方与误差均方相比较，从而确认或否认某些因素或交互作用的重要性。用公式概括为：

总变异＝组间变异+组内变异

其中：组间变异由各因素所引起，组内变异由个体差异所引起的，或者说由误差引起的。

常用的方差分析法有以下4种：

（1）完全随机设计资料的方差分析(单因素方差分析)

（2）随机区组设计资料的方差分析(二因素方差分析)

（3）拉丁方设计资料的方差分析(三因素方差分析)

（4）R*C析因设计资料的方差分析(有交互因素的方差分析)

SAS系统中，ANOVA过程可以处理以上情形的方差分析，但它要求每个分类因子的组合观察数相等，即数据是均衡的。若不均衡，就要求用GLM过程进行处理。在只考虑组间变异和误差变异时，称为单向方差分析。此时ANOVA会自动处理均衡和非均衡数据。

在方差分析中，每次只研究1个指标时，称之为一元方差分析(简称ANOVA)，同时考虑多个指标时，称之为多元方差分析(MANOVA)。在这一节里，我们还将讨论协方差分析。

9.2.1 均衡数据的方差分析(ANOVA过程)

1．过程格式

PROC ANOVA 选项

CLASS 变量表；

MODEL 因变量表=效应；

MEANS 效应[／选择项]；

2．使用说明

(1)程序中，CLASS语句和MODEL语句是必需的，而且，CLASS语句必须出现在MODEL语句之前。

(2)CLASS语句中的变量是分类变量，可以是数值型，也可以是字符型。

(3)MODEL语句指明因变量和自变量(因子变量)效应。效应是分类变量的各种组合，效应可以是主效应、交互效应、嵌套效应和混合效应。对应的效应模型如下：

·主效应模型MODEL y=a b c；

模型中，a，b，c是主效应，y是因变量。下同。

·交互模型MODEL y=a b c a*b a*c b*c a*b*c；

模型中，a*b，a*c,b*c，a*b*c是交互效应。

·嵌套效应模型MODEL y=a b c(a b)；

模型中，c(a b)是嵌套效应。

·混合效应模型MODEL y=a b(a) c(a) b*c(a)；

(4)MEANS语句是选择语句，计算并输出所列的效应对应的因变量均值，若指明了选择项，则将进行主效应均值间的检验。常用的选择项如下：BON、DUNCAN、LSD、REGWF、REGWO、SNK(Q检验)、SCHEFFE、SIDAK、SMM(GT2)、TUKEY、WALLER。以上选择项在实际应用中，一般选择一种或两种方法即可。

ALPHA=p确定检验的显著性水平。缺省值是0.05。

3．举例

(1)完全随机设计资料的方差分析(单因素方差分析)

某劳动卫生研究所研究棉布、府绸、的确凉、尼龙四种衣料吸附十硼氢量。每种衣料各做五次测量，所得数据如下表。试检验各种衣料吸附十硼氢量有没有显著差别?

各种衣料间棉花吸附十硼氢量

棉布府绸的确凉尼龙

2.33 2.48

3.06

4.00

2.00 2.34

3.06 5.13

2.93 2.68

3.00

4.61

2.73 2.34 2.66 2.80

2.33 2.22

3.06 3.60

程序： (yp115.sas)

编程说明：数据步中，用循环控制变量a做分类变量，其水平数是4，分别代表不同的衣料。过程步中，用CLASS语句指明一个因素a，用MODEL 语句反映出该因素的效果模型。

结果说明：

在输出中，找CLASS语句指出的变量的P值。此例中，P≤0.0003，可得出各衣料组间有非常显著差异。说明各种衣料间吸附十硼氢量是不同的。R-Square(R平方)对单向方差分析时，描述组间变异占总变异的比例，它越接近1，说明变异越归因于组间变异。

(2)随机区组设计资料的方差分析(两因素方差分析)

用4种不同方法治疗8名患者，其血浆凝固时间的资料如下表，试分析影响血浆凝固的因素。

程序：(yp116.sas)

编程说明：

数据步中，变量a代表不同治疗方法，其水平数是4，变量b代表区

组因素不同患者，其水平数是8。过程步中，用CLASS语句指明两个因素a和b，用MODEL语句指明二因素的效果模型。

结果说明：

总处理间：F＝14.03，P≤O.0001，故总体有非常显著差异。

A因素：F=6.62，P≤O.0025，故认为因素A(治疗方法)对血浆凝固时间影响很大。

B因素：F=17.20，P≤0.0001，故认为因素B(不同患者)对血浆凝固时间影响很大。

如果想对4种治疗方法进行两两比较，可在MODEL语句后增加一MEANS语句。例如，用SNK法比较，语句为：

MEANS A／SNK；(yp117.sas)

(3)拉丁方设计资料的方差分析(三因素方差分析)

若试验中涉及到3个因素，各因素间不存在交互作用，或交互作用很小可忽略不计，同时各因素的水平数又都相同，此种资料称拉丁方设计资料。

例：五种防护服，由五人各在不同的五天中穿着测定脉博数，如下表所示。试比较五种防护服对脉搏数有无不同。

程序：(yp118.sas)

编程说明：

数据步中，变量a代表试验时间因素，变量b代表受试者因素，变量C代表防护服因素。过程步中，用CLASS指明三个因素，MODEL语句反映三个因素的效果模型。

结果说明：

因F=6.80 P≤0.0011，故总体有非常显著差异，其中只有一个区组B 因素(受试者)影响极大(F=16.27 P≤0.001)，而区组A因素(试验时间)和处理因素(防护服)都无显著差异。

(4)R*C析因设计资料的方差分析(有交互因素的方差分析)

治疗缺铁性贫血病人12例，分为4组给予不同治疗，一个月后观察红细胞增加数(百万／mm)，资料如下表。试分析两种药物对红细胞增加的影响。

程序：(yp119.sas)

编程说明：

数据步中。变量a 代表甲药因素．变量b 代表乙药因素。过程步中，用CLASS 语句指明两个因素，用MODEL 语句反映主效应(a b)和交互效应(a*b)模型。

结果说明：

因F=98.75， P=0.0001<0.01，故总体有非常显著差异，其中A 因素(甲药)、B 因素(乙药)、a*b 因素(a 和b 的交互作用)都对红细胞增加数有非常大的影响。

为进一步分析其影响的强弱，可用MEANS 语句：

MEANS a b a*b ： (yp120.sas)

(5)用Dunnet 法进行均值间的比较(对照组与多个实验组的比较检验)。

Dunnett 方法在比较若干“实验水平”与“对照水平”时使用一个临界差距（T ）：

n

ms q T 2=，其中ms 为将方差分析表中交互效应归入残差平方和后的均方，即：

df df ss

ss ms 残差交互效应残差交互效应++=

n 为样本大小, q 为来自Dunnett 分布的分位数，可根据样本组个数r ，自由度df(交互效应df+残差df)，以及比较的显著水平α，查表确定。

Dunnett 多重比较仅适用各“实验水平”与“对照水平”之间差异的比较，不适用于各“实验水平”间的相互差异比较。

例：有资料如下表，试比较白血病鼠与正常鼠脾中DNA平均含量(mg ／g)差别有无显著意义?

表；白血病鼠与正常鼠脾中DNA平均含量(mg／g)

脾中DNA平均含量(mg／g)

A 正常鼠（对照组）12.3 13.2 13.7 15.2 15.4 15.8 16.9

B1 自发性白血病鼠10.8 11.6 12.3 12.7 13.5 13.5 14.8

B2 移植性白血病鼠9.8 10.3 11.1 11.7 11.7 12.0 12.3 程序：（yp121.sas）

编程说明：

数据步中，变量a代表因素。过程步中，用MODEL语句指明一因素的效果模型，用MEANS语句进行a因素水平间的Dunnett法比较。

结果说明：

因F=11.42，P≤0.0006，故总体或A因素各水平间有非常显著差异。当A因素有意义时，可进一步分析组间的差异。

Dunnet法检验时，用“***”指明有显著性，由此可以看出，两组均值(A=2，A=3)与对照组(A=1)相比较，差别有显著意义。若改换对照组，假设是第2组(A=2)。应修改选择项DUNNETT(‘2’)。输出中给出DUNNET t的临界值为2.399、A3与A1均值之差为-1.9，其95％置信限的下限和上限分别为-3.597和-O.203。

(6)用SNK法(Q检验)作均值间的两两比较。

为研究酵解作用对血糖浓度的影响，从8名健康人中抽取了血液并制备成血滤液。每一个受试者的血滤液又分成4份，然后随机地把各份血滤液分别放置0、45、90、135分钟后测定其中血糖浓度(资料见下表)。试比较放置时间对血糖浓度有无影响。

程序：(yp122.sas)

编程说明：

数据步中，变量a代表受试者因素，变量b代表放置时间因素。过程步中，作二因素的方差分析，用CLASS指明二个因素a和b，用MODEL语

句指明二因素的效果模型，用MEANS语句对因素b作SNK检验(Q检验)。

结果说明：

因F=43.66,P≤0.0001，故总体有非常显著差异。

因F=28.74, P≤O.0001，说明因素A(受试者)对血糖浓度有很大影响。

因F=78.49,P≤0.000l，说明因素B(放置时间)对血糖浓度有很大影响。

对B因素各组均值作SNK法比较检验:

经SNK检验，可以看出：由于O和45分钟的均值间用同一字母A分组，故无显著差别(Means with the same letter are not significantly different)，其余都有显著差别。

9.2.2.非平衡数据的方差分析(GLM过程)

在SAS／STAT中GLM(General Linear Models)过程分析功能最多。它包括：

简单回归(一元回归)

加权回归

多重回归及多元回归

多项式回归

方差分析(尤其对不平衡设计资料更为有效)

偏相关分析

协方差分析

多元方差分析

反应面模型分析

重复测量方差分析

这里，我们先介绍有关GLM过程在方差分析中的应用。

GLM过程里，采用MODEL语句反映因变量与自变量(即效应)的模型，其形式是：

模型说明模型类型

MODEL Y=A B C 主效应

MODEL Y=A B A*B 交互效应

MODEL Y=A B A(B) 嵌套效应

MODEL Y1 Y2=A B 多元方差分析

MODEL Y=A X 协方差分析

其中A，B，C是分类变量；X，Y是连续型变量。

前面介绍的ANOVA过程只能用于平衡设计资料(每个分类因子的组合观察数相等)的方差分析，当不平衡时，只能用GLM过程。

1．过程格式

PROC GLM；

CLASS 变量表；

MOIDEL 因变量一效应；

MEANS 效应／[选择项]；

RUN；

2．使用说明

(1)程序中，CLASS语句和MODEL语句是必需的，而且CLASS语句必须出现在MODEL语句之前。

(2)CLASS语句中的变量是分类变量，可以是数值型，也可以是字符型。

(3)MODEL语句指明因变量和自变量(因子变量)效应，其效应可以是主效应、交互效应、嵌套效应和混合效应。对应的效应模型同ANOVA过程中效应模式相同。

(4)MEANS语句是选择语句，计算所列的效应对应的因变量均值，若指明了选择项，则将进行主效应均值间的检验。常用的选择项如下：BON、DUNCAN、LSD、REGWF、REGWQ、SNK(Q检验)、SCHEF、SIDAK、SMM (GT2)、TUKEY、WALLER、DUNNETT(1-n检验)

以上选择项在实际应用中，一般选择一种或两种方法即可。

ALPHA=p给出检验的显著性水平。缺省值是0.05。

3．举例

(1)单因素方差分析

健康男子各年龄组淋巴细胞转化率(％)如下表，问各组淋巴细胞转化率的均值之间的差异是否显著?

表：健康男子各年龄组淋巴细胞转化率(％)

11～20岁58 61 61 62 63 68 70 70 74 78

41～50岁54 57 57 58 60 60 63 64 66

61～75岁43 52 55 56 60

程序：(yp123.sas)

编程说明：

变量a代表年龄组因子，有3个水平值；变量n代表各组例数。过程步中，用CLASS语句指明因素为分类变量a，用MODEL语句指明单因素的效果模型。

结果说明：

因F=9.77, P=0.001<0.05，故总体或A因素各水平间有非常显著差异。即说明各年龄组的淋巴细胞转化率有显著差异。

（2）二因素方差方析

程序：（yp124.sas）

编程说明：

数据步中，变量a代表因素a，其水平数是2，变量b代表因素b，其水平数是3，变量n是每组的重复次数。过程步中，用CLASS语句指明两个因素a和b，用MODEL语句指明带有交互影响二因素的效果模型。

结果说明：

通常GLM过程按两种方法输出各因素所引起的离差平方和(Type I SS，Type Ⅲ SS)。Type I SS按累积效应(有交互效应)输出，Type Ⅲ SS 按独立效应(只有主效应)输出。即若有充分的理由将所有因素按主次顺序(交互效应在主效应之后)排列在MODEL语句中，应选择Type I SS的结果。否则选择Type Ⅲ SS的结果。这里我们选择Type I SS的输出。

讨论：因F=2.50，O.05

结论：因F＝4.22，P≤0.0127，故总体有显著差异，而且A因素、B 因素的F值和P值分别为F=4.70,P≤0.0378，F=3.98,P≤O.0287，故认为两因素对分析指标都有显著影响。

9.2.3协方差分析

协方差分析法是把直线回归法与方差分析法结合起来的一种方法。它利用回归的关系消除自变量不等的影响后，再进行方差分析。如在营养实验中，用直线回归把食物消耗量(或初始重量)与所增体重的关系找出来，然后求出当食物消耗量(或初始重量)化为相等时，各饲料组动物所增体重的调整均值，用协方差分析检验调整均值间的差异。协方差分析适合于完全随机化设计资料、随机化区组设计资料、拉丁方资料、析因设计资料等。 1．过程格式

PROC GLM；

CLASS 效应；

MODEL 因变量＝协变量表效应／SOLUTION；

LSMEANS 效应[／选择项]；

RUN；

2．使用说明

(1)CLASS语句指明分类变量(因素)，也叫方差变量。

(2)MODEL语句指明因变量与协变量、分类变量的模型。选择项SOLUTION要求输出回归方程中各参数的估计值。

(3)LSMEANS语句是可选语句。要求输出每个列在该语句中效应的最小二乘方均值。常用的选择项有：

E 计算最小平方均值的可估计函数

STDERR 输出最小平方均值标准误差

PDIFF 输出假设HO

E=effects 指定模型中的一个效应项用作误差项

SINGULAR=n 调整可估计性检验

3．举例

实验中将15只白鼠用随机抽样法分为两组，第一组10只白鼠饲以人工母乳，第二组5只白鼠饲以奶粉。现将两组白鼠九周内之食物消耗量(X，克)及所增体重(Y，克)列表下，问：摄取两种不同饲料的白鼠所增体重的均值有无显著的差别?

程序： (yp126.sas)

编程说明：

数据步中，变量a代表组别，是离散变量，其水平数是2；变量X是协变量，用来读取食物消耗量的数据；变量Y是因变量，用来读取所增体重的数据，变量n是每组的样本数。过程步中，用CLASS语句指明分类变量，用MODEL语句指明因变量与协变量的模型。

结果说明：

因F＝0.08, 0.05

耗量作了调整后所增体重无显著差异。

根据Type皿SS对应的结果，可以看出：因F=67.39, P≤0.0001，故认为x对y的影响非常显著，即X，Y间存在直线关系。

调整后的直线方程Y=-23.73589044+0.26919749X

两组调整均值分别是Y(1)=125.993129, Y(2)=124.453741

因P=0.7829>0.05，故认为两调整均值间无显著差异。

第三节非参数检验(NPARlWAY过程)

前面介绍的t检验和方差分析属于参数的假设检验，其前提条件是：总体的分布为正态分布。但在许多实际问题中，总体的分布形式很难确定，或者总体的分布为偏态的，在这种情况下，参数的假设检验就不适用了。一般可选用非参数统计方法。

参数检验和非参数检验

单样本z和 t检验（proc means,proc

ttest）卡方检验、符号检验、Wilcoxon符号秩检验、K_S检验；

（proc freq, proc univariate）

成对样本

两独立样本Wilcoxon符号秩检验，

Mann Whitney U检验、

K_S检验，（proc

npar1way ）

多个独立样本检验（多重比较）方差分析（ proc anova, proc

glm ）

Kruskal-Willis Test、

Friedman Test（ proc

npar1way ）

相关性Pearson 相关（ proc corr）Spearman、Kendall秩

相关（ proc corr ）在SAS／STAT中，NPARlWAY过程专门用来处理单因素的非参数性检验。它基于经验分布函数和跨过单向分类的因变量的秩得分，计算出几个统计量，用以检验变量的分布在跨过不同组时有相同的位置参数。秩得分包括Wilcoxon得分、中位数得分、Savage得分和Van der Waerden得分。

1．过程格式

PROC NPARlWAY [选择项]；

CLASS 变量；

VAR 变量表；

OUTPUT OUT=数据集名统计关键字=名字；

BY 变量；

2．说明

(1)PROC NPARlWAY语句中的选择项

WILCOXON 要求处理数据的秩分析或Wilcoxon得分。对于两样本，进行Wilcoxon秩和检验(U检验)，对于多样本，给出Kruskal-Wallis 检验(H检验)。

ANOVA对原始数据进行方差分析（多样本）

NOPRINT抑制OUTPUT窗口输出；

(2)CLASS语句是必备语句，用来指定一个且仅一个分类变量。

(3)VAR语句指定用于分析的因变量。若缺省，则分析所有数值型变量。

3．举例

(1)两样本的差异显著性检验

某地测定急性克山病患者与克山病区健康人的血磷值如下表所示，问克山病患者与健康人血磷值有无显著差异?

程序：（yp129.sas）

编程说明：

数据步中，变量a代表因素，其水平数是2，变量n代表每组的样本数。过程步中，用CLASS语句指明单因素变量a，用VAR语句指明分析变量，并作Wilcoxon秩和检验。

结果说明：

这是两样本的检验，只需选择Wilcoxon检验的结果。因Z=2.23347, P=O.0255

(2)多样本问的检验

把主诉无月经患者，根据其基础代谢量分为甲状腺功能亢进、正常及低下三组，测定尿中总雌激素排泄量如下表所示，问三组间有无显著差异?用Kruskal—Wallis的H检验法(多个样本的比较)检验并作出结论。

表：尿中总雌激素排泄量

甲状腺功能尿中总雌激素排泄量

亢进13 6.8 7.8 15.5 11.4

正常14.1 7.9 6 6.7 5.2 10.7 3.8 9.4 13.9

低下9.2 11.6 4.2

程序： (yp130.sas)

编程说明：

数据步中，变量a代表因素，其水平数是3，变量n代表每组的样本数，过程步中，用CLASS语句指明单因素变量a，用VAR语句指明分析变量，并作Wilcoxon秩和检验。

结果说明：

因CHISQ=1.3455，P=0.5103>0.05故差异不显著，即甲状腺功能不同的无月经患者尿中雌激素的排泄量并无显著差异。

附录：多重比较的问题

问题：相对于两两比较，多重比较本质上是实施了多重的两两比较，在单个两两比较的显著水平α给定的情况下，m 重两两比较的第一类错误为()m α-1，（假定各重比较是两两独立的），因此，多重比较中一个真实的原假设被拒绝的可能性就越大，除非我们对这一类错误加以控制。

由概率论中Bonferroni 不等式，如果m 次比较时的m 个置信区间的置信水平分别为m A A A ,,,21 ，那么这m 个区间合在一起的置信水平大于等

于∑=--m

i i A 1)1(1，误差率下界为：m er )1(1αα--≥，这里α就是每次两两比较

的显著水平。

（比如，m=2，05.0=α，则两重联合检验的显著水平应等于

1-（1-0.05）*（1-0.05）=0.0975≈0.1，差不多是分开检验时显著水平0.05的两倍！）

由此，当任意给定k 组样本，检验他们中心是否相等的两两比较的所有可能性应该有m=2k C 种。当k 稍有增大时，误差率的增大是惊人的：

05.0=α，当k=3时，? m=3, 14.0=er α

05.0=α, 当k=4时, ? m=6, 26.0=er α

05.0=α, 当k=5时, ? m=10, 40.0=er α

误差率如此高的增长，对含多重比较的数据分析来说是灾难性的，必须加以控制。