《生物统计学教案(七)》由会员分享，可在线阅读，更多相关《生物统计学教案(七)（6页珍藏版）》请在金锄头文库上搜索。
1、生物统计学教案（7）生物统计学教案第七章拟合优度检验教学时间：2学时教学方法：课堂板书讲授教学目的：重点掌握二项分布的检验、正态性的检验，掌握独立性检验，了解X2的可加性。讲授难点：正态性的检验、二项分布的检验7.1 拟合优度检验的一般原理7.1.1 什么是拟合优度检验用来检验实际观测数与依照某种假设或模型计算出来的理论数之间的一致性的方法。可分为两种类型：(1)拟合优度检验：检验观测数与理论数之间的一致性。(2)独立性检验：通过检验实际观测数与理论数之间的一致性来判断事件之间的独立性。7.1.2拟合优度检验的统计量例黄圆豌豆与绿皱豌豆杂交，第二代分离数目如下：黄圆黄皱绿圆绿皱总计实测数(Oi
2、)315(01)101(02)108(03)32(04)556理论数(Ti)312.75(T1)104.25(T2)104.25(T3)34.75(T4)556拟合优度的一般做法是：(1)将观测值分为k种不同类别，如四种类型豌豆。(2)共获得n个独立观测值，第i类观测值的数目为Oi。如Oi。4,他们的和等于n。(3)第i类的概率为pi,如上述四类豌豆的概率分别为9/16、3/16、3/16、1/16,概率之和等于1。(4)第i类的理论数Ti=np,k个理论数之和等于n。如上例中的一T4,它们的和等于n(5) Oi与Ti不符合程度的计算：求k个OiTi之和，显然它们包等于00 求k个(OiTi)
3、2之和，得不出相对的不符合程度。Oi=9、Ti=6,OiTi=3;Oi=49、Ti=46,OiTi=3。前者的不符合程度远大于后者。 求k个(OiTi)/Ti2之和，但仍有问题。如：Oi=8、Ti=5以及Oi=80、Ti=50时OiTi/Ti都等于0.6。 为了解决上述问题，以Ti为权求加权值。k、T ii =1Oi-Ti、/(Oi-T)2=乙Ti)i=iTi由上式所定义的统计量也称为Xo近似服从，分布，可由，分布表中查出临界值。 条件：当理论数小于5和df=1时，上式与/分布偏离较大，因此：1)当理论数小于5时，应与相邻项合并直到5。2)当df=1时应做连续型矫正，矫正方法,一一、21(|O
4、i-TiI-0.5)=L1 =1Ti如下：J的自由度为：df=k-1a当理论数已经给定或计算理论数时所用的参数已知时a=0o若总体参数没有给出，需由样本数据估计，这时aw0,a为需由样本估计的参数的个数。7.2拟合优度检验7.2.1 一般程序1、对数据进行分组2、根据总体参数计算理论数Ti。这时df=k-13、由样本数据估计参数并理论数Ti。这时df=k-1-a0a为所估计参数的个数。4、合并理论数小于5的各组，记合并后的组数为ko5、零假设是观测数与理论数符合，拟合优度，检验为非参数统计，零假设可形象地记为：H0:OT=0。6、计算出，值，与临界值比较，当/双2时拒绝H0O7.2.2 参数日
5、已知时，二项分布的检验例检验上一节给出的例子。理论数均大于5,df1,小已知，Ho：OT=0,a=0.05。将数据代入公式。2 JOi-Ti2f=Z1 =1Ti4_2_2222 315-312.75101-104.25108104.25 x = 0.926+2.778 = 3.70432-34.75312.75104.25104.2534.75=0.0160.1010.1350.218=0.470从附表中查出X3,0.05=7.815,衣：0.05。结论是接受H0,杂交结果符合9:3:3:1的分离比。例用正常翅的野生型果蝇与残翅果蝇杂交，F1代均表现为正常翅。F1代自交，在F2代中包含311个
6、正常翅，和81个残翅。问这一分离比是否符合孟德尔3：1的理论比？解正常翅残翅实际观测数31181392理论数29498392O-T（未矫正）1717(O-T)2289289(O-T)2/T0.9832.9492一一X=0.983+2.949=3.932H。:OT=0,a=0.05,df=1,(0.05=3.841,J,0.05结论：正常翅与残翅的分离比不符合3：1以上的计算是在df=1但未作矫正时所得结果，下面计算矫正后的x解正常翅残翅|OT|0.516.516.5(|O-T|-0.5)2272.25272.25(|O-T|-0.5)2/T0.9262.778Ho：OT=0,a=0.05,df
7、=1,病.05=3.841,羯名.05结论：正常翅与残翅的分离比符合3：1从以上的结果可以看出，同一个问题矫正与不矫正所得结论不同。因为矫正后的结果比矫正前的低，若未矫正已经接受了H。，可以不再矫正；若未矫正时拒绝H0,则一定要矫正。7.3独立性检验7.3.1 列联表演检验列联表2检验属独立性检验。例下表给出不同给药方式与给药效果后效(A)A后效率口服(B)58409859.2%注射(B)64319567.4%总数12271193列联表X2检验的原理现在要考虑的是给药方式与给药效果有无关联，如果有关联，即不同的给药方式产生不同的效果；反之，如果无关联，即不同的给药方式的治疗效果没有不同。从另一
8、个角度讲，我们要考虑的是不同的给药方式与给药效果之间是否相互独立，因此列联表J检验又称为独立性检验。上表中没有理论数，需要由实测数推出。计算理论数的前提是假设给药方式与给药效果并无关联。在无关联的假设下，计算出理论数，然后比较实测数与理论数之间是否吻合。如果吻合，说明不同给药方式与给药效果是无关联的，不同的给药方式并不影响治疗效果。如果不吻合，说明不同的给药方式，其治疗效果不同。列联表2检验的步骤1、提出零假设：假设实测数与给药方式和给药效果并无关联的前提下所计算出的理论数之间无差异。即H0：OT=0。2、计算理论数：若事件A和事件B是相互独立的，则P(AB)=P(A)P(B)。在给药方式和效
9、果之间是相互独立的前提下，计算口服(事件B)有效(事件A)的概率P(BA)=P(B)P(A)=(98/193)(122/193)。其理论数=(98/193)(122/193)(193)=(98)(122)/93=61.15。以类似的方法可以计算出另外三个理论数。3、计算，化如同吻合度检验那样计算/值，若左丸，则则接受Ho；若/，a,则才!绝Ho。与吻合度检验一样，理论数不得小于5,当理论数小于5时应使用另外的方法。4、确定自由度：因为每一行的各理论数受该行总数约束，每一列的各理论数受该列总数约束，所以自由度df=(r-1)(c-1)。上例的计算结果如下:有效无效口服Oi=58。2=4098Ti
10、=(98)(122)/193=61.95T2=(98)(71)/193=36.05注射O3=64O4=3195T3=(95)(122)/193=60.05T4=(95)(71)/193=34.95总数122711932JOiTi2二2-i=1Ti_2_2_2_25861.954036.0564-60.0531-34.95=+61.9536.0560.0534.95=0.2520.4330.2600.466=1.391H0:O一T=0,=0.05,df=2-12-1=1,o.05 = 3.841 ,2 :二2 0.05结论：接受H。，不同给药方式的治疗效果没有显著不同。本例的df=1应当矫正，矫正后的，值更小，不会影响结论，可以不再矫正。5、rxc列联表：上面的例子称为2X2列联表，对于行、列大于2的情况称为rxc列联表。其理论数的计算与2X2列联表相同：Tij=(i行总数)(j列总数)/总数。df=(r-1)(c-1)0推断过程也类似，不再举例。
