七章卡平方测验.ppt

《七章卡平方测验.ppt》由会员分享，可在线阅读，更多相关《七章卡平方测验.ppt（69页珍藏版）》请在三一文库上搜索。

1、第七章 卡平方( )测验,第一节 卡平方( )的定义和分布 第二节 在方差同质性测验中的应用 第三节 适合性测验 第四节 独立性测验 第五节 的可加性和联合分析,第一节 卡平方( )的定义和分布,所谓 ，是指相互独立的多个正态离差平方值的总和，即：,(71),其中，yi 服从正态分布 ， 为标准正态离差。,yi不一定来自同一个正态总体，即 及 可以是不同正态分布的参数。若通常所研究的对象属同一个总体，则 ， ，从而,(72),抽样分布的密度函数为,累积分布函数为,分布的自由度为独立的正态离差的个数，此处 v =n ，其分布图形为一组具不同自由度 v 值的曲线(图7.1)。 值最小为0，最大为+

2、，因而在坐标轴的右面。 自由度小时呈偏态，随着自由度增加，偏度降低，至+时，呈对称分布。 该分布的平均数为 v ，方差为2v。,图7.1 不同自由度的 分布曲线,若所研究的总体 不知，而以样本 代替，则,(73),此时独立的正态离差个数为n1个，故 v =n1。,与u、t、F统计数的比较：,按定义 ，当只有1个正态离差时 ,，当s的自由度无限增大时 ，,此时 的 v =1。,，当 的自由度无限增大时 ，,v 为s12的自由度。,K.Pearson(1900)根据 的上述定义从属性性状的分布推导出用于次数资料(亦称计数资料)分析的 公式：,(74),上式中O为观察次数，E为理论次数，i=1，k为

3、计数资料的分组数，自由度为 v ，依分组数及其相互独立的程度决定，这种形式的 分布图形与图7.1相同。,值是多项 ui2 或 (OE)2/E 之和， 具有可加性。,第二节 在方差同质性测验中的应用,一、一个样本方差与给定总体方差比较的假设测验,可用来测验单个样本方差s2其所代表的总体方差和给定的总体方差值C是否有显著差异，简称为一个样本与给定总体方差的比较。 在作两尾测验时有 ，对 。其显著大于和小于C的值是 和 ,此时，H0在 显著水平上被否定。,例7.1 硫酸铵施于水田表层试验，得4个小区的稻谷产量为517、492、514、522(kg)，计得样本方差为175.6(kg)2。现要测验H0:

4、 对HA: ，采用显著水平 =0.05。,据 可算得：,查附表6，在 v =n1=3时， /2和(1 /2)水平的 临界值为： ， 。现 ，大于 ，在0.229.35范围外，符合H0的概率小于0.05，H0被否定。,结论:这一样本并非从 的总体中所抽取的。,若测验该样本总体方差是否小于某给定总体方差C，则作一尾测验，即H0： C对HA： C , 如果算得的 ，则否定H0，否则接受H0；这里应用分布的右边一尾。 如果测验其是否大于C，则H0： C对HA： C，若算得的 ，则否定H0；这是应用分布的左边一尾。,例7.2 试审查例7.1试验结果的总体方差是否真大于某一定值，如50(kg)2？,这里试

5、验的表面结果方差175.6(kg)2大于50(kg)2，要问其总体方差是否真正大，抑或并不大，甚至小于50(kg)2,测验假设H0： 50对HA： 50。 取5%为显著水平。 查附表6，这一测验的 临界值为 ，,而计算的 ，,因10.540.35，所以H0应被接受，即总体方差并不小于50(kg)2。,根据 ，可应用 分布由样本s2 给出一个总体 置信区间, 在此区间内包括有总体 的概率为( )， 即,(75),从而有：,(76A),已知 ，故(76A)又可记为：,(76B),例7.3 求例7.1资料总体 的95%置信限。,因为 ， ， ，且已知 s2 =175.6，故对总体方差 的95%置信限

