第四章协方差分析精选PPT.ppt

《第四章协方差分析精选PPT.ppt》由会员分享，可在线阅读，更多相关《第四章协方差分析精选PPT.ppt（34页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第四章 协方差分析第1页，此课件共34页哦 方方差差分分析析中中，所所接接触触到到的的各各种种处处理理多多数数都都是是人人为为控控制制的的。有有时时一一些些变变量量很很难难或或者者不不可可能能人人为为控控制制，对对于于这这种种情情况况则则不不能能用用第第三三章章所所述述的的方方差差分分析析方方法法，推推断处理之间的差异而应用协方差分析的方法做推断。断处理之间的差异而应用协方差分析的方法做推断。例如，在研究不同饲养条件下，动物的增重情况例如，在研究不同饲养条件下，动物的增重情况时，由于动物的原体重不同，如果只考虑饲料对增重时，由于动物的原体重不同，如果只考虑饲料对增重的作用而不考虑原体重时增重的

2、影响，显然是不全面的作用而不考虑原体重时增重的影响，显然是不全面的。动物的增重不仅与饲料有关，而且还与原体重有的。动物的增重不仅与饲料有关，而且还与原体重有关。一般来讲，原体重高的增童较多，原体重低的增关。一般来讲，原体重高的增童较多，原体重低的增重较少，增重与原体重之间有回归关系，为了得到正重较少，增重与原体重之间有回归关系，为了得到正确的结论，就要在排除原体重对增重的影响之后，推确的结论，就要在排除原体重对增重的影响之后，推断处理的效应。这就是本章所要讨论的协方差分析。断处理的效应。这就是本章所要讨论的协方差分析。第2页，此课件共34页哦在协方差分析中，通常将动物增重称为反应变量在协方差分

3、析中，通常将动物增重称为反应变量(responsevariable)Y，而与，而与Y 有线性回归关系的另一有线性回归关系的另一变量变量(原始体重原始体重)称为称为伴随变量伴随变量(concomitantvariable)或或协变协变量量(covariate)X。协方差分析就是通过反应变量与协变。协方差分析就是通过反应变量与协变量之间的线性关系，调整观察的反应变量。假若不做这量之间的线性关系，调整观察的反应变量。假若不做这种调整，由于协变量的存在会使误差平方和加大，其结种调整，由于协变量的存在会使误差平方和加大，其结果有可能检验不出反应变量在不同处理之间所存在的真果有可能检验不出反应变量在不同处

4、理之间所存在的真正差异。因此，正差异。因此，协方差分析协方差分析是一种调整无法控制又影响是一种调整无法控制又影响效应的变量的方差分析方法，是方差分析与回归分析的结效应的变量的方差分析方法，是方差分析与回归分析的结合。在上述动物重实验中，很难将全部实验动物都选为具合。在上述动物重实验中，很难将全部实验动物都选为具有相同体重的个体，动物的原体重是无法控制的变量，而有相同体重的个体，动物的原体重是无法控制的变量，而原体重对增重又有明显的影响，使用协方差分析，便可排原体重对增重又有明显的影响，使用协方差分析，便可排除原体重除原体重X对增重对增重Y的影响，使结果更为可靠的影响，使结果更为可靠。第3页，此

5、课件共34页哦第第一一节节 具具一一个个协协变变量量的的一一种种方方式式分分组组的的协协方方差差分分 析析（one-way classification with asinglecovariate）本本节节讨讨论论协协方方差差分分析析中中最最简简单单的的一一种种情情况况，即即以以只只有有一一个个协协变变量量的的单单因因素素实实验验设设计计的的协协方方差差分分析析为为例例，说说明明协协方方差差分分析析的的基基本本原原理理和和计计算算方方法法。单单因因素协方差分析的统计模型素协方差分析的统计模型是：是：第4页，此课件共34页哦其其中中yij是是第第i 次次处处理理所所得得到到的的反反应应变变量量的

6、的第第j 次次观观察察值值。c cij是是相相当当于于yij的的协协变变量量值值。c c是是c cij的的平平均均数数，m m是是总总平平均均数数，a ai是是第第i次次处处理理效效应应，b b是是yij在在c cij上上的的线线性性回回归归系系数数，e eij是是随随机机误误差差成成份份。做做协协方方差差分分析析，需需要要满满足足以以下下几几个个条条件件：e eij是是服服从从正正态态分分布布的的独独立立随随机机变变量量；b b0，即即yij与与c cij之之间间存存在在线线性性关关系系；各各处处理理的的回回归归系系数数都都相相同同；处处理理效效应应之之和和等等于于零零（a ai0）以以及及

