《协方差分析理论与案例》由会员分享，可在线阅读，更多相关《协方差分析理论与案例（12页珍藏版）》请在金锄头文库上搜索。
1、1协方差分析理论与案例假设我们有 N 个个体的 K 个属性在 T 个不同时期的样本观测值，用 ,ity，N,t=1，T,k=1，K 表示。一般假定 的观测值是某随机实验的itx y结果，该实验结果在属性向量 和参数向量 下的条件概率分布为 。使x(,)fyx用面板数据的最终目标之一就是利用获取的信息对参数 进行统计推断，譬如常假设假定的 是关于 的线性函数的简单模型。协方差分析检验是识别样本y波动源时广泛采用的方法。方差分析：常指一类特殊的线性假设，这类假设假定随机变量 y 的期望值仅与所考察个体所属的类（该类由一个或多个因素决定）有关，但不包括与回归有关的检验。而协方差分析模型具有混合特征，
2、既像回归模型一样包含真正的外生变量，同时又像通常的方差一样允许每个个体的真实关系依赖个体所属的类。常用来分析定量因素和定性因素影响的线性模型为： *,1,itititityxuNtT从两个方面对回归系数估计量进行检验：首先，回归斜率系数的同质性；其次，回归截距系数的同质性。检验过程主要有三步：(1) 检验各个个体在不同时期的斜率和截距是否都相等；(2) 检验（各个体或各时期的）回归斜率（向量）是否都相等；(3) 检验各回归截距是否都相等。显然，如果接受完全同同质性假设（1） ，则检验步骤中止。但如果拒绝了完全同质性性假设，则（2）将确定回归斜率是否相同。如果没有拒绝斜率系数的同质性假设，则（3
3、）确定回归截距是否相等。 （1）是从（2） 、 （3）分离出来的。基本思想：在作两组或多组均数 ， ， 的假设检验前，用线性回1y2ky归分析方法找出协变量 X与各组 Y之间的数量关系，求得在假定 X相等时修定均数 ， ， 然后用方差分析比较修正均数间的差别，这就是协方差1y2ky分析的基本思想。协方差分析的应用条件：要求各组资料都来自正态总体，且各组的方差相等；（ 检验或方差分析的条件）各组的总体回归系数 相等，且都不等t i于 0（回归方程检验） 。因此，应用协方差分析前，要对资料进行方差齐性检验和回归系数的假设检验（斜率同质性检验） ，只有满足上述两个条件之后才能应用，否则不宜使用。各比
4、较组协变量 X 与分析指标 Y 存在线性关系（按直线回归分析方法进行判断） 。2各比较组的总体回归系数 相等，即各直线平行(绘出回归直线，看是否i平行) 。 协方差分析适用的资料：完全随机设计、随机区组设计、拉丁方设计、析因设计等资料；协变量 X 可以仅有一个，称一元协方差分析；协变量也可以有多个，称多元协方差分析。 协方差计算公式：相关系数： 22()()xyr将公式右端的分子分母同除以自由度(n-1)，得： 22()/1()xynr其中：是 x 的均方 MSx，它是 x 的方差 的无偏估计量；2()1n2x是 y 的均方 MSy，它是 y 的方差 的无偏估计量；2() 2y称为 x 与 y
5、 的平均的离均差的乘积和，简称均积，记为 MPxy，()1xn即 ()()=11xy xyyynMPn与均积相应的总体参数叫协方差（covariance） ，记为 COV(x,y)或 。统计xy学证明了，均积 MPxy是总体协方差 COV(x,y)的无偏估计量，即 EMPxy= COV(x,y)。于是，样本相关系数 r 可用均方 MSx、MS y，均积 MPxy表示为：yxMPS相应的总体相关系数 可用 x 与 y 的总体标准差 、 ，总体协方差xyCOV(x,y)或 表示如下：xy(,)xyxyCOV3均积与均方具有相似的形式，也有相似的性质。在方差分析中，一个变量的总平方和与自由度可按变异
6、来源进行剖分，从而求得相应的均方。统计学已证明：两个变量的总乘积和与自由度也可按变异来源进行剖分而获得相应的均积。这种把两个变量的总乘积和与自由度按变异来源进行剖分并获得相应均积的方法亦称为协方差分析。1.协方差分析是将线性回归与方差分析相结合的一种分析方法；2.把对反应变量 Y 有影响的因素 X 看作协变量，建立 Y 对 X 的线性回归，利用回归关系把 X 值；3.化为相等，再进行各组 Y 的修正均数间比较。修正均数是假设各协变量取值固定在其总均数时的反应变量 Y 的均数。其实质是从 Y 的总离均差平方和 中，扣除协变量 X 对 Y 的回归2()平方和 ，对离回归平方和 作进一步分解后再进行
7、方差分析。2()方差分析的前提是除随机误差外，水平变量是影响观测值的唯一变量，方差分析数据结构： i ijijuteY协方差分析将方差分析与回归分析结合了起来，协方差分析数据结构： ()y e xi ij ijijutXuY第 i 组第 j 个观测值 一般均值第 i 组的组效应随机误差回归系数协变量效应4协方差案例：设有 k个处理、 n次重复的双变量试验资料，每处理组内皆有 n对观测值 x、 y，则该资料为具 kn对 x、 y观测值的单向分组资料，其数据一般模式如表101所示。表 1 kn对观测值x、y的单向分组资料的一般形式处 理 处理1 处理2 处理i 处理k观测指标 x y x y x
8、y x y观测值xij、y ij(i=1,2,kj=1,2,n)x11x12x1jx1ny11y12y1jy1nx21x22x2jx2ny21y22y2jy2nxi1xi2xijxinyi1yi2yijyinxk1xk2xkjxknyk1yk2ykjykn总 和 x1. y1. x2. y2. xi. yi. xk. yk.平均数 表1的x和y变量的自由度和平方和的剖分参见单因素试验资料的方差分析方法一节。其乘积和的剖分则为：总变异的乘积和 是 与 和 与 的离均差乘积之和，即：TSPijx.ijy.(1)knyxSPkinjijkinj ijijT .).)( 11 =kn-1 Tdf其中，
9、 。 knyknxyxkiikii .,.,.,. 11 处理间的乘积和 是 与 和 与 的离均差乘积之和乘以n，即：tSP.i.i(10-6)ki ikiiit kyxyxn11.).)(dft处理内的乘积和 是 与 和 与 的离均差乘积之和，即：eSPijx.iijy.i(10-7)kinj ki tTkiinjijiijije SPyxyx1 11.).)(=k(n-1) edf5以上是各处理重复数n相等时的计算公式，若各处理重复数n 不相等，分别为n1、n 2、n k，其和为 ，则各项乘积和与自由度的计算公式为：kin1kijkiiijTi yxSP11.= -1 (10-8)Tdfk
10、in1 kikt nyxyxSP121 .kdft- =SPT-SPtinjijeiyxSP1 knyxyxn.21= -k =dfT-dft (10-9)edfi1有了上述SP和df ，再加上x和y的相应SS ，就可进行协方差分析。【例10.1】 为了寻找一种较好的哺乳仔猪食欲增进剂，以增进食欲，提高断奶重，对哺乳仔猪做了以下试验：试验设对照、配方1、配方2、配方3共四个处理，重复12次，选择初始条件尽量相近的长白种母猪的哺乳仔猪48头，完全随机分为4组进行试验，结果见表102，试作分析。此例， =18.25+15.40+15.65+13.85=63.15. 4321xx=141.80+13
11、0.10+144.80+133.80=550.50yyk=4，n=12，kn=412=48表102 不同食欲增进剂仔猪生长情况表 （单位：kg）处 理 对照 配方1 配方2 配方3观 测指 标初生重x50日龄重 y初生重x50日龄重y初生重x50日龄重y初生重 x50日龄重y1.50 12.40 1.35 10.20 1.15 10.00 1.20 12.401.85 12.00 1.20 9.40 1.10 10.60 1.00 9.801.35 10.80 1.45 12.20 1.10 10.40 1.15 11.601.45 10.00 1.20 10.30 1.05 9.20 1.1
12、0 10.601.40 11.00 1.40 11.30 1.40 13.00 1.00 9.201.45 11.80 1.30 11.40 1.45 13.50 1.45 13.90观 察 值xij,yij1.50 12.50 1.15 12.80 1.30 13.00 1.35 12.8061.55 13.40 1.30 10.90 1.70 14.80 1.15 9.301.40 11.20 1.35 11.60 1.40 12.30 1.10 9.601.50 11.60 1.15 8.50 1.45 13.20 1.20 12.401.60 12.60 1.35 12.20 1.25
13、 12.00 1.05 11.201.70 12.50 1.20 9.30 1.30 12.80 1.10 11.00总 和 x i.,yi. 18.25 141.80 15.40 130.80 15.65 144.80 13.85 133.80平 均 ,1.52 11.82 1.28 10.84 1.30 12.07 1.15 1.15协方差分析的计算步骤如下：(一)求x变量的各项平方和与自由度1、总平方和及自由度 75.148.6325.841.63)0.85.10.(. 22222)(knxSijxT=kn-1=412-1=47)(df2、处理间平方和与自由度83.0415.6)8.13
14、65.40.125.8(.122222)( knxSkixt=k-1=4-1=3)(xtdf3、处理内平方和与自由度= - =1.75-0.83=0.92)(xeS)(xT)(xtS= - =47-3=44 )(df)()(tdf(二)求y变量各项平方和与自由度1、总平方和与自由度 76.9485.031.6485.0).10.24.1( 2222.)(knySijyT=kn-1=412-1=47)(df2、处理间平方和与自由度 68.1450.)8.130.48.130.4(12.1 22222.2)(knySiyt=k-1=4-1=3)(tdf3、处理内平方和与自由度= - =96.76-11.68=85.08)(yeS)(yT)(ytS= - =47-3=44)(df)()(tdf(三)求x和y 两变量的各项离均差乘积和与自由度1、总乘积和与自由度7knyxSPkinjijT.125.8140.5.630.72 12450.630.1. =kn-1=412-1=47),(yxTdf2、
