生物统计学：《生物统计附试验设计》课程教材（第三版）第十章协方差分析.doc_大学文库

第十章协方差分析第一节协方差分析的意义协方差分析有二个意义,一是对试验进行统计控制,二是对协方差组分进行估计,现分述如下。、对试验进行统计控制为了提高试验的精确性和准确性,对处理以外的一切条件都需要采取有效措施严加控制,使它们在各处理间尽量一致,这叫试验控制。但在有些情况下,即使作出很大努力也难以使试验控制达到预期目的。例如:研究几种配合饲料对猪的增重效果,希望试验仔猪的初始重相同,因为仔猪的初始重不同,将影响到猪的增重。经研究发现:增重与初始重之间存在线性回归关系。但是,在实际试验中很难满足试验仔猪初始重相同这一要求。这时可利用仔猪的初始重(记为x)与其增重(记为y)的回归关系,将仔猪增重都矫正为初始重相同时的增重,于是初始重不同对仔猪增重的影响就消除了。由于矫正后的增重是应用统计方法将初始重控制一致而得到的,故叫统计控制。统计控制是试验控制的一种辅助手段。经过这种矫正, 试验误差将减小,对试验处理效应估计更为准确。若y的变异主要由x的不同造成(处理没有显著效应),则各矫正后的y’间将没有显著差异(但原y的差异可能是显著的)。若y的变异除掉x不同的影响外,尚存在不同处理的显著效应,则可期望各y’间将有显著差异(但原y间差异可能是不显著的)。此外,矫正后的y’和原y的大小次序也常不一致。所以,处理平均数的回归矫正和矫正平均数的显著性检验,能够提高试验的准确性和精确性,从而更真实地反映试验实际。这种将回归分析与方差分析结合在一起,对试验数据进行分析的方法,叫做协方差分析( analysis of covariance) 、估计协方差组分在第八章曾介绍过表示两个相关变量线性相关性质与程度的相关系数的计算公式若将公式右端的分子分母同除以自由度(m-1),得 ∑(x-x)y-y)(n-1) (10-1) ∑(x-x) (n-1) 其中 ∑ 是x的均方MS,它是x方差G2的无偏估计量:

197 第十章协方差分析第一节协方差分析的意义协方差分析有二个意义，一是对试验进行统计控制，二是对协方差组分进行估计，现分述如下。一、对试验进行统计控制为了提高试验的精确性和准确性，对处理以外的一切条件都需要采取有效措施严加控制，使它们在各处理间尽量一致，这叫试验控制。但在有些情况下，即使作出很大努力也难以使试验控制达到预期目的。例如：研究几种配合饲料对猪的增重效果，希望试验仔猪的初始重相同，因为仔猪的初始重不同，将影响到猪的增重。经研究发现：增重与初始重之间存在线性回归关系。但是，在实际试验中很难满足试验仔猪初始重相同这一要求。这时可利用仔猪的初始重(记为x)与其增重(记为y)的回归关系，将仔猪增重都矫正为初始重相同时的增重，于是初始重不同对仔猪增重的影响就消除了。由于矫正后的增重是应用统计方法将初始重控制一致而得到的，故叫统计控制。统计控制是试验控制的一种辅助手段。经过这种矫正，试验误差将减小，对试验处理效应估计更为准确。若y的变异主要由x的不同造成(处理没有显著效应)，则各矫正后的 y  间将没有显著差异(但原y间的差异可能是显著的)。若y的变异除掉x不同的影响外，尚存在不同处理的显著效应，则可期望各 y  间将有显著差异(但原y间差异可能是不显著的)。此外，矫正后的 y  和原y的大小次序也常不一致。所以，处理平均数的回归矫正和矫正平均数的显著性检验，能够提高试验的准确性和精确性，从而更真实地反映试验实际。这种将回归分析与方差分析结合在一起，对试验数据进行分析的方法，叫做协方差分析(analysis of covariance)。二、估计协方差组分在第八章曾介绍过表示两个相关变量线性相关性质与程度的相关系数的计算公式：    − − − − = 2 2 ( ) ( ) ( )( ) x x y y x x y y r 若将公式右端的分子分母同除以自由度(n-1)，得       − −       − − − − − =    ( 1) ( ) ( 1) ( ) ( )( ) /( 1) 2 2 n y y n x x x x y y n r （10-1）其中 1 ( ) 2 −  − n x x 是x的均方MSx，它是x的方差 2  x 的无偏估计量；