6、的下限L1和上限L2为：,于是95%的置信限为：,注: 这一置信限并不对称，即从L1到s2的距离不等于s2到L2的距离。,利用置信限也可做显著性测验，例7.1中给定总体的 ，在56.32394.5范围外，故亦推断两者非同一总体。标准差的置信限可进而算出为：,即,本例因 较小，故方差置信限的区间甚大。,一般n30时，单个样本方差用 分布来测验和推断置信区间；n30时， 分布近似对称， 近似服从N(0，1)分布，因此，用u测验并进行区间估计。,二、几个样本方差的同质性测验,假定有3个或3个以上样本，每一样本均可估得一方差, 则由 可测验各样本方差是否来自相同方差总体的假设，这称为方差的同质性测验(

7、 test for homogeneity among variances )，可写为H0： (k为样本数) 对HA： 不全相等。这一测验方法由Bartlett氏(1937)提出，故又称为Bartlett测验( Bartlett test )，是一种近似的 测验。,假如有k个独立的方差估计值：, ,各具 个自由度，那么合并的方差 为：,(77),由此，Bartlett 值为：,(78),(79),上式的 ，ni为样本容量，而C为矫正数：,(710),如采用常用对数，则(79)可写为,(711),上述(78)如不用C进行矫正，亦近似地作 分布，具有 ； 若所得 值不显著，则不必再作矫正，应接受H

8、0； 若 值与 接近，应作矫正。 如果算得的 值 ，便否定H0，表明这些样本所属总体方差不是同质的。,例7.4 假定有3个样本方差s12=4.2, s22=6.0, s32=3.1，各具有自由度 ， ， ，试测验其是否同质。,假设H0： 对HA：3个方差不全相等(这里的HA不能用不等号表示，因为如H0被否定，只能推论3者不相等而并不能确定属于 、 、 等情况的哪一种)。 然后，在表7.1进行同质性测验的计算：,表7.1 3个方差同质性测验的计算,由表7.1可得：,查附表6，当时 ， 0.744的概率在0.500.75之间，符合H0的概率不小，因此说明本例的3个方差估计值是同质性的。,实际应用上

9、本例可不需再作C矫正，因为 =27.9496027.14452=0.80508明显很小，直观已可判断不会显著。,第三节 适合性测验,一、适合性 测验的方法 二、各种遗传分离比例的适合性测验 三、次数分布的适合性测验,一、适合性 测验的方法,例:玉米花粉粒中形成淀粉粒或糊精是一对相对性状。淀粉粒遇碘呈蓝色反应，因而可以用碘试法直接观察花粉粒的分离现象。某项实验观察淀粉质与非淀粉质玉米杂交的F1代花粉粒，经碘处理后有3437粒呈蓝色反应，3482粒呈非蓝色反应。根据遗传学理论可假设玉米花粉粒碘反应为11，由此可以计得3437+3482=6916粒花粉中，蓝色反应与非蓝色反应的理论次数应各为3459

10、.5粒。设以O代表观察次数，E代表理论次数，可将上列结果列成表7.2。,表7.2 玉米花粉粒碘反应观察次数与理论次数,此处要推论是否符合11分离，只要看观察次数与理论次数是否一致，故可用 测验，可分为四个步骤：,（1）设立无效假设，即假设观察次数与理论次数的差异由抽样误差所引起，即H0：花粉粒碘反应比例为11与HA：花粉粒碘反应比例不成11。,（2）确定显著水平 =0.05。,（3）在无效假设为正确的假定下，计算超过观察 值的概率，这可由 计得 值后，按自由度查附表6得到。试验观察的 值愈大，观察次数与理论次数之间相差程度也愈大，两者相符的概率就愈小。,（4）依所得概率值的大小，接受或否定无效

11、假设 在实际应用时，往往并不需要计算具体的概率值。 若实得 时，则H0发生的概率小于等于 ，属小概率事件，H0便被否定； 若实得 时，则H0被接受。,例如表7.2资料， 查附表6，当 时 =3.84 ，实得 =0.2926小于 ，所以接受H0。即认为观察次数和理论次数相符，接受该玉米F1代花粉粒碘反应比率为11的假设。,然而按 的定义 分布是连续性的，而次数资料则是间断性的。由间断性资料算得的 值有偏大的趋势(尤其在 时)，需作连续性矫正。其方法是：在度量观察次数相对于理论次数的偏差时，将各偏差的绝对值都减1/2，即|OE|1/2。矫正后的 用 表示，即,(712),如表7.2资料的 值为：,

