第九章体育统计学.docx

1、第九章体育统计学第九章 检 验检验是一种用途较广的显著性检验方法，它既适用于参数估计，也适用于非参数估计，同方差分析F检验一样，检验是一种能同时对多种资料进行检验的方法。在本章将阐述如何将统计假设检验的概念加以推广，使之适用于样本来自元素可分成三个或更多个类目的总体的情形。这里只介绍用检验在检验两个以上样本率或构成比之间差异时有无显著性的方法。第一节 分布与表在第七章叙述总体方差的假设检验时，曾简单提到分布。定义：设从标准正态分布N（0，1）中，取出n 个样品u i ( i =1,2, )，它们的平方和记作，即= u 12 + u 22 + + u n2 = （91）一般称为服从参数为 n 的

2、的分布。 （卡方分布）如果从一般正态分布N（，）中，抽出n 个样品x i , 因为u i = ，所以公式（91）转化为： = = （92）有时在实用中，用 代替 ，则有 = = 这同第七章中提过的相同。如果将 的各种不同数值连同相应的相对出现频数列出，就得到 = 的抽样分布，这就是自由度为n 的分布。图91 与几种自由度相对应的分布由图91可看到的样本分布形式有两个特点：（1）正态，而是正偏态分布，右侧无限延长以横轴为渐近线；（2）不同自由度的取样分布不同。随自由度增大，曲线趋于“对称”、“正态”分布。实际上，在自由度为30时，曲线基本对称且近似正态分布。分布的实际重要性是基本于这样的事实，即

3、对大样本而言，量X 2 = （93）的分布接近于的分布。式中：Ai = 观察频数，即“实际数”；Ti = 预期频数，即“理论数”；K = 给定情形下实际数与理论数的配对数。所以近似地有= （94）公式（94）是下面进行检验的常用公式。检验主要是根据样本算出的值与查值表所得的临界值作比较。若，P（H0），拒绝H0 ，则称在显著性水平上差异有显著性意义；否则称在水平上差异无显著性意义。分布的尾端面积或概率P值在附表“值表”中按不同的自由度n分别列出。至于自由度n是直接由资料划分类别的多少所决定的，与样本的大小无关。确定n的原则是：互相排斥的类别数k ，减去所受的限制数m ，即n= km。第二节 两

4、个样本率差别的检验当两个样率不同时，可能由两个原因造成：一是差异由抽样误差所致；二是它们来自不同质的总体，确有不同。为区分二者，则要做检验。检验的步骤通过下例说明。例91 （资料同712）为研究游泳与患慢性鼻炎有无关系，随机抽测游泳与田径两专业的学生进行对比，抽测资料如下：组别患病人数未患病人数合计患病率游泳专业20(13)60(67)8025%田径专业6(13)74(67)807. 5%合计2613416016. 25%问患慢性鼻炎是否与专业有关？（= 0. 05） 上表中，20 60这四个格子内的数字是整个表的基本数 6 74字，故这种资料称为“四格表资料”，或称为“22列联表”。解：（1

5、）无效假设 H 0 ： = （两专业患病率相同，均为合计的患病率16. 25%）（2）计算值：= 实际数A已知，理论数T（表中括号内的数字）是根据H 0 推算出来的：先由合计的160人患病26人知理论上的患病率为 =16. 25%，则游泳专业80人理论上的患病人数为T11 = 80 = 13由此可知理论数的计算公式为：T r c = （95）式中：r 行号, c 列号, n r r 行的合计数, n c c 列的合计数, n 总例数。这样，T12 = = 67T21 = = 13T22 = = 67= = = 3. 7692 + 0. 7313 + 3. 7692 + 0. 7313 = 9.

6、 001（3）取= 0. 05；“行列表”资料的自由度为： =（行数1）（列数1）=（r1）（c1）（96）本例= 21）（21）= 1查值表得到= 3. 84（4）判断结果：= 9. 001= 3. 84P（H0） 0. 05，拒绝H0 ，由数据上看，认为患慢性鼻炎与游泳专业有关。此结论与例711用U检验去作的结论一致。我们还可以用四格表专用公式计算值：= （97）式中a、b、c、d 分别为四格表资料实际数，总例数n = a + b + c + d。 把上例数据 20( a ) 60( b ) 代入公式（97）得： 6( c ) 74( d ) = = 9. 001结果与前法相同。另外需要注