7、协协变变量量不不受受处理效应的影响等。处理效应的影响等。第5页，此课件共34页哦 从（从（41）式中可以看出：协方差分析的模型是方）式中可以看出：协方差分析的模型是方差分析和回归分析线性模型的结合。差分析和回归分析线性模型的结合。a ai是单因素方差是单因素方差分析中的处理效应，分析中的处理效应，b b是回归分析中的回归系数。在是回归分析中的回归系数。在（41）式中的协变量是以（）式中的协变量是以（c cijc c）的形式而不是以）的形式而不是以c cij 的形式表示的，因此参量的形式表示的，因此参量m m是总的平均值。统计是总的平均值。统计模型的另一种形式是：模型的另一种形式是：其中其中m

8、m并不等于总平均值。在这个模型中，总平均值并不等于总平均值。在这个模型中，总平均值为为m mb b c c。经常使用的模型是（。经常使用的模型是（41）式。）式。第6页，此课件共34页哦 协方差分析需要计算以下一些量：协方差分析需要计算以下一些量：第7页，此课件共34页哦以上各式的符号：以上各式的符号：S、T 和和E 分别表示总的、处理的误差分别表示总的、处理的误差的平方和及交叉乘积和。它们之间的关系可用通式的平方和及交叉乘积和。它们之间的关系可用通式STE表示。仔细阅读并分析以上各式，弄清楚各式的意义，表示。仔细阅读并分析以上各式，弄清楚各式的意义，对下面的学习是有益的。对下面的学习是有益的

9、。第8页，此课件共34页哦协方差分析的核心协方差分析的核心，就是通过协变量调整反应变量。，就是通过协变量调整反应变量。下面讨论如何做这种调整。在统计模型（下面讨论如何做这种调整。在统计模型（41）中，）中，m m的的估计值为估计值为c c，b b的估计值为的估计值为b*，处理效应，处理效应a ai 的估计值的估计值为为c cic cb*（c cic c）。其中，）。其中，在这个模型中，误差平方和在这个模型中，误差平方和SSe为：为：具具a(n1)1自由度。自由度。第9页，此课件共34页哦实验误差均方由下式估计，实验误差均方由下式估计，若实验不存在处理效应，模型（若实验不存在处理效应，模型（41

10、）将变为：）将变为：具具an2自由度，其中自由度，其中S2XYSXX是由于是由于Y和和X的回归的回归所产生的平方和。所产生的平方和。第10页，此课件共34页哦如如果果实实验验本本身身存存在在处处理理效效应应，但但却却按按不不存存在在处处理理效效应应对对待待，这这时时所所计计算算出出来来的的误误差差平平方方和和SSe要要大大于于按按存存在在处处理理效效应应计计算算所所得得到到的的误误差差平平方方和和SSe。两两者者的的差差（SSeSSe）是是由由于于处处理理效效应应a ai 所所产产生生的的平平方方回回归归系系数数，具具a1自自由由度度。可可用用F 检检验验不不存存在在处处理理效效应的假设。应的

11、假设。若若FFa1,a(n1),a a，则接受，则接受H0：a ai0；若；若FFa1,a(n1),a a，则拒绝，则拒绝H0：a ai0。第11页，此课件共34页哦将以上结果列在表将以上结果列在表41中。中。表表41协方差分析调整的方差分析协方差分析调整的方差分析变差来源变差来源平方平方和和自自由由度度均均方方F回归回归处理处理误差误差S2XYSXXSSeSSe(SYYS2XYSXX)(EYYE2XYEXX)SSeEYYE2XYEXX1aa(n1)1(SSeSSe)(a1)MSeSSea(n1)1(SSeSSe)(a-1)MSe总和总和SYYan1第12页，此课件共34页哦表表41与方差分析