12、=0.2798仍然小于 =3.84，结论与前相同。,这是因样本较大，故 与 值的相差不大。,一般 的样本，尤其是小样本，在计算 值时必须作连续性矫正，否则所得 值偏大，容易达到显著水平。对 2的样本，都可以不作连续性矫正。,当 30时， 分布已近于对称，而 的分布是正态的，具平均数 和标准差1。 因而，当 30时可采用正态离差u测验代替 测验 ，即,如u1.64，即表示实得 值有显著性。,二、各种遗传分离比例的适合性测验,例7.5 大豆花色一对等位基因的遗传研究，在F2获得表7.3所列分离株数。问这一资料的实际观察比例是否符合于31的理论比值。,表7.3 大豆花色一对等位基因遗传的适合性测验,

13、H0：大豆花色F2分离符合31比率；HA：不符合31比率。 显著水平 =0.05。 由于该资料只有k=2组， ，故在计算 值时需作连续性矫正。,由 可得：,查附表6， 。现 故应接受H0，说明大豆花色这对性状是符合31比率，即符合一对等位基因的表型分离比例。,分离比例一类的适合性测验计算 时，也可以不经过计算理论次数，而直接得出,(713),其中，A和a分别为显性组和隐性组的实际观察次数；n=A+a，即总次数。本例资料代入(713)有：,与(712)算得的 值相同。,对于仅划分为两组(如显性与隐性)的资料，如测验其与某种理论比率的适合性，则其 值皆可用类似(713)的简式求出。这些简式列于表7

14、.4。,表7.4 测验两组资料与某种理论比率符合度的 值公式,例7.6 两对等位基因遗传试验，如基因为独立分配，则F2代的四种表现型在理论上应有9331的比率。有一水稻遗传试验，以稃尖有色非糯品种与稃尖无色糯性品种杂交，其F2代得表7.5结果。试检查实际结果是否符合9331的理论比率。,表7.5 F2代表型的观察次数和根据9331算出的理论次数,首先，按9331的理论比率算得各种表现型的理论次数E， 如稃尖有色非糯稻 E=743(9/16)=417.94， 稃尖有色糯稻 E=743(3/16)=139.31，。 H0：稃尖和糯性性状在F2的分离符合9331； HA：不符合9331。 显著水平：

15、 =0.05。 然后计算 值,因本例共有k=4组，故 =k-1=3。查附表6， ,现实得 ,所以否定H0，接受HA，即该水稻稃尖和糯性性状在F2的实际结果不符合9331的理论比率。,这一情况表明，该两对等位基因并非独立遗传，而可能为连锁遗传。,测验实际结果与9331理论比率的适合性，也可不经过计算理论次数而直接用以下简式,(714),上式中的a1、a2、a3、a4分别为9331比率中各项表现型的实际观察次数，n为总次数。,如本例，可由(714)算得：,前面的 =92.696，与此 =92.706略有差异，系前者有较大计算误差之故。,实际资料多于两组的 值通式则为：,(715),上式的mi为各项

16、理论比率，ai为其对应的观察次数。,如本例，亦可由(715)算得,与此一致。,三、次数分布的适合性测验,适合性测验还经常用来测验试验数据的次数分布是否和某种理论分布(如二项分布、正态分布等)相符，以推断实际的次数分布究竟属于哪一种分布类型。,例7.7 在大豆品种Richland田间考察单株粒重的变异是否符合正态分布。考查数据归成次数分布表列于表7.6，组距为5g，该分布的次数n、平均数 、标准差s 均列于表基部。,表7.6 大豆单株粒重观察分布与理论正态分布的适合性测验 (摘自Steel and Torrie，1980)(单位：g),测验的假设为H0：观察分布符合理论分布，HA：观察分布不符合