198 1 ( ) 2 −  − n y y 是y的均方MSy，它是y的方差 2  y 的无偏估计量； 1 ( )( ) −  − − n x x y y 称为x与y的平均的离均差的乘积和，简称均积，记为MPxy，即 MPxy= 1 ( )( ) −  − − n x x y y = 1 ( )( ) −  −   n n x y xy （10-2）与均积相应的总体参数叫协方差（covariance），记为COV(x,y)或  xy 。统计学证明了，均积MPxy是总体协方差COV(x,y)的无偏估计量，即 EMPxy= COV(x,y)。于是，样本相关系数r可用均方MSx、MSy，均积MPxy表示为： x y xy MS MS MP r = （10-3）相应的总体相关系数  可用x与y的总体标准差  x 、 y ，总体协方差COV(x,y)或  xy 表示如下： x y xy x y COV x y       = = ( , ) （10-4）均积与均方具有相似的形式，也有相似的性质。在方差分析中，一个变量的总平方和与自由度可按变异来源进行剖分，从而求得相应的均方。统计学已证明：两个变量的总乘积和与自由度也可按变异来源进行剖分而获得相应的均积。这种把两个变量的总乘积和与自由度按变异来源进行剖分并获得获得相应均积的方法亦称为协方差分析。在随机模型的方差分析中，根据均方MS和期望均方EMS的关系，可以得到不同变异来源的方差组分的估计值。同样，在随机模型的协方差分析中，根据均积MP和期望均积EMP的关系，可得到不同变异来源的协方差组分的估计值。有了这些估计值，就可进行相应的总体相关分析。这些分析在遗传、育种和生态、环保的研究上是很有用处的。由于篇幅限制，本章只介绍对试验进行统控制的协方差分析。第二节单因素试验资料的协方差分析设有k个处理、n次重复的双变量试验资料，每处理组内皆有n对观测值x、y，则该资料为具kn对x、y观测值的单向分组资料，其数据一般模式如表10—1所示。表10—1 kn对观测值x、y的单向分组资料的一般形式处理处理1 处理2 … 处理i … 处理k 观测指标 x y x y … x y … x y 观测值 xij、yij (i=1,2,…k j=1,2,…n) x11 x12 … x1j … x1n y11 y12 … y1j … y1n x21 x22 … x2j … x2n y21 y22 … y2j … y2n … … … … … … xi1 xi2 … xij … xin yi1 yi2 … yij … yin … … … … … xk1 xk2 … xkj … xkn yk1 yk2 … ykj … ykn 总和 x1. y1. x2. y2. … xi. yi. … xk. yk. 平均数 . 1 x . 1 y . 2 x . 2 y … . i x . i y … . k x . k y

199 表10—1的x和y变量的自由度和平方和的剖分参见单因素试验资料的方差分析方法一节。其乘积和的剖分则为：总变异的乘积和 SPT 是 ij x 与 x.. 和 ij y 与 y.. 的离均差乘积之和，即： k n x y SP x x y y x y k i n j i j i j k i n j T i j i j .. .. ( ..)( ..) 1 1 1 1 =  − − =  − = = = = (10-5) T df =kn-1 其中， k n y y k n x x x y y x k i i k i i .. , .. .. .. ., .. ., .. 1 1 = = = = = = 。处理间的乘积和 SPt 是 . i x 与 x.. 和 . i y 与 y.. 的离均差乘积之和乘以n，即：   = = = − − = − k i i i k i t i i i i k n x y x y n SP n x x y y 1 1 . . . . 1 ( . ..)( . ..) (10-6) dft = k −1 处理内的乘积和 SPe 是 ij x 与 . i x 和 ij y 与 . i y 的离均差乘积之和，即：    = = = = = = − − = − = − k i n j k i T t k i i i n j e i j i i j i i j i j x y SP SP n SP x x y y x y 1 1 1 1 1 . . 1 ( .)( .) (10-7) e df =k(n-1) 以上是各处理重复数n相等时的计算公式，若各处理重复数n不相等，分别为n1、n2、…、 nk，其和为 = k i i n 1 ，则各项乘积和与自由度的计算公式为：   = = = = − k i n j k i i i i T ij ij i n x y SP x y 1 1 1 . . T df == k i i n 1 -1 (10-8) = = + + + − k i i k k k t n x y n x y n x y n x y SP 1 2 2 2 1 1 1 . . .. .. ... . . . . dft = k −1 = = = k i n j e ij ij i SP x y 1 1 -       + + + k k k n x y n x y n x y . . ... . . . . 2 2 2 1 1 1 =SPT-SPt e df == k i i n 1 -k =dfT-dft (10-9) 有了上述SP和df，再加上x和y的相应SS，就可进行协方差分析。【例10.1】为了寻找一种较好的哺乳仔猪食欲增进剂，以增进食欲，提高断奶重，对哺乳仔猪做了以下试验：试验设对照、配方1、配方2、配方3共四个处理，重复12次，选择初始条件尽量相近的长白种母猪的哺乳仔猪48头，完全随机分为4组进行试验，结果见表 10—2，试作分析。此例， .. . . . . 1 2 3 4 x = x + x + x + x =18.25+15.40+15.65+13.85=63.15 .. . . . . 1 2 3 4 y = y + y + y + y =141.80+130.10+144.80+133.80=550.50