12、表基本上是一致的。所不同的与方差分析表基本上是一致的。所不同的是没一项平方和都是经过调整的。因此，是没一项平方和都是经过调整的。因此，协方差分析协方差分析又称为调整的方差分析又称为调整的方差分析。在变差来源一列中，总的变。在变差来源一列中，总的变差是由具差是由具an1自由度的自由度的SYY度量的；回归的变差由度量的；回归的变差由具具1自由度的平方和自由度的平方和S2XYSXX度量。假若不存在协变度量。假若不存在协变量，则量，则SXYSXXEXYEXX0。误差平方和将简化为。误差平方和将简化为EYY，处理平方和为，处理平方和为SYYEYYTYY，成为一种方式分，成为一种方式分组的方差分析。然而由

13、于存在协变量，我们必须通过组的方差分析。然而由于存在协变量，我们必须通过Y在在X上的回归，调整上的回归，调整SYY和和EYY（见表（见表41）。因为在）。因为在调整平方和时，用了另一个参量调整平方和时，用了另一个参量b，所以调整的误差，所以调整的误差平方和具平方和具a(n1)1自由度，而不是自由度，而不是a(n1)自由度。自由度。第13页，此课件共34页哦 通常将协方差分析结果纳成协方差分析表。通常将协方差分析结果纳成协方差分析表。表表42具一个协变量的一种方式分组实验的协方差分析表具一个协变量的一种方式分组实验的协方差分析表变差变差来源来源自由度自由度平方和与平方和与交叉乘积和交叉乘积和因因

14、回回归归所所做做的的调调整整XXYYY自由度自由度均均方方处理处理误差误差总和总和a1a(n1)an1TXXEXXSXXTXYEXYSXYTYYEYYSYYSSeEYYE2XYEXXSSeSYYS2XYSXXa(n1)1an2MSeSSea(n1)1调整调整的处的处理响理响SSeSSea1(SSeSSe)(a1)第14页，此课件共34页哦在协方差分析表中，除列入检验假设所需要在协方差分析表中，除列入检验假设所需要的处理效应平方和之外，还列入了全部平方和及的处理效应平方和之外，还列入了全部平方和及交叉乘积和。协方差分析的结果，不论零假设是交叉乘积和。协方差分析的结果，不论零假设是否可以接受，都需

15、对处理平方数否可以接受，都需对处理平方数c ci给予解释。由给予解释。由于于c ci包括处理效应和在协变量上的回归效应，因包括处理效应和在协变量上的回归效应，因此对平方数也要做相应的调整。调整的方法如下：此对平方数也要做相应的调整。调整的方法如下：第15页，此课件共34页哦根据模型（根据模型（42），协方差分析需满足以下要求：），协方差分析需满足以下要求：即：各处理的方差应具备齐性，它们都是从具有同即：各处理的方差应具备齐性，它们都是从具有同一方差的正态总体中的来的；个处理的回归系数一方差的正态总体中的来的；个处理的回归系数b bi均均等于等于b b以及反应变量与协变量之间的回归系数以及反应变

16、量与协变量之间的回归系数b b0。因此，。因此，在对一组数据做协方差分析时，首先要对以上各个条件做在对一组数据做协方差分析时，首先要对以上各个条件做检验。只有以上条件得到满足时，才能做协方差分析。检验。只有以上条件得到满足时，才能做协方差分析。第16页，此课件共34页哦第二节第二节 协方差分析的计算方法协方差分析的计算方法 例例41 比比较较三三种种猪猪饲饲料料A1，A2，A3对对猪猪增增重重的的影影响响，测测得得每每头头猪猪的的增增重重（Y）和和出出生生重重（X），数数据据列列在在表表43中中。问问三三种种饲饲料料对对猪猪增增重重是是否否有有显显著著不不同同的的效效果果？表表43不同饲料对猪

17、增重的影响不同饲料对猪增重的影响A1XY16851383116512761280169114841790X113.750Y181.750A2XY17971690181001895211032210619991894X218.625Y298.000A3XY228924912083239525100271023010532110X325.375Y396.875第17页，此课件共34页哦在在这这个个问问题题中中，若若不不考考虑虑出出生生重重，则则是是一一个个单单因因素素方方差差分分析析的的问问题题；若若不不同同饲饲料料的的增增重重效效果果没没有有显显著著差差异异，则则成成为为增增重重对对出出生生重重