17、理论分布。 按理论分布计算出各组的理论次数(E)，此例中正态分布下的理论次数可先计算出各组限的正态离差及其理论频率(P)，乘以总观察次数(n)便得到各组的理论次数。 例如第1组,第2组 P(5.5y10.5)=P(2.065u1.674) =0.04710.0195=0.0276,相应的理论次数E，第一组为0.0195229=4.5； 第二组为0.0276229=6.3,其他各组按同法计算后均列入表7.6。,自由度 =1412=11，因扣去组数的自由度1个，估计2个参数 和 的自由度2个。 查附表6， 为11时 =10.47的概率P在0.250.50范围内，观察分布与理论分布无显著差异，因而接

18、受H0，说明大豆单株粒重的分布符合正态分布。,用于进行次数分布的适合性测验时有一定的近似性，为使这类测验更确切，一般应注意以下几点：,(1)总观察次数n应较大，一般不少于50。 (2)分组数最好在5组以上。 (3)每组理论次数不宜太少，至少为5，尤其首尾各组。若组理论次数少于5，最好将相邻组的次数合并为一组。但Cochran认为头尾二组最小理论次数在0.5或1时也可不合并。,例7.7中尾端三组理论次数均较少，若将后三组合并，则,P(55.5y70.5)=P(1.841u3.013) =0.9887-0.9671=0.0316,该组理论次数为 0.0316229=7.27, (OE)2/E=(1

19、07.27)2/7.27=1.025 =123=9，查附表6， 为9时 =10.425的概率P在0.250.50范围内，结论同前。,第四节 独立性测验,应用 进行独立性测验的无效假设是： H0：两个变数相互独立，对HA：两个变数彼此相关。,计算过程: (1)将所得次数资料按两个变数作两向分组，排列成相依表； (2)根据两个变数相互独立的假设，算出每一组格的理论次数； (3)由 算得 值。,这个 的自由度随两个变数各自的分组数而不同，设横行分r组，纵行分c组，则 =(r1)(c1)。 当观察的 时，便接受H0，即两个变数相互独立；当观察的 时，便否定H0，接受HA，即两个变数相关。,独立性测验方

20、法的各种类型 一、22表的独立性测验 二、2C表的独立性测验 三、rc表的独立性测验,一、22表的独立性测验,一、22表的独立性测验,22相依表是指横行和纵行皆分为两组的资料。在作独立性测验时，其 =(21)(21)=1，故计算 值时需作连续性矫正。,例7.8 调查经过种子灭菌处理与未经种子灭菌处理的小麦发生散黑穗病的穗数，得相依表7.7，试分析种子灭菌与否和散黑穗病穗多少是否有关。,表7.7 防治小麦散黑穗病的观察结果,假设H0：两变数相互独立，即种子灭菌与否和散黑穗病病穗多少无关；HA：两变数彼此相关。 显著水平 =0.05。 根据两变数相互独立的假定，算得各组格的理论次数。 如种子灭菌项

21、的发病穗数O1=26，其理论次数E1=(21076)/460=34.7，即该组格的横行总和乘以纵行总和再除以观察总次数(下同)；同样可算得 O2=50 的 E2=(25076)/460=41.3； O3=184的E3=(210384)/460=175.3； O4=200的E4=(250384)/460=208.7。 以上各个E值填于表7.7括号内。,以上各个E值代入 有,这里 =(21)(21)=1，查附表6， ，现实得 ，故P0.05，应否定H0。即种子灭菌与否和散黑穗病发病高低有相关，种子灭菌对防治小麦散黑穗病有一定效果。,22表的独立性测验也可不经过计算理论次数而直接得到 值。22表的一

22、般化形式如表7.8。,按表中的符号,表7.8 22表的一般化形式,(716),如本例各观察次数代入(716)可得：,与前面结果相同。,二、2C表的独立性测验,2C表是指横行分为两组，纵行分为C3组的相依表资料。 在作独立性测验时，其 =(21)(c1)=c1。由于c3,故不需作连续性矫正。,例7.9 进行大豆等位酶Aph的电泳分析，193份野生大豆、223份栽培大豆等位基因型的次数列于表7.9，试分析大豆Aph等位酶的等位基因型频率是否因物种而不同。,表7.9 野生大豆和栽培大豆Aph等位酶的等位基因型次数分布,H0：等位基因型频率与物种无关；HA：两者有关，不同物种等位基因型频率不同。 显著