1、误差项回归关系的分析误差项回归关系分析的意义是要从剔除处理间差异的影响的误差变异中找出50日龄重)与初生重(x)之间是否存在线性回归关系。计算出误差项的回归系数并对线性回归关系进行显著性检验,若显著则说明两者间存在回归关系。这时就可应用线性回归关系来校正y值(50日龄重)以消去仔猪初生重(x)不同对它的影响。然后根据校正后的y值(校正50日龄重)来进行方差分析。如线性回归关系不显著,则无需继续进行分回归分析的步骤如下 (1)计算误差项回归系数,回归平方和,离回归平方和与相应的自由度从误差项的平方和与乘积和求误差项回归系数: 7.1848 (10-10) 0.92 误差项回归平方和与自由度 47.49 (10-11) 0.92 误差项离回归平方和与自由度 Sar)-SS)=8508-4749=3759 dre)=dfely-dRle)=44-1=43 (2)检验回归关系的显著性(表10-5) 表105哺乳仔猪50日龄重与初生重的回归关系显著性检验表变异来源 F F 误差回归 4749 4749 54.32** 7255 误差离回归 3759 43 0.874 误差总和 85.08 F检验表明,误差项回归关系极显著,表明哺乳仔猪50日龄重与初生重间存在极显著的线性回归关系。因此,可以利用线性回归关系来校正y,并对校正后的y进行方差分析 2、对校正后的50日龄重作方差分析 (1)求校正后的50日龄重的各项平方和及自由度利用线性回归关系对50日龄重作校正,并由校正后的50日龄重计算各项平方和是相当麻烦的,统计学已证明,校正后的总平方和、误差平方和及自由度等于其相应变异项的离回归平方和及自由度,因此,其各项平方和及自由度可直接由下述公式计算 ①校正50日龄重的总平方和与自由度,即总离回归平方和与自由度 8.25 175 dr=dr()-可fR(y)=47-1=46 ②校正50日龄重的误差项平方和与自由度,即误差离回归平方和与自由度 SSe=SSe(r)-SSR(e)=SSelr)-" 0.92 dJ=de()-d(R)=441=43 上述回归自由度均为1,因仅有一个自变量x

202 1、误差项回归关系的分析误差项回归关系分析的意义是要从剔除处理间差异的影响的误差变异中找出50日龄重(y)与初生重(x)之间是否存在线性回归关系。计算出误差项的回归系数并对线性回归关系进行显著性检验，若显著则说明两者间存在回归关系。这时就可应用线性回归关系来校正y值(50日龄重)以消去仔猪初生重(x)不同对它的影响。然后根据校正后的y值(校正50日龄重)来进行方差分析。如线性回归关系不显著，则无需继续进行分析。回归分析的步骤如下： (1)计算误差项回归系数，回归平方和，离回归平方和与相应的自由度从误差项的平方和与乘积和求误差项回归系数： 7.1848 0.92 6.61 ( ) ( ) = = = e x e yx e SS SP b (10-10) 误差项回归平方和与自由度 47.49 0.92 6.612 ( ) 2 ( ) = = = e x e R e SS SP SS (10-11) dfR(e)=1 误差项离回归平方和与自由度 SSr(e) = SSe( y) - SSR(e) =85.08-47.49=37.59 (10-12) dfr(e) = df e( y) - df R(e) =44-1=43 (2)检验回归关系的显著性(表10—5) 表10—5 哺乳仔猪50日龄重与初生重的回归关系显著性检验表变异来源 SS df MS F F0.01 误差回归 47.49 1 47.49 54.32** 7.255 误差离回归 37.59 43 0.8742 误差总和 85.08 44 F检验表明，误差项回归关系极显著，表明哺乳仔猪50 日龄重与初生重间存在极显著的线性回归关系。因此，可以利用线性回归关系来校正y，并对校正后的y进行方差分析。 2、对校正后的50日龄重作方差分析 (1)求校正后的50日龄重的各项平方和及自由度利用线性回归关系对50日龄重作校正，并由校正后的50日龄重计算各项平方和是相当麻烦的，统计学已证明，校正后的总平方和、误差平方和及自由度等于其相应变异项的离回归平方和及自由度，因此，其各项平方和及自由度可直接由下述公式计算。 ①校正50日龄重的总平方和与自由度，即总离回归平方和与自由度 57.85 1.75 8.25 96.76 2 ( ) 2 ( ) ( ) ( ) ' = − = − = − = T x T T T y R y T y SS SP SS SS SS SS (10-13) ' T df = dfT ( y) - df R( y) =47-1=46 ②校正50日龄重的误差项平方和与自由度，即误差离回归平方和与自由度 37.59 0.92 6.61 85.08 2 ( ) 2 ( ) ( ) ( ) ' = − = − = − = e x e e e y R e e y SS SP SS SS SS SS (10-14) ' e df = df e( y) - df e(R) =44-1=43 上述回归自由度均为1，因仅有一个自变量x