18、的的一一元元回回归归问问题题。实实际际上上，在在研研究究动动物物增增重重问问题题时时，不不同同饲饲料料和和出出生生重重对对增增重重的的影影响响都都要要考考虑虑。因因此此，在在推推断断不不同同饲饲料料的的增增重重效效应应时时，为为了了排排除除出出生生重重的的影影响响，应应使使用用方方差差分分析析与与回回归归分分析析相相结结合合的的方方法法，即即以以协协方方差差分分析析的方法做推断。的方法做推断。首首先先检检验验e eij：NID(0,s s2)，b bib b及及b b0是是否否可可以以得得到到满足。检验的方法如下：满足。检验的方法如下：分别计算三种饲料饲养猪的出生重与增重间的分别计算三种饲料饲

19、养猪的出生重与增重间的回归系数回归系数bi，并列出回归方程，并列出回归方程 第18页，此课件共34页哦将三组数据合并（只需将三组数据的平方和与交将三组数据合并（只需将三组数据的平方和与交叉乘积合并即可），计算公共的回归系数叉乘积合并即可），计算公共的回归系数b*。用公共的回归系数代替用公共的回归系数代替bi，可以得到三条平行的回归，可以得到三条平行的回归线：线：第19页，此课件共34页哦 用每一处理的各自回归系数用每一处理的各自回归系数bi，计算剩余平方和：，计算剩余平方和：将三组剩余平方和相加，得到组内剩余平方和。将三组剩余平方和相加，得到组内剩余平方和。相应的自由度为相应的自由度为a(n2

20、)18。用公共回归系数用公共回归系数b*计算三条平行回归线的剩余平方和，计算三条平行回归线的剩余平方和，它们的和就是误差平方和它们的和就是误差平方和见公式见公式(49)(413)，相，相应的自由度为应的自由度为a(n1)120。第20页，此课件共34页哦检验方差齐性：检验方差齐性：中已计算出各处理的剩余平方和，中已计算出各处理的剩余平方和，各具各具8116自由度。各处理的均方分别为：自由度。各处理的均方分别为：检验方差齐性的一个简便方法是用两个差异最大的均检验方差齐性的一个简便方法是用两个差异最大的均方做方做F 检验。检验。F6,6,0.054.82，FF0.05，可可以以认认为为各各组组方方

21、差差具具备备齐齐性。性。第21页，此课件共34页哦 检检验验回回归归线线是是否否平平行行(b bib b)：在在中中已已经经计计算算出出组组内内剩剩余余平平方方和和SSeG和和用用公公共共回回归归系系数数b*计计算算得得到到的的误误差差平平方方和和SSe。SSeG完完全全是是由由随随机机因因素素造造成成的的：三三条条回回归归线线用用同同一一b*计计算算出出的的误误差差平平方方和和SSe，包包括括由由于于随随机机误误差差及及回回归归系系数数两两种种变变差差所所产产生生的的平平方方和和，因因而而回归系数平方和，回归系数平方和，回归系数自由度为误差自由度与组内误差自由度之差，回归系数自由度为误差自由

22、度与组内误差自由度之差，第22页，此课件共34页哦然然后后用用MSeG对对MS回回归归系系数数最最检检验验。若若两两者者差差异异不不显显著著，说说明明由由回回归归系系数数所所产产生生的的方方差差与与由由于于随随机机因因素素所所引引起起的的方方差差差差别别不不大大，即即可可证证明明原原来来的的三三条条回回归线是平行的。用归线是平行的。用F 检验，检验，代入数值，代入数值，F2,18,0.053.55，FF0.05，因此三条回归线是平行的。，因此三条回归线是平行的。第23页，此课件共34页哦 也也可可以以将将三三条条回回归归线线做做图图。从从图图上上观观察察它它们们是是否否近于平行。若接近平行的话

23、，也可以不做以上的检验。近于平行。若接近平行的话，也可以不做以上的检验。检验回归是否显著：利用方差分析做检验。方差检验回归是否显著：利用方差分析做检验。方差分析表如下：分析表如下：表表44检验回归显著性的方差分析表检验回归显著性的方差分析表 变差变差 来源来源平平 方方 和和 自自 由由 度度均均 方方F 回归回归 剩余剩余SSR=E2XY/EXX=1010.76SSe=EYYE2XY/EXX=227.615 1a(n1)1=20 MSR=SSR=1010.76MSe=SSe/a(n-1)-1=11.38 F=MSR/MSe=88.8 总和总和EYY=1238.375a(n-1)=21F1,2

