在一个实验中设k (k>2)个处理组,其中k-1个实验组,一个共同的对照组。实验结果经方差分析各组均数间有差别,进一步作两两比较时,若只需将实验组与对照组比较,而各实验组间不作相互比较时,可用Dunnett检验(C. W. Dunnett,1955)。它与 “多个样本均数间的两两比较”类似,但标准误的计算方法不同。
Dunnett检验的方法步骤如下:
(1) 列出各组均数及方差分析表,将各组 (包括对照组)均数按大小顺序排列,列出对比组及对比范围内的组数a。如将表3的均数排列为:
各对比组为C与甲、C与乙、C与丙,相应的a为2、3、4。
(2) 计算各实验组与对照组均数之差T-C。
(3) 计算T-C的标准误ST-C。若对比的实验组例数NT与对照组例数NC相等时用式(1),不等时用式(2)。
式中MS组内为方差分析表中的组内均方,式(1)中n为每组例数。
(4) 计算各对比组的统计量q′值,
(5) 以组内均方的自由度v和组数a查q′界值表。由于各对比组的a不全相同,故q′界值也不同,是为多重标准。亦有人一律用a=k,得一公用的q′界值,是为单一标准,但界值较大。双侧检验查表1、单侧检验查表2得P值,按所取检验水准作出推断结论。
C. W. Dunnett提出: 当实验组与对照组的方差不齐时,除sT-C的算法与上法不同外, 表1中的界值须作调整;当对照组例数比任何实验组例数均多时,表1中的界值亦须调整,详见表1引文。
表1 Dunnett检验用q′界值表(双侧)
上行:p = 0.05,下行: P=0.01
v | 组 数,a | 2 | 3 | 4 | 5 | 6 | ||||
5 | 2.57 | 3.03 | 3.29 | 3.48 | 3.62 | |||||
6 | 2.45 | 2.86 | 3.10 | 3.26 | 3.39 | |||||
7 | 2.36 | 2.75 | 2.97 | 3.12 | 3.24 | |||||
8 | 2.31 | 2.67 | 2.88 | 3.02 | 3.13 | |||||
9 | 2.26 | 2.61 | 2.81 | 2.95 | 3.05 | |||||
10 | 2.23 | 2.57 | 2.76 | 2.89 | 2.99 | |||||
11 | 2.20 | 2.53 | 2.72 | 2.84 | 2.94 | |||||
12 | 2.18 | 2.50 | 2.68 | 2.81 | 2.90 | |||||
13 | 2.16 | 2.48 | 2.65 | 2.78 | 2.87 | |||||
14 | 2.14 | 2.46 | 2.63 | 2.75 | 2.84 | |||||
15 | 2.13 | 2.44 | 2.61 | 2.73 | 2.82 | |||||
16 | 2.12 | 2.42 | 2.59 | 2.71 | 2.80 | |||||
17 | 2.11 | 2.41 | 2.58 | 2.69 | 2.78 | |||||
18 | 2.10 | 2.40 | 2.56 | 2.68 | 2.76 | |||||
19 | 2.09 | 2.39 | 2.55 | 2.66 | 2.75 | |||||
20 | 2.09 | 2.38 | 2.54 | 2.65 | 2.73 | |||||
24 | 2.06 | 2.35 | 2.51 | 2.61 | 2.70 | |||||
30 | 2.04 | 2.32 | 2.47 | 2.58 | 2.66 | |||||
40 | 2.02 | 2.29 | 2.44 | 2.54 | 2.62 | |||||
60 | 2.00 | 2.27 | 2.41 | 2.51 | 2.58 | |||||
120 | 1.98 | 2.24 | 2.38 | 2.47 | 2.55 | |||||
∞ | 1.96 | 2.21 | 2.35 | 2.44 | 2.51 |
摘自 Dunnett CW: New tables for multiple compa-risons with a control,Biometrics,20: 482,1964
表2 Dunnett检验用q′界值表(单侧)
上行:P=0.05,下行: P= 0.01
v | 组 数, a | 2 | 3 | 4 | 5 | 6 | ||||
5 | 2.02 | 2.44 | 2.68 | 2.85 | 2.98 | |||||
6 | 1.94 | 2.34 | 2.56 | 2.71 | 2.83 | |||||
7 | 1.89 | 2.27 | 2.48 | 2.62 | 2.73 | |||||
8 | 1.86 | 2.22 | 2.42 | 2.55 | 2.66 | |||||
9 | 1.83 | 2.18 | 2.37 | 2.50 | 2.60 | |||||
10 | 1.81 | 2.15 | 2.34 | 2.47 | 2.56 | |||||
11 | 1.80 | 2.13 | 2.31 | 2.44 | 2.53 | |||||
12 | 1.78 | 2.11 | 2.29 | 2.41 | 2.50 | |||||
13 | 1.77 | 2.09 | 2.27 | 2.39 | 2.48 | |||||
14 | 1.76 | 2.08 | 2.25 | 2.37 | 2.46 | |||||
15 | 1.75 | 2.07 | 2.24 | 2.36 | 2.44 | |||||
16 | 1.75 | 2.06 | 2.23 | 2.34 | 2.43 | |||||
17 | 1.74 | 2.05 | 2.22 | 2.33 | 2.42 | |||||
18 | 1.73 | 2.04 | 2.21 | 2.32 | 2.41 | |||||
19 | 1.73 | 2.03 | 2.20 | 2.31 | 2.40 | |||||
20 | 1.72 | 2.03 | 2.19 | 2.30 | 2.39 | |||||
24 | 1.71 | 2.01 | 2.17 | 2.28 | 2.36 | |||||
30 | 1.70 | 1.99 | 2.15 | 2.25 | 2.33 | |||||
40 | 1.68 | 1.97 | 2.13 | 2.23 | 2.31 | |||||
60 | 1.67 | 1.95 | 2.10 | 2.21 | 2.28 | |||||
120 | 1.66 | 1.93 | 2.08 | 2.18 | 2.26 | |||||
∞ | 1.64 | 1.92 | 2.06 | 2.16 | 2.23 |
摘自 Dunnett CW:A multiple comparison prece-dure for comparing several treatments witha control,J. Amer. Statist. Assoc.,50: 1096,1955
例 试验三菱莪术的抑癌作用,将小白鼠分四组,先致癌,然后对甲、乙、丙组依次注药0.5ml、1.0ml、1.5ml,对照组不用药。表3是各鼠的肿瘤重量。问各用药组的瘤重是否减轻?
表3 三个实验组与一个对照组小白鼠的肿瘤重量(g)
对照组 | 实 验 组 | 甲 | 乙 | 丙 | |||
X | 3.6 5.6 | 3.0 3.7 | 0.4 1.3 | 3.3 3.2 | |||
n | 10 | 10 | 10 | 10 | |||
| 4.66 | 2.50 | 2.46 | 1.87 |
先就表3资料作方差分析。
H0: μ0=μr1=μr2=μr3,
H1:四个总体均数不等或不全相等。
α=0.05。
结果见表4,按α=0.05水准拒绝H0,接受H1,可认为四个总体均数不等。将三个实验组均数分别与对照组均数作比较。
表4 表3资料的方差分析
变异来源 | SS | v | MS | F | P |
总 | 86.74 | 39 |
|
|
|
组 间 | 45.09 | 3 | 15.0300 | 12.99 | <0.01 |
H0: μr=μC,
H1:μT<μC。
单侧α=0.05。
四组均数的大小顺序见表3末行,由此确定各对比的实验组,两均数之差|R-c|及组数α,见表5第(1)、(2)、(3)栏。按式(1)得
按式(3)得各对比组q′值见表5第(4)栏。按v=36及第 (3)栏α,查表2得各P值,见第(5)栏。按α=0.05水准,三种比较均拒绝H0,接受H1,故可认为三个实验组的肿瘤平均重量都低于对照组,即三种剂量的三菱莪术对于小白鼠均有抑癌作用。
表5 三个实验组与对照组均数的比较
对比的实验组 | |-| | a | q' | P |
甲 组 | 2.16 | 2 | 4.49 | <0.01 |
如按单一标准(ν=36,α=4),上述三种比较仍得P<0.01,结果相同。