203 ③校正50日龄重的处理间平方和与自由度 SSt SST SSe  =  −  =57.87-37.59=20.28 (10-15) t T e df  = df  − df  =k-1=4-1=3 (2)列出协方差分析表，对校正后的50日龄重进行方差分析(表10—6) 查F表： F0.01(3,43) =4.275(由线性内插法计算)，由于F=7.63＞ F0.01(3,43) ，P＜0.01，表明对于校正后的50日龄重不同食欲添加剂配方间存在极显著的差异。故须进一步检验不同处理间的差异显著性，即进行多重比较。表10—6 表10-2资料的协方差分析表变异来源 df SSx SSy SPxy b 校正50日龄重的方差分析 F df' SS' MS 处理间(t) 3 0.83 11.68 1.64 机误(e) 44 0.92 85.08 6.61 7.1848 43 37.59 0.8742 总和(T) 47 1.75 96.76 8.25 46 57.87 校正处理间 3 20.28 6.76 7.63** 3、根据线性回归关系计算各处理的校正50日龄平均重误差项的回归系数 byx(e) 表示初生重对50日龄重影响的性质和程度，且不包含处理间差异的影响，于是可用 byx(e) 根据平均初生重的不同来校正每一处理的50日龄平均重。校正50 日龄平均重计算公式如下： . . ( . ..) ( ) y y b x x i  = i − yx e i − (10-16) 公式中： . i y  为第i处理校正50日龄平均重； . i y 为第i处理实际50日龄平均重(见表10—2)； . i x 为第i处理实际平均初生重(见表10—2)； x.. 为全试验的平均数， 1.3156 48 .. 63.15 .. = = = kn x x byx(e) 为误差回归系数， byx(e) =7.1848 将所需要的各数值代入(10—16)式中，即可计算出各处理的校正50日龄平均重(见表 10—7)。表10—7 各处理的校正50日龄平均重计算表处理 . i x - x.. byx(e) ( . i x - x.. ) 实际50日龄平均重校正50日龄平均重 . i y - byx(e) ( . i x - x.. ) 对照 1.52-1.3156=0.2044 7.1848×0.2044=1.4686 11.82 11.82-1.1686=10.3514 配方1 1.28-1.3156=-0.0356 7.1848×（-0.0356）=-0.2588 10.84 10.84+0.2558=12.0758 配方2 1.30-1.3156=-0.0156 7.1848×（-0.0156）=-0.1121 12.07 12.07+0.1121=12.1821 配方3 1.15-1.3156=-0.1656 7.1848×（-0.1656）=-1.1898 11.15 11.15+1.1898=12.3398 4、各处理校正50日龄平均重间的多重比较各处理校正50日龄平均重间的多重比较，即各种食欲添加剂的效果比较。 (1)t检验检验两个处理校正平均数间的差异显著性，可应用t检验法： . . . . i j y y i j S y y t  −   −  = (10-17)