23、水平 =0.05。 根据H0算得各观察次数的相应理论次数： 如观察次数29的E=(19351)/416=23.66， 观察次数22的E=(22351)/416=27.34，； 将其填于表7.9的括号内。,再代入 可得：,此处 =(21)(31)=2。查附表6， ，现 ，P0.05，应否定H0，接受HA。即不同物种Aph等位基因型频率有显著相关，或者说不同物种的Aph等位基因型频率有显著差别。,2C表独立性测验的值，也可直接由下式得到。,2C表的一般化形式如表7.10。,(717),( i=1，2，3，c ),表7.10 2C表的一般化形式,三、rc表的独立性测验,若横行分r 组，纵行分c 组，

24、且r3，c3，则为rc相依表。 对rc表作独立性测验时，其 =(r1)(c1)，计求 值不需要连续性矫正。,例7.10 表7.11为不同灌溉方式下水稻叶片衰老情况的调查资料。试测验稻叶衰老情况是否与灌溉方式有关。,表7.11 水稻在不同灌溉方式下叶片的衰老情况,H0：稻叶衰老情况与灌溉方式无关；HA：稻叶衰老情况与灌溉方式有关。 取 =0.05。 根据H0的假定，计算各组格观察次数的相应理论次数：如与146相应的E=(481160)/547=140.69， 与183相应的E=(481205)/547=180.26， 所得结果填于表7.11括号内。,根据 可得,本例 =(31)(31)=4，查附

25、表6， ，现 ,P0.05，故应接受H0，即不同灌溉方式对水稻叶片的衰老情况没有显著影响。,rc表的一般化形式如表7.12。,表7.12 rc表的一般化形式,由表7.12直接计算值的公式：,(718),( i=1，2，r；j=1，2，3，c ),将表7.11资料，代入(718)有,第五节 的可加性和联合分析,例7.11 表7.13给出三个大豆组合F3家系世代对豆秆黑潜蝇抗性家系与感性家系的分离数据，每一家系由1个F2单株衍生，抗性家系中包括有全抗家系及抗感分离的家系。经对三个组合分别的 测验，均符合3抗1感理论分离比例。现要求进一步检测三组合综合起来是否符合31分离比例，三组合间是否一致符合3

26、1分离比例，或三组合是否具同质性。,表7.13 三个大豆组合F3家系世代对豆秆黑潜蝇抗性的分离数据(理论分离比为3抗1感）,H0：三组合综合起来符合3抗1感分离比例，HA：综合群体不符合31分离比例；及H0：三组合的分离比表现同质，一致为31，HA：三组合分离比例不同质。,要测验上列假设，必须计算出相应的 值。表7.13中列出有多种 值。,（1）各组合分别的 及 已用于测验各组合与理论分离比例31的相符性。这里不仅列出 值用于各测验；同时列出 值，因为 不具可加性，只有 值具有可加性。,（2）三个组合综合为一群体时的 值，或称为 =3.06，亦具1个自由度。这一值可用以测验第一个无效假设，根据

27、其概率为0.050.10，可推论三合一的群体总的分离比例亦符合31。 （3）三组合各 的总和 =3.66，具有3个自由度。若将这3个自由度分解，1个归属于三组合间的共性，2个归属于三组合间的个性，它们相应的 值为 和 。 已在（2）中进行过测验，剩下 具2个自由度可用以测验第二个无效假设，三个组合的同质性。此处 =0.60, 时P=0.500.75。说明符合同质性假设的概率甚大，接受此假设，因而三个组合表现一致的31分离比例是确实的。,根据表7.13的数据，三组合的亲本表现确实的抗、感差异，F2衍生的F3家系表现出抗性为显性并一致符合3抗(抗+感）1感的家系间分离比例，因而可推论大豆对豆秆黑潜蝇的抗性是由1对显性基因控制的，组合间表现出一致的结果。,本例中因试验结果很一致，因而引出了共同的结论。若各个 的结果出入较大， 与个别组合的结果不一致， 表现出显著性，那么将着重分析各组合间的非同质性及各组合的特异性。,