24、0,0.01=8.1，FF0.01，回归是极显著的。这一步检验的许，回归是极显著的。这一步检验的许多计算，是在协方差分析过程中得出来的。所以这一步检验多计算，是在协方差分析过程中得出来的。所以这一步检验也可放在协方差分析的最后进行。也可放在协方差分析的最后进行。第24页，此课件共34页哦在在大大多多数数生生物物学学问问题题中中，以以上上几几点点要要求求基基本本上上都都可可以以满满足足。在在实实际际应应用用时时，只只要要根根据据bi所所做做出出的的回回归归线线是是否否平平行行。若若近近于于平平行行，则则可可不不必必做做繁繁琐琐的的检验，直接进行以下的分析。检验，直接进行以下的分析。将表将表43中

25、的数据编码，每一个中的数据编码，每一个Y 都减去都减去90，每，每一个一个X 都减去都减去20，列成表，列成表45。第25页，此课件共34页哦表表45协方差分析计算表协方差分析计算表X1X2X3Y1Y2Y3X1Y1X2Y2X3Y3 -5 -7 -9 -8 -8 -4 -6 -3 -3-4-2-2 1 2-1-2 2 4 0 3 5 7 10 12 -5-7-25-14-10 1-6 0 7 0 10 5 13 16 9 4 -1 1-7 5 10 12 15 20 25 49 225 112 80 -4 36 0 -210 -20 -10 13 32 -9 -8 -240 15 50 84 1

26、50 240 和和和和和和和和和的平方和的平方Xij平方的和平方的和-502500 344-11121 43 431849 347-184470 734 和和和的平方和的平方Yij平方的和平方的和 -6643561032 644096 696 553025 945 5311477 2673 和和523-235411041 第26页，此课件共34页哦计算下列各值：计算下列各值：第27页，此课件共34页哦由公式（由公式（415）求出：）求出：具具an2(3)(8)222自自由由度度。并并由由(413)式式计计算出算出具具a(n1)13(81)120自由度。自由度。第28页，此课件共34页哦 将上述

27、结果列成协方差分析表将上述结果列成协方差分析表(表表46)。表的最后。表的最后一行，是为了检验假设一行，是为了检验假设H0：a ai0所计算的平方和，所计算的平方和，具具a1312自由度。自由度。表表46协协方方差差分分析析表表变差来源变差来源自由度自由度平平 方方 和和 与与 交交 叉叉 乘乘 积积 调调 整整 平方和平方和调调 整整自由度自由度均均 方方FXXYY饲饲 料料 误误 差差 2 21545.250175.250 59.875 20.8751317.5831238.375227.615 2011.381 总总 和和23720.5001080.7502555.958934.8332

28、2调整饲料调整饲料707.2182353.609 31.07*第29页，此课件共34页哦检验不同饲料的增重效果，即检验检验不同饲料的增重效果，即检验H0：a ai0，根据，根据(416)式：式：F2,22,0.015.72，FF0.01，结论是不同饲料的增重效，结论是不同饲料的增重效果极显著。果极显著。第30页，此课件共34页哦公共回归系数，公共回归系数，回归系数的显著性检验，回归系数的显著性检验，H0：b b0，F1,20,0.018.1，FF0.01，结论是回归极显著。，结论是回归极显著。第31页，此课件共34页哦 最后还应计算出调整平均数，以便判断哪一种饲料最后还应计算出调整平均数，以便判断哪一种饲料最好。根据最好。根据(417)式：式：得到：得到：第32页，此课件共34页哦第33页，此课件共34页哦调整的平均数与未调整的平均数之间，存在极大的调整的平均数与未调整的平均数之间，存在极大的不同：调整前的顺序是不同：调整前的顺序是y2 y3 y1，第一种饲料的增，第一种饲料的增重效果最差。经过调整后的为重效果最差。经过调整后的为y2 y1 y3 ，第三种饲，第三种饲料的增重效果最差。这样的结论是真实的。由此可以料的增重效果最差。这样的结论是真实的。由此可以看出处理这类问题时，协方差分析是多么必要！看出处理这类问题时，协方差分析是多么必要！第34页，此课件共34页哦