7、意： 只是一种近似。在1T5，且n 40时，对于自由度为1 的四格表来说为减小偏差，需要计算校正的值。此时：= （98）= （99）公式（98）为公式（94）的校正；公式（99）为公式（97）的校正。第三节 多个样本率（行列表）的检验四格表是行列表中最简单的形式。至于行列表（亦称r c 表）的理论与计算步骤均与四格表类同。用公式= 去计算行列表的值较为麻烦，为简便可用下面公式，由实际数直接计算出值。= （910）例92 上例的结论认为患慢性鼻炎与游泳专业有关。为进一步了解患慢性鼻炎的发病率与游泳专项运动年限有无关系，收集调查资料如下表中。解：（1）无效假设H 0 ：发病率与运动年限无关。（2）

8、计算值：= = 281 ( + + + + + + + 1 ) 10 . 36（3）取= 0. 05，=（r1）（c1）=（41）（21）=3。查值表（书后附表6）知= 7. 81（4）判断结果：，拒绝H 0 ，认为不同运动年限的游泳运动员鼻炎发病率的差别的显著性意义。行列表检验的注意事项：（1）在行列表中如果有五分之一以上的格子的理论数小于5，或有一个格子的理论数小于1时，则应使理论数小于5的格子与邻组合并以增大理论数。否则易导致错误的结论。但合并时要注意是否合理，不同质的资料不可并组！此时只有增加观察例数再作统计分析。（2）做行列表检验得到P（H 0 ）0. 05，则拒绝无效假设H 0 ，

9、说明被比较的几个样本率之间总的差别有显著性，但不能作出任何两组间差别都有显著性的结论！要比较某两个样本率之间的差别如何，需另作检验。分布检验程序一、程序功能检验常用于符合度分析。应用本程序时如按程序所问输入行元素的个数、列元素的个数和每个元素值，程序就能自动计算出其值、自由度数及其结尾值（检验水平）。使用者将结尾值与给定的显著性水平比较即能得出检验的结果。二、程序中使用的主要符号说明Z：TOTAL XSQUARE VALUE 总计值；V：DEGREE OF FREEDOM 自由度数；X2：TAILEND VALUE 结尾值。三、程序所依据的理论计算公式Z = 其中：foi 为每一组的观察频数；

10、fei 为相应的理论频数。四、程序名称：XSQUARE . TES5REM XSquare. TES10PRINT “xsquare. Tes”20PRINT30DIM V1(25), V2(5), A(5), M(4)40INPUT “number of rows =”, R50PRINT “number of rows =”; R60INPUT “number of columns =”, C70PRINT “number of columns =”; C80PRINT “Contingency table : ”90FOR I = 1 TO R100PRINT “ROW” ; I110F

11、OR J = 1 TO C120PRINT “Element”; J;125PRINT “Element”, J;130INPUT E131EA = (I1)*C + J132V1(EA) = E135PRINT E136PRINT140NEXT J145PRINT150NEXT I160PRINT170L = 0180M = 1190FOR I = 1 TO R200FOR J = 1 TO C210A(I) = A(I) + V1(M)220M = M + 1230NEXT J240L = L + A(I)250NEXT I260N = R * C270FOR I = 1 TO C280F

12、OR J = 1 TO N STEP C290V2(I) = V2(I) + V1(J)300NEXT J310NEXT I320Z = 0330PRINT “Observed value” ; “ ” ; “Expected value” ; “ ” ; “XSQUARE CONTRIBUTION”340FOR I = 1 TO C350PRINT “Column” ; I360FOR J = 1 TO R370P = A(J) * V2(I) / L380X = I + (J1) * C390Y = (V1(X)P) 2 / P400Z = Z + Y410PTINT “ ”V1(X),

13、“ ”P, “ ” ; Y420NEXT J430 NEXT I440V = (C1) * (R1)450PRINT “Xsquare =” ; z460PRINT “Degree of freedom =” ; v470PRINT “Tial end value =” ;480R = 1490FOT I = V TO 2 STEP 2500R = R * I510NEXT I520K = Z (INT (V + 1) / 2) * EXP (Z / 2) / R530IF INT (V / 2) = V / 2 THEN 560540J = SQR (2 / Z / 3. 141592653

14、6)550GOTO 570560J = 1570L = 1580M = 1590V = V + 2600M = M * Z / V610IF M . 0000001 THEN 640620L = L + M630GOTO 590640X2 = 1J * K * L650PRINT X2660END五、例 题用A和B两种方法，对同一物体进行处理的试验结果如下表：类 别 有 效 无 效 A 44 10 B 81 35试有检验法检验这两种方法的处理结果有无显著性差别？（= 0. 05）解：按程序所问输入数据，计算得到结果。RUNxsquare. Tesnumber of rows=2number of columns=2Contingency table:Row 1Element 1 44Element 2 10ROW 2Element 1 81Element 2 35Observed value Expected value xSquare ContributionColumn 144 3