生物统计学教案(7).doc

生物统计学教案第七章 拟合优度检验教学时间：2学时教学方法：课堂板书讲授教学目的：重点掌握二项分布的检验、正态性的检验，掌握独立性检验，了解X2的可加性。讲授难点：正态性的检验、二项分布的检验7.1 拟合优度检验的一般原理7.1.1 什么是拟合优度检验用来检验实际观测数与依照某种假设或模型计算出来的理论数之间的一致性的方法。可分为两种类型：（1）拟合优度检验：检验观测数与理论数之间的一致性。（2）独立性检验：通过检验实际观测数与理论数之间的一致性来判断事件之间的独立性。7.1.2 拟合优度检验的统计量例 黄圆豌豆与绿皱豌豆杂交，第二代分离数目如下： 黄圆 黄皱 绿圆 绿皱 总计实测数(Oi) 315(O1) 101（O2） 108(O3) 32(O4) 556理论数(Ti) 312.75(T1) 104.25(T2) 104.25(T3) 34.75(T4) 556拟合优度的一般做法是： （1）将观测值分为k种不同类别，如四种类型豌豆。（2）共获得n个独立观测值，第i类观测值的数目为Oi。如O1O4，他们的和等于n。（3）第i类的概率为pi，如上述四类豌豆的概率分别为9/16、3/16、3/16、1/16，概率之和等于1。（4）第i类的理论数Ti = npi, k个理论数之和等于n。如上例中的T1T4,它们的和等于n。（5）Oi与Ti不符合程度的计算： 求k个OiTi之和，显然它们恒等于0。 求k个(OiTi)2之和，得不出相对的不符合程度。Oi9、Ti6，OiTi3；Oi49、Ti46，OiTi3。前者的不符合程度远大于后者。 求k个(OiTi)/Ti2之和，但仍有问题。如：Oi8、Ti5以及Oi80、Ti50时OiTi/Ti都等于0.6。 为了解决上述问题，以Ti为权求加权值。由上式所定义的统计量也称为2。近似服从2分布，可由2分布表中查出临界值。 条件：当理论数小于5和df1时，上式与2分布偏离较大，因此：1) 当理论数小于5时，应与相邻项合并直到5。2)当df1时应做连续型矫正，矫正方法如下： 2的自由度为：dfk1a当理论数已经给定或计算理论数时所用的参数已知时a0。若总体参数没有给出，需由样本数据估计，这时a 0，a为需由样本估计的参数的个数。7.2 拟合优度检验7.2.1 一般程序 1、对数据进行分组2、根据总体参数计算理论数Ti。这时dfk13、由样本数据估计参数并理论数Ti。这时dfk1a。a为所估计参数的个数。4、合并理论数小于5的各组，记合并后的组数为k。5、零假设是观测数与理论数符合，拟合优度2检验为非参数统计，零假设可形象地记为：H0：OT0。6、计算出2值，与临界值比较，当22时拒绝H0。7.2.2 参数已知时，二项分布的检验例 检验上一节给出的例子。理论数均大于5，df1，已知，H0：OT = 0， 0.05。将数据代入公式。从附表中查出23, 0.057.815，220.05 结论：正常翅与残翅的分离比不符合31以上的计算是在df1但未作矫正时所得结果，下面计算矫正后的2。解 正常翅 残翅 |OT|0.5 16.5 16.5 (|OT|0.5)2 272.25 272.25 (|OT|0.5)2/T 0.926 2.778 2 = 0.926+2.778 = 3.704 H0：OT0，0.05，df1，20.053.841，220.05 结论：正常翅与残翅的分离比符合31从以上的结果可以看出，同一个问题矫正与不矫正所得结论不同。因为矫正后的结果比矫正前的低，若未矫正已经接受了H0，可以不再矫正；若未矫正时拒绝H0，则一定要矫正。7.3 独立性检验7.3.1 列联表2检验列联表2检验属独立性检验。例 下表给出不同给药方式与给药效果 给药方式 有效（A） 无效（） 总数 有效率 口服(B) 58 40 98 59.2 注射() 64 31 95 67.4 总 数 122 71 193 列联表2检验的原理现在要考虑的是给药方式与给药效果有无关联，如果有关联，即不同的给药方式产生不同的效果；反之，如果无关联，即不同的给药方式的治疗效果没有不同。从另一个角度讲，我们要考虑的是不同的给药方式与给药效果之间是否相互独立，因此列联表2检验又称为独立性检验。上表中没有理论数，需要由实测数推出。计算理论数的前提是假设给药方式与给药效果并无关联。在无关联的假设下，计算出理论数，然后比较实测数与理论数之间是否吻合。如果吻合，说明不同给药方式与给药效果是无关联的，不同的给药方式并不影响治疗效果。如果不吻合，说明不同的给药方式，其治疗效果不同。 列联表2检验的步骤1、提出零假设：假设实测数与给药方式和给药效果并无关联的前提下所计算出的理论数之间无差异。即H0：OT0。2、计算理论数：若事件A和事件B是相互独立的,则P(AB)=P(A)P(B)。在给药方式和效果之间是相互独立的前提下，计算口服（事件B）有效（事件A）的概率P(BA)P(B)P(A) = (98/193) (122/193)。其理论数T1(98/193)(122/193)(193) = (98)(122)/93 = 61.15。以类似的方法可以计算出另外三个理论数。3、计算2值：如同吻合度检验那样计算2值，若22，则拒绝H0。与吻合度检验一样，理论数不得小于5，当理论数小于5时应使用另外的方法。 4、确定自由度：因为每一行的各理论数受该行总数约束，每一列的各理论数受该列总数约束，所以自由度df(r-1)(c-1)。上例的计算结果如下： 有 效 无 效 总 数 口服 O1=58 O2=40 98 T1=(98)(122)/193=61.95 T2=(98)(71)/193=36.05 注射 O3=64 O4=31 95 T3=(95)(122)/193=60.05 T4=(95)(71)/193=34.95 总数 122 71 193结论：接受H0，不同给药方式的治疗效果没有显著不同。本例的df=1应当矫正