204         −  =  + −  ( ) 2 . . 2 ( . .) e x i j e y y SS x x n S MS i j (10-18) 式中， . . i j y  − y  为两个处理校正平均数间的差异；  −  . . i j y y S 为两个处理校正平均数差数标准误； MSe  为误差离回归均方； n为各处理的重复数； . i x 为处理i的x变量的平均数； . j x 为处理j的x变量的平均数； SSe(x)为x变量的误差平方和例如，检验食欲添加剂配方1与对照校正50日龄平均重间的差异显著性： . . 1 2 y  − y  =10.3514-12.0758=-1.7244 MSe  =37.59/43=0.8742 n=12 . 1 x =1.52, . 2 x =1.28, SSe(x)=0.92 将上面各数值代入(10—18)式得： . . 1 2  −  y y S = 0.4477 0.92 (1.52 1.28) 12 2 0.8742 2 =         −  + 于是 3.85 0.4477 10.3514 12.0758 = − − t = 查t值表，当自由度为43时(见表10—6误差自由度)，t0.01(43)=2.70(利用线性内插法计算)， |t| ＞t0.01(43)，P＜0.01，表明对照与食欲添加剂1号配方校正50日龄平均重间存在着极显著的差异，这里表现为1号配方的校正50日龄平均重极显著高于对照。其余的每两处理间的比较都须另行算出  −  . . i j y y S ,再进行t检验。（2）最小显著差数法利用t检验法进行多重比较，每一次比较都要算出各自的  −  . . i j y y S ，比较麻烦。当误差项自由度在 20以上，x变量的变异不甚大(即x变量各处理平均数间差异不显著)，为简便起见，可计算一个平均的  −  . . i j y y S 采用最小显著差数法进行多重比较。  −  . . i j y y S 的计算公式如下：  −  . . i j y y S =         − + ( 1) 1 2 ( ) ( ) ' SS k SS n MS e x e t x (10-19) 公式中SSt(x)为x变量的处理间平方和。然后按误差自由度查临界t值，计算出最小显著差数： (dfe) LSD t   =  −  . . i j y y S (10-20) 本例x变量处理平均数间差异极显著，不满足“x变量的变异不甚大”这一条件，不应采用此处所介绍的最小显著差数法进行多重比较。为了便于读者熟悉该方法，仍以本例的数据说明之。此时  −  . . i j y y S = 0.4354 0.92 (4 1) 0.83 1 12 2 0.8742 =        − +  由 ' e df =43,查临界t值得：t0.05(43)=2.017，t0.01(43）=2.70 于是 LSD0.05=2.017×0.4353=0.878 LSD0.01=2.70×0.4353=1.175

205 不同食欲添加剂配方与对照校正50日龄平均重比较结果见表10—8。表10—8 不同食欲添加剂配方与对照间的效果比较表食欲添加剂配方校正50日龄平均重对照校正50日龄平均重差数 1 12.0758 10.3514 1.7244** 2 12.1821 10.3514 1.8307** 3 12.3398 10.3514 1.9884** 多重比较结果表明：食欲添加剂配方1、2、3号与对照比较，其校正50 日龄平均重间均存在极显著的差异，这里表现为配方1、2、3号的校正50日龄平均重均极显著高于对照。 (3)最小显著极差法当误差自由度在20以上，x变量的变异不甚大，还可以计算出平均的平均数校正标准误 y S ，利用LSR 法进行多重比较。 y S 的计算公式如下： y S =         − + ( 1) 1 ( ) ( ) ' SS k SS n MS e x e t x (10-21) 然后由误差自由度 e df  和秩次距k查SSR表（或q表），计算最小显著极差： LSR = SSRS y （10-22）对于【例10.1】资料，由于不满足“x变量的变异不甚大”这一条件，不应采用此处所介绍的LSR法进行多重比较。为了便于读者熟悉该方法，仍以【例10.1】的数据说明之。此时 ' MSe =0.8742，n=12，SSt(x)=0.83，SSe(x)=0.92，k=4，代入(10—21)式可计算得： y S = 0.3078 0.92 (4 1) 0.83 1 12 0.8742 =        − + SSR值与LSR值见表10—9。表10—9 SSR值与LSR值表秩次距k 2 3 4 SSR0.05 2.86 3.01 3.10 SSR0.01 3.82 3.99 4.10 LSR0.05 0.883 0.929 0.957 LSR0.01 1.179 1.232 1.266 各处理校正50日龄平均重多重比较结果见表10—10。表10—10 各处理校正50日龄平均重多重比较表（SSR法）处理 ' i. y ' i. y -10.3514 ' i. y -12.0758 ' i. y -12.1821 配方3 12.3398 1.9884** 0.2640 0.1577 配方2 12.1821 1.8307** 0.1063 配方1 12.0758 1.7244** 对照 10.3514 多重比较结果表明：食欲添加剂配方3、2、1号的哺乳仔猪校正50 日龄平均重极显著高于对照，不同食欲添加剂配方间哺乳仔猪校正50日龄平均重差异不显著

生物统计学：《生物统计附试验设计》课程教材（第三版）第十章 协方差分析

生物统计学：《生物统计附试验设计》课程教材（第三版）第十章协方差分析