第十一章非参数检验 前面有关章节讨论的参数检验都要求总体服从一定的分布,对总体参数的检验是建立在 这种分布基础上的。例如,两样本平均数比较的t检验和多个样本平均数比较的F检验,都 要求总体服从正态分布,推断两个或多个总体平均数是否相等。本章引入另一类检验—一—非 参数检验(non- parametric test)。非参数检验是一种与总体分布状况无关的检验方法,它不 依赖于总体分布的形式,应用时可以不考虑被研究的对象为何种分布以及分布是否已知。非 参数检验主要是利用样本数据之间的大小比较及大小顺序,对两个或多个样本所属总体是否 相同进行检验,而不对总体分布的参数如平均数、标准差等进行统计推断。当样本观测值的 总体分布类型未知或知之甚少,无法肯定其性质,特别是观测值明显偏离正态分布,不具备 参数检验的应用条件时,常用非参数检验。非参数检验具有计算简便、直观,易于掌握,检 验速度较快等优点 非参数检验法从实质上讲,只是检验总体分布的位置(中位数)是否相同,所以对于总 体分布已知的样本也可以采用非参数检验法,但是由于它不能充分利用样本内所有的数量信 息,检验的效率一般要低于参数检验方法。例如,非配对资料的秩和检验,其效率为检验 的864%,就是说以相同概率判断出差异显著,t检验所需的样本个数要少13.6%。非参数 检验内容很多,本章只介绍常用的符号检验( sign test),秩和检验(rank- sum test)和等级 相关分析( rank correlation analysis)三种。 第一节符号检验 、配对资料的符号检验 )配对资料符号检验的意义配对资料符号检验是根据样本各对数据之差的 正负符号多少来检验两个总体分布位置的异同,而不去考虑差值的大小。每对数据之差为正 值用“+”表示,负值用“一”表示。可以设想如果两个总体分布位置相同,则正或负出现 的次数应该相等。若不完全相等,至少不应相差过大,否则超过一定的临界值就认为两个样 本所来自的两个总体差异显著,分布的位置不同。显然这种检验比较的是中位数而不是平均 数,当分布对称时,中位数与平均数相等。 (二)配对资料符号检验的基本步骤 1、提出无效假设与备择假设 HO:甲、乙两处理差值d总体中位数=0; H4:甲、乙两处理差值d总体中位数≠0。 此时进行两尾检验。若将HA中的“≠”改为“<”或“>”,则进行一尾检验。 2、计算差值并赋予符号求甲、乙两个处理的配对数据的差值d,d>0者记为“+”, d<0者记为“-”,d=0记为“0”。统计“+”、“-”、“0”的个数,分别记为n,n-,n0,令 n=n+n。检验的统计量为K,等于n、n中的较小者,即K=mn{m4,n}。 207
207 第十一章 非参数检验 前面有关章节讨论的参数检验都要求总体服从一定的分布,对总体参数的检验是建立在 这种分布基础上的。例如,两样本平均数比较的 t 检验和多个样本平均数比较的 F 检验,都 要求总体服从正态分布,推断两个或多个总体平均数是否相等。本章引入另一类检验——非 参数检验(non-parametric test)。非参数检验是一种与总体分布状况无关的检验方法,它不 依赖于总体分布的形式,应用时可以不考虑被研究的对象为何种分布以及分布是否已知。非 参数检验主要是利用样本数据之间的大小比较及大小顺序,对两个或多个样本所属总体是否 相同进行检验,而不对总体分布的参数如平均数、标准差等进行统计推断。当样本观测值的 总体分布类型未知或知之甚少,无法肯定其性质,特别是观测值明显偏离正态分布,不具备 参数检验的应用条件时,常用非参数检验。非参数检验具有计算简便、直观,易于掌握,检 验速度较快等优点。 非参数检验法从实质上讲,只是检验总体分布的位置(中位数)是否相同,所以对于总 体分布已知的样本也可以采用非参数检验法,但是由于它不能充分利用样本内所有的数量信 息,检验的效率一般要低于参数检验方法。例如,非配对资料的秩和检验,其效率为 t 检验 的 86.4%,就是说以相同概率判断出差异显著,t 检验所需的样本个数要少 13.6%。非参数 检验内容很多,本章只介绍常用的符号检验(sign test),秩和检验(rank-sum test)和等级 相关分析(rank correlation analysis)三种。 第一节 符号检验 一、配对资料的符号检验 (一)配对资料符号检验的意义 配对资料符号检验是根据样本各对数据之差的 正负符号多少来检验两个总体分布位置的异同,而不去考虑差值的大小。每对数据之差为正 值用“+”表示,负值用“-”表示。可以设想如果两个总体分布位置相同,则正或负出现 的次数应该相等。若不完全相等,至少不应相差过大,否则超过一定的临界值就认为两个样 本所来自的两个总体差异显著,分布的位置不同。显然这种检验比较的是中位数而不是平均 数,当分布对称时,中位数与平均数相等。 (二)配对资料符号检验的基本步骤 1、提出无效假设与备择假设 HO:甲、乙两处理差值 d 总体中位数=0; HA:甲、乙两处理差值 d 总体中位数≠0。 此时进行两尾检验。若将 HA 中的“≠”改为“<”或“>”,则进行一尾检验。 2、计算差值并赋予符号 求甲、乙两个处理的配对数据的差值 d,d>0 者记为“+”, d<0 者记为“-”,d=0 记为“0”。统计“+”、“-”、“0”的个数,分别记为 0 n+ ,n− ,n ,令 n = n+ + n− 。检验的统计量为 K,等于 n+ 、 n− 中的较小者,即 min{ , } K = n+ n−
3、统计推断由n查附表11符号检验用K临界值表(表中Pa表示两尾概率,用 于两尾检验,Pu表示一尾概率,用于一尾检验)得临界值K0sm,Kaom。如果K> Ko.05m P>0.05,则不能否定Ho,表明两个试验处理差异不显著:如果K0mm<K≤Kasm,0.01 <P≤005,则否定H,接受H4,表明两个试验处理差异显著;如果K≤K0m,P≤0.01, 则否定Ho,接受H4,表明两个试验处理差异极显著(注意:当K恰好等于临界K值时 其确切概率常小于附表11中列出的相应概率)。 【例11.1】某研究测定了噪声刺激前后15头猪的心率,结果见表11-1。问噪声对猪的 心率有无影响? 表11-1猪噪声刺激前后的心率(次/分钟 刺激前617068738581656272847660807971 刺激后757985778487887674818578888084 差值1491741623-14-239-1881-13 符号 这是一个配对资料两尾检验的问题 1、提出无效假设与备择假设 Ho:噪声刺激前后猪的心率差值d总体中位数=0; H1:噪声刺激前后猪的心率差值d总体中位数≠0 2、计算差值并赋予符号噪声刺激前后的差值及符号列于表11第4行和第5行, 从而得n4=2、n=13,n=n+n=2+13=15,K=mm{n+,n}=n4=2。 3、统计推断当m=15时,查附表11得临界值 Ko.05/15=3, Koalas=2,因为K=2=Kn0sy, P≤0.01,表明噪声刺激对猪的心率影响极显著。 值得注意的是,虽然符号检验方法简单,但是,由于利用的信息较少,所以效率较低, 且样本的配对数少于6时,不能检验出差别,在7-—12时也不敏感,配对数在20以上时符 号检验才较为有用。 二、样本中位数与总体中位数比较的符号检验 为了判断一个样本是否来自某已知中位数的总体,即样本所在总体的中位数是否等于某 已知总体的中位数,就需要进行样本中位数与总体中位数的差异显著性检验。其符号检验 的基本步骤为: 1、提出无效假设与备择假设 Ho:样本所在的总体中位数=已知总体中位数 H4:样本所在的总体中位数≠已知总体中位数 此时进行两尾检验。如果将备择假设H中的“≠”改为“<”或“>”,则进行一尾检 2、计算差值、确定符号及其个数将样本各观测值中大于已知总体中位数者记为 “+”,小于者记为“一”,等于者记为“0”。统计“+”、“一”、“0”的个数,分别记为 n、n、n,令n=n4+n。假设检验的统计量K为n,、n中的较小者,即K=min{m+,n} 3、统计推断由n查附表11符号检验用K临界值表,得临界值 Ko.osini,Knmm。如
208 3、统计推断 由 n 查附表 11 符号检验用 K 临界值表(表中 P(2)表示两尾概率,用 于两尾检验,P(1)表示一尾概率,用于一尾检验)得临界值 K0.05(n),K0.01(n)。如果 K>K0.05(n) , P>0.05,则不能否定 HO,表明两个试验处理差异不显著;如果 K0.01(n) <K≤K0.05(n) ,0.01 <P≤0.05,则否定 HO,接受 HA,表明两个试验处理差异显著;如果 K≤K0.01(n),P≤0.01, 则否定 HO,接受 HA,表明两个试验处理差异极显著(注意:当 K 恰好等于临界 K 值时, 其确切概率常小于附表 11 中列出的相应概率)。 【例 11.1】某研究测定了噪声刺激前后 15 头猪的心率,结果见表 11-1。问噪声对猪的 心率有无影响? 表 11-1 猪噪声刺激前后的心率(次/分钟) 猪 号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 刺激前 61 70 68 73 85 81 65 62 72 84 76 60 80 79 71 刺激后 75 79 85 77 84 87 88 76 74 81 85 78 88 80 84 差 值 -14 -9 -17 -4 1 -6 -23 -14 -2 3 -9 -18 -8 -1 -13 符 号 - - - - + - - - - + - - - - - 这是一个配对资料两尾检验的问题。 1、提出无效假设与备择假设 HO :噪声刺激前后猪的心率差值 d 总体中位数=0; HA :噪声刺激前后猪的心率差值 d 总体中位数≠0。 2、计算差值并赋予符号 噪声刺激前后的差值及符号列于表 11-1 第 4 行和第 5 行, 从而得 n+ = 2、 n− =13 , n = n+ + n− = 2 +13 =15 , K = min{n+ ,n− } = n+ = 2。 3、统计推断 当 n=15 时,查附表 11 得临界值 K0.05(15)=3,K0.01(15)=2,因为 K=2= K0.01(15), P≤0.01,表明噪声刺激对猪的心率影响极显著。 值得注意的是,虽然符号检验方法简单,但是,由于利用的信息较少,所以效率较低, 且样本的配对数少于 6 时,不能检验出差别,在 7—12 时也不敏感,配对数在 20 以上时符 号检验才较为有用。 二、样本中位数与总体中位数比较的符号检验 为了判断一个样本是否来自某已知中位数的总体,即样本所在总体的中位数是否等于某 一已知总体的中位数,就需要进行样本中位数与总体中位数的差异显著性检验。其符号检验 的基本步骤为: 1、提出无效假设与备择假设 HO:样本所在的总体中位数=已知总体中位数; HA:样本所在的总体中位数≠已知总体中位数。 此时进行两尾检验。如果将备择假设 HA 中的“≠”改为“<”或“>”,则进行一尾检 验。 2、计算差值、确定符号及其个数 将样本各观测值中大于已知总体中位数者记为 “+”,小于者记为“-”,等于者记为“0”。统计“+”、 “-”、 “0” 的个数,分别记为 n+ 、n−、 0 n ,令 n = n+ + n− 。假设检验的统计量 K 为 n+ 、n− 中的较小者,即 min{ , } K = n+ n− 。 3、统计推断 由 n 查附表 11 符号检验用 K 临界值表,得临界值 K0.05(n),K0.01(n)。如
果K>K0sm),P>0.05,则不能否定Ho,表明样本中位数与已知总体中位数差异不显著 如果Km)<K≤Ksm,001<P≤005,则否定Ho,接受H,表明样本中位数与已知总 体中位数差异显著:如果K≤K0m,P≤0.01,则否定Ho,接受HA,表明样本中位数与己 知总体中位数差异极显著。 【例11.2】已知某品种成年公黄牛胸围平均数为140厘米,今在某地随机抽取10头该 品种成年公黄牛,测得一组胸围数字:128.1,144.4,150.3,146.2,140.6,139.7,134.1, 124.3,147.9,143.0(cm)。问该地成年公黄牛胸围与该品种胸围平均数是否有显著差异? 表11-2成年公黄牛胸围测定值符号检验表 牛号1 6 7 10 胸围12811441503146214061397134112431479143 差值119446.3620.6 0.35.9-15.779 符号 1、提出无效假设与备择假设 Ho:该地成年公黄牛胸围的平均数=140厘米 H1:该地成年公黄牛胸围的平均数≠140厘米 2、计算差值、确定符号及其个数样本各观测值与总体平均数的差值及其符号列 于表11-3,并由此得n4=6,n=4,n=6+4=10,K=mm{n,n}=n=4。 3、统计推断由m=10,查附表1,得K0mo=1,K>kasn0,P>0.05,不能否定 Ho,表明样本平均数与总体平均数差异不显著,可以认为该地成年公黄牛胸围的平均数与 该品种胸围总体平均数相同。 第二节秩和检验 秩和检验也叫做符号秩和检验( signed rank- sum test),是一种经过改进的符号检验 或称 Wilcoxon检验,其统计效率远较符号检验为高。因为它除了比较各对数据差值的符号 外,还要比较各对数据差值大小的秩次高低。方法是通过将观测值按由小到大的次序排列, 编定秩次,求出秩和进行假设检验。秩和检验与符号检验法不同,要求差数来自某些对称分 布的总体,但并不要求每一差数来自相同的分布 、配对试验资料的符号秩和检验( Wilcoxon配对法) (一)基本步骤 1、提出无效假设与备择假设 :差值d总体的中位数=0 HA:差值d总体的中位数≠0 此时进行两尾检验。若将H1中的“≠”改为“<”或“>”,则进行一尾检验。 2、编秩次、定符号先求配对数据的差值d,然后按d的绝对值从小到大编秩次 再根据原差值正负在各秩次前标上正负号,若差值d=0,则舍去不记,若有若干个差值d的 绝对值相等,则取其平均秩次 209
209 果 K>K0.05(n) ,P>0.05,则不能否定 HO,表明样本中位数与已知总体中位数差异不显著; 如果 K0.01(n) <K≤K0.05(n) ,0.01<P≤0.05,则否定 HO,接受 HA,表明样本中位数与已知总 体中位数差异显著;如果 K≤K0.01(n),P≤0.01,则否定 HO,接受 HA,表明样本中位数与已 知总体中位数差异极显著。 【例 11.2】已知某品种成年公黄牛胸围平均数为 140 厘米,今在某地随机抽取 10 头该 品种成年公黄牛,测得一组胸围数字:128.1, 144.4, 150.3, 146.2, 140.6, 139.7, 134.1, 124.3, 147.9, 143.0(cm)。 问该地成年公黄牛胸围与该品种胸围平均数是否有显著差异? 表 11-2 成年公黄牛胸围测定值符号检验表 牛号 1 2 3 4 5 6 7 8 9 10 胸围 128.1 144.4 150.3 146.2 140.6 139.7 134.1 124.3 147.9 143 差值 -11.9 4.4 6.3 6.2 0.6 -0.3 -5.9 -15.7 7.9 3 符号 - + + + + - - - + + 1、提出无效假设与备择假设 HO :该地成年公黄牛胸围的平均数=140 厘米, HA :该地成年公黄牛胸围的平均数≠140 厘米。 2、计算差值、确定符号及其个数 样本各观测值与总体平均数的差值及其符号列 于表 11-3,并由此得 = 6, = 4, n+ n− n = 6 + 4 = 10, K = min{n+ ,n− } = n− = 4 。 3、统计推断 由 n=10,查附表 11,得 K0.05(10)=1,K>K0.05(10) ,P>0.05,不能否定 HO,表明样本平均数与总体平均数差异不显著,可以认为该地成年公黄牛胸围的平均数与 该品种胸围总体平均数相同。 第二节 秩和检验 秩和检验也叫做符号秩和检验(signed rank-sum test),是一种经过改进的符号检验, 或称 Wilcoxon 检验,其统计效率远较符号检验为高。因为它除了比较各对数据差值的符号 外,还要比较各对数据差值大小的秩次高低。方法是通过将观测值按由小到大的次序排列, 编定秩次,求出秩和进行假设检验。秩和检验与符号检验法不同,要求差数来自某些对称分 布的总体,但并不要求每一差数来自相同的分布。 一、配对试验资料的符号秩和检验(Wilcoxon 配对法) (一)基本步骤 1、提出无效假设与备择假设 HO:差值 d 总体的中位数=0; HA:差值 d 总体的中位数≠0。 此时进行两尾检验。若将 HA 中的“≠”改为“<”或“>”,则进行一尾检验。 2、编秩次、定符号 先求配对数据的差值 d,然后按 d 的绝对值从小到大编秩次。 再根据原差值正负在各秩次前标上正负号,若差值 d=0,则舍去不记,若有若干个差值 d 的 绝对值相等,则取其平均秩次
3、确定统计量T分别计算正秩次及负秩次的和,并以绝对值较小的秩和绝对值为 检验的统计量T。 4、统计推断记正、负差值的总个数为n,根据n查附表101)符号秩和检验用T 临界值表,得T0msm,Taom。如果T>sm,P>O05,则不能否定Ho,表明两个试验处 理差异不显著:如果Tmm<T≤7sm,0.01<P≤0.05,则否定Ho,接受H,表明两个 试验处理差异显著:如果T≤0m,P≤001,则否定Ho,接受HA,表明两个试验处理差 异极显著(注意:当T恰好等于临界T值时,其确切概率常小于附表10()中列出的相应概 率)。 【例11.3】某试验用大白鼠研究饲料维生素E缺乏与肝脏中维生素A含量的关系,先 将大白鼠按性别、月龄、体重等配为10对,再把每对中的两只大白鼠随机分配到正常饲料 组和维生素E缺乏饲料组,试验结束后测定大白鼠肝中维生素A的含量如表11-4。试检验 两组大白鼠肝中维生素A的含量是否有显著差异。 表11-3不同饲料鼠肝维生素A含量资料(国际单位/克) 鼠对别 E常饲料组3550200031003000395038003620375034503050 维生素E缺乏组2450240031001800320032503620270027001750 差值d 1200750550 10507501300 秩次 7 3.5 1、提出无效假设与备择假设 HO:差值d总体的中位数=0 H4:差值d总体的中位数≠0。 2、编秩次、定符号计算表113中配对数据差值d,将d=0的舍去,共有差值n=8 个。按绝对值从小到大排列秩次并标上相应的符号,差值绝对值为750的有两个,它们的秩 次为3和4,所以其平均秩次为(3+4)/2=3.5,结果见表11-3。 3、确定统计量T此例,正号有7个,其秩次为2,3.5,35,5,6,7,8,秩次 和为:2+3.5+3.5+5+6+7=35 负号只有1个,其秩次为1,秩次和等于1 负号秩次和较小,所以7=1 4、统计推断由=8查附表10(1)得,TD08=3,Taom=0,因为8<7<700(p, 0.01<P<0.05,否定Ho,接受H4,表明两个试验处理差异显著 二、非配对试验资料的秩和检验( Wilcoxon非配对法) 非配对试验资料的秩和检验是关于分别抽自两个总体的两个独立样本之间秩和的成组 比较,它比配对资料的秩和检验的应用更为普遍。 (一)基本步骤 1、提出无效假设与备择假设 HO:甲样本所在的总体的中位数=乙样本所在的总体的中位数 HA:甲样本所在的总体的中位数≠乙样本所在的总体的中位数
210 3、确定统计量 T 分别计算正秩次及负秩次的和,并以绝对值较小的秩和绝对值为 检验的统计量 T。 4、统计推断 记正、负差值的总个数为 n,根据 n 查附表 10(1)符号秩和检验用 T 临界值表,得 T0.05(n),T0.01(n)。如果 T>T0.05(n) ,P>0.05,则不能否定 HO,表明两个试验处 理差异不显著;如果 T0.01(n) <T≤T0.05(n) ,0.01<P≤0.05,则否定 HO,接受 HA,表明两个 试验处理差异显著;如果 T≤T0.01(n),P≤0.01,则否定 HO,接受 HA,表明两个试验处理差 异极显著(注意:当 T 恰好等于临界 T 值时,其确切概率常小于附表 10(1)中列出的相应概 率)。 【例 11.3】某试验用大白鼠研究饲料维生素 E 缺乏与肝脏中维生素 A 含量的关系,先 将大白鼠按性别、月龄、体重等配为 10 对,再把每对中的两只大白鼠随机分配到正常饲料 组和维生素 E 缺乏饲料组,试验结束后测定大白鼠肝中维生素 A 的含量如表 11-4。试检验 两组大白鼠肝中维生素 A 的含量是否有显著差异。 表 11-3 不同饲料鼠肝维生素 A 含量资料(国际单位/克) 鼠对别 1 2 3 4 5 6 7 8 9 10 正常饲料组 3550 2000 3100 3000 3950 3800 3620 3750 3450 3050 维生素E缺乏组 2450 2400 3100 1800 3200 3250 3620 2700 2700 1750 差值di 1100 -400 0 1200 750 550 0 1050 750 1300 秩次 +6 -1 +7 +3.5 +2 +5 +3.5 +8 1、提出无效假设与备择假设 HO:差值 d 总体的中位数=0; HA:差值 d 总体的中位数≠0。 2、编秩次、定符号 计算表 11-3 中配对数据差值 di,将 d =0 的舍去,共有差值 n=8 个。按绝对值从小到大排列秩次并标上相应的符号,差值绝对值为 750 的有两个,它们的秩 次为 3 和 4,所以其平均秩次为(3+4)/2=3.5,结果见表 11-3。 3、确定统计量 T 此例,正号有 7 个,其秩次为 2,3.5,3.5,5,6,7,8,秩次 和为:2+3.5+3.5+5+6+7=35; 负号只有 1 个,其秩次为 1,秩次和等于 1。 负号秩次和较小,所以 T=1。 4、统计推断 由 n=8 查附表 10(1)得,T0.05(8)=3,T0.01(n)=0,因为 T0.01(8) <T<T0.05(8) , 0.01<P<0.05,否定 HO,接受 HA,表明两个试验处理差异显著。 二、非配对试验资料的秩和检验(Wilcoxon 非配对法) 非配对试验资料的秩和检验是关于分别抽自两个总体的两个独立样本之间秩和的成组 比较,它比配对资料的秩和检验的应用更为普遍。 (一)基本步骤 1、提出无效假设与备择假设 HO:甲样本所在的总体的中位数=乙样本所在的总体的中位数; HA:甲样本所在的总体的中位数≠乙样本所在的总体的中位数
此时进行两尾检验。若将H4中的“≠”改为“<”或“>”,则进行一尾检验 2、求两个样本合并数据的秩次假设两个样本的含量分别为m和n2,则将两样 本的观测值合并后,总的数据为n+m2个。将合并后的数据按从小到大的顺序排列,与每个 数据对应的序号即为该数据的秩次,最小数值的秩次为“1”,最大数值的秩次为“m+n2”。 遇不同样本的相同观测值时,其秩次取原秩次的平均秩次,但是同一样本内遇相同的观测值 时则不必求平均秩次,秩次孰先孰后都可以。 3、确定统计量T将两个样本重新分开,并计算各自的秩和。将较小的那个样本 含量作为m,其秩和作为检验的统计量T。若n=m2,则任取一组的秩和为T 4、统计推断由n、(n-m)查附表10(3),得接受区域705-105,Tio-7o。若 T在T05-105之内,P>005,则不能否定Ho,表明两个试验处理差异不显著;若T在 705-705之外但在T1-70之内,0.01<P≤0.05,则否定Ho,接受H4,表明两个试验 处理差异显著;若T在T0-70之外,P<0.01,则否定Ho,接受H,表明两个试验处理 差异极显著。 【例11.4】研究两种不同能量水平饲料对5-6周龄肉仔鸡增重(克)的影响,资料 如表11-4所示。问两种不同能量水平的饲料对肉仔鸡增重的影响有无差异? 表11-4两种不同能量水平饲料的肉仔鸡增重及秩和检验 肉仔鸡增重(g) 高能量603585598620617650 秩次128.51114 T1=73.5 低能量489457512567512585591531467n2=9 秩次3 8.5106 2T2=46.5 1、提出无效假设与备择假设 Ho:髙能量饲料增重总体的中位数=低能量饲料增重总体的中位数 H:高能量饲料增重总体的中位数≠低能量饲料增重总体的中位数 2、编秩次将两组数据混合从小到大排列为秩次。在低能量组有两个“512”,不求 平均秩次,其秩次分别为4和5;在高、低两组有一对数据为“585”,需求它们的平均秩次: (8+9)2=85。结果见表11-4 3、确定统计量T以较小样本的秩次和为统计量T,即T=73.5 4统计推断由n=6,nm=96=3查附表103)得,o05-To5为31-65,7o1-7o 为26-70。7=73.5在01-To,即26-70之外,P<0.01,否定hb,接受H,表明饲料 能量高低对肉仔鸡増重的影响差异极显著 、多个样本比较的秩和检验( Kruskal-wallis法,H法) 多个样本比较的秩和检验的 Kruskal- Wallis法,又称H检验法。该法的前提是假设抽 样总体是连续的和相同的,利用多个样本的秩和来推断它们分别代表的总体之分布位置是否 相同,检验的基本步骤是: 1、提出无效假设与备择假设 211
211 此时进行两尾检验。若将 HA 中的“≠”改为“<”或“>”,则进行一尾检验。 2、求两个样本合并数据的秩次 假设两个样本的含量分别为 n1 和 n2,则将两样 本的观测值合并后,总的数据为 n1+n2 个。将合并后的数据按从小到大的顺序排列,与每个 数据对应的序号即为该数据的秩次,最小数值的秩次为“1”,最大数值的秩次为“n1+n2”。 遇不同样本的相同观测值时,其秩次取原秩次的平均秩次,但是同一样本内遇相同的观测值 时则不必求平均秩次,秩次孰先孰后都可以。 3、确定统计量 T 将两个样本重新分开,并计算各自的秩和。将较小的那个样本 含量作为 n1,其秩和作为检验的统计量 T。若 n1=n2,则任取一组的秩和为 T。 4、统计推断 由 n1、(n2–n1)查附表 10(3),得接受区域 0.05 ' T0.05 − T , 0.01 ' T0.01 −T 。若 T 在 0.05 ' T0.05 − T 之内,P>0.05,则不能否定 HO,表明两个试验处理差异不显著;若 T 在 0.05 ' T0.05 −T 之外但在 0.01 ' T0.01 −T 之内,0.01<P≤0.05,则否定 HO,接受 HA,表明两个试验 处理差异显著;若 T 在 0.01 ' T0.01 −T 之外,P<0.01,则否定 HO,接受 HA,表明两个试验处理 差异极显著。 【例 11.4】 研究两种不同能量水平饲料对 5-6 周龄肉仔鸡增重(克)的影响,资料 如表 11-4 所示。问两种不同能量水平的饲料对肉仔鸡增重的影响有无差异? 表 11-4 两种不同能量水平饲料的肉仔鸡增重及秩和检验 饲 料 肉仔鸡增重(g) 高能量 秩 次 603 12 585 8.5 598 11 620 14 617 13 650 15 n1=6 T1=73.5 低能量 秩 次 489 3 457 1 512 4 567 7 512 5 585 8.5 591 10 531 6 467 2 n2=9 T2=46.5 1、提出无效假设与备择假设 HO:高能量饲料增重总体的中位数=低能量饲料增重总体的中位数; HA:高能量饲料增重总体的中位数≠低能量饲料增重总体的中位数。 2、编秩次 将两组数据混合从小到大排列为秩次。在低能量组有两个“512”,不求 平均秩次,其秩次分别为 4 和 5;在高、低两组有一对数据为“585”,需求它们的平均秩次: (8+9)/2=8.5。结果见表 11-4。 3、确定统计量 T 以较小样本的秩次和为统计量 T,即 T= 73.5。 4、统计推断 由n1=6, n2-n1=9-6=3查附表10(3)得, 0.05 ' T0.05 −T 为31—65, 0.01 ' T0.01 −T 为 26—70。T=73.5 在 0.01 ' T0.01 −T ,即 26—70 之外,P<0.01, 否定 HO,接受 HA,表明饲料 能量高低对肉仔鸡增重的影响差异极显著。 三、多个样本比较的秩和检验(Kruskal-Wallis 法,H 法) 多个样本比较的秩和检验的 Kruskal-Wallis 法,又称 H 检验法。该法的前提是假设抽 样总体是连续的和相同的,利用多个样本的秩和来推断它们分别代表的总体之分布位置是否 相同,检验的基本步骤是: 1、提出无效假设与备择假设