0协方差分析最新版

易读文库协方差 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 第一节协方差分析的意义协方差分析有二个意义，一是对试验进行统计控制，二是对协方差组分进行估计，现分述如下。一、对试验进行统计控制为了提高试验的精确性和准确性，对处理以外的一切条件都需要采取有效措施严加控制，使它们在各处理间尽量一致，这叫试验控制。但在有些情况下，即使作出很大努力也难以使试验控制达到预期目的。例如：研究几种配合饲料对猪的增重效果，希望试验仔猪的初始重相同，因为仔猪的初始重不同，将影响到猪的增重。经研究发现：增重与初始重之间存在线性回归关系。但是，在实际试验中很难满足试验仔猪初始重相同这一要求。这时可利用仔猪的初始重（记为X）与其增重（记为y）的回归关系，将仔猪增重都矫正为初始重相同时的增重，于是初始重不同对仔猪增重的影响就消除了。由于矫正后的增重是应用统计 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 将初始重控制一致而得到的，故叫统计控制。统计控制是试验控制的一种辅助手段。经过这种矫正，试验误差将减小，对试验处理效应估计更为准确。若y的变异主要由x的不同造成（处理没有显著效应），则各矫正后的y•间将没有显著差异（但原y间的差异可能是显著的）。若y的变异除掉x不同的影响外，尚存在不同处理的显著效应，则可期望各y•间将有显著差异（但原y间差异可能是不显著的）。此外，矫正后的y和原y的大小次序也常不一致。所以，处理平均数的回归矫正和矫正平均数的显著性检验，能够提高试验的准确性和精确性，从而更真实地反映试验实际。这种将回归分析与方差分析结合在一起，对试验数据进行分析的方法，叫做协方差分析（analysisofcovariance）。、估计协方差组分在第八章曾介绍过 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 示两个相关变量线性相关性质与程度的相关系数的计算公式:为（X-x）（y-V）:;二（x—x）2^（y—y）2若将公式右端的分子分母同除以自由度（n-1），得T)(10-1)'（x-x）（y-y）/（n-1）其中'（X—X）2、、、2n-1是x的均方MSx，它是x的方差二%的无偏估计量;'(y-y)2n_1是y的均方MSy，它是y的方差匚2的无偏估计量;、(X_x)(y_y)n-1称为x与y的平均的离均差的乘积和，简称均积，记为MPxy='(x_X)(y_y)n-1丁(瓦x)(瓦y)肖一nn—1MPxy,即卩(10-2)与均积相应的总体参数叫协方差(covariance)，记为COV(x,y)或二xy。统计学证明了,均积MPxy是总体协方差COV(x,y)的无偏估计量，即EMPxy=COV(x,y)。于是，样本相关系数r可用均方MSx、MSy，均积MPxy表示为：MPxy/、r—y(10-3).MSxMSy相应的总体相关系数示如下：”可用x与y的总体 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 差cx、二y，总体协方差COV(x,y)或;「xy表、_COV(x,y)_匚xy二x二y(10-4)均积与均方具有相似的形式，也有相似的性质。在方差分析中，一个变量的总平方和与自由度可按变异来源进行剖分，从而求得相应的均方。统计学已证明：两个变量的总乘积和与自由度也可按变异来源进行剖分而获得相应的均积。这种把两个变量的总乘积和与自由度按变异来源进行剖分并获得获得相应均积的方法亦称为协方差分析。在随机模型的方差分析中，根据均方MS和期望均方EMS的关系，可以得到不同变异来源的方差组分的估计值。同样，在随机模型的协方差分析中，根据均积皿味口期望均积EM的关系，可得到不同变异来源的协方差组分的估计值。有了这些估计值，就可进行相应的总体相关分析。这些分析在遗传、育种和生态、环保的研究上是很有用处的。由于篇幅限制，本章只介绍对试验进行统控制的协方差分析。第二节单因素试验资料的协方差分析设有k个处理、n次重复的双变量试验资料，每处理组内皆有n对观测值x、y,则该资料为具kn对x、y观测值的单向分组资料，其数据一般模式如表10—1所示。表10—1kn对观测值x、y的单向分组资料的一般形式处理处理1处理2处理i处理k观测指标xyXyxyXyxnynX21y21Xi1yi1Xk1yk1观测值X12y12X22y22Xi2yi2Xk2yk2xj、yij(i=1,2,X1jyuX2jy?XijyijXkjykjj=1,2,-n)X1ny1nX2ny2nXinyinXknykn总和X1.y1.X2・y2.Xi.yi.Xk.yk.平均数x1-Y1-X2-Y2-Xi.yi-Xk.Yk.表10—1的x和y变量的自由度和平方和的剖分参见单因素试验资料的方差分析方法一节。其乘积和的剖分则为：总变异的乘积和SPT是xij与x..和yij与y..的离均差乘积之和，即：knknSFT八v(Xj-X..)(yj-y.J八vXjyji4ji4j4x..y..kn(10-5)dft=kn-1kk其中，X..=7Xi.,y..=7yiVi4o处理间的乘积和SR是Xi.与X..和yi.与y..的离均差乘积之和乘以n,即:kSR=n》(XiiA1J…x..)(yi.…y..)；ni二Xi.yi.Xi.yi.kn(10-6)dft二k-1处理内的乘积和SR是xj与xi.和yij与yi.的离均差乘积之和，即：knkn,k__1SF^='、(xij-Xj.)(yj-yi.)=''xijyij-xi.yi.=SRT-SR(10-7)i土j土iajdniddfe=k(n-1)以上是各处理重复数n相等时的计算公式，若各处理重复数n不相等，分别为m、n2、…、m,k其和为7ni,则各项乘积和与自由度的计算公式为：i吕kmSR一寸寸为.yi——Xijyij_~ki1(10-8)kdfT八ni-1i=1x1.y1.x2.y2.sp+.n1n2dft-k-1kmor>ppl|X1.y1.SRXijyiji=1j_n1kdfe=SLni-k=dfT-dfti=1Xk.yk.x..y..一-knktn-niid十竺2+...+3J=SRt-srn2nk(10-9)有了上述SP和df,再加上x和y的相应SS,就可进行协方差分析。【例10.1】为了寻找一种较好的哺乳仔猪食欲增进剂，以增进食欲，提高断奶重，对哺乳仔猪做了以下试验：试验设对照、配方1、配方2、配方3共四个处理，重复12次，选择初始条件尽量相近的长白种母猪的哺乳仔猪48头，完全随机分为4组进行试验，结果见表10—2,试作分析。此例，X..X2.X3.X4.=18.25+15.40+15.65+13.85=63.15y..二y1.y2.y3.y4.=141.80+130.10+144.80+133.80=550.50k=4,n=12,kn=4X12=48表10—2不同食欲增进剂仔猪生长情况表（单位：kg）处理对照配方1配方2配方3观测指标初生重x50日龄重y初生重x50日龄重y初生重x50日龄重y初生重x50日龄重y1.5012.401.3510.201.1510.001.2012.401.8512.001.209.401.1010.601.009.801.3510.801.4512.201.1010.401.1511.601.4510.001.2010.301.059.201.1010.601.4011.001.4011.301.4013.001.009.20观察值1.4511.801.3011.401.4513.501.4513.90xij,yij1.5012.501.1512.801.3013.001.3512.801.5513.401.3010.901.7014.801.159.301.4011.201.3511.601.4012.301.109.601.5011.601.158.501.4513.201.2012.401.6012.601.3512.201.2512.001.0511.201.7012.501.209.301.3012.801.1011.00总和Xj.,yi.18.25141.8015.40130.8015.65144.8013.85133.80平均xi.,Vi.1.5211.821.2810.841.3012.071.151.15协方差分析的计算步骤如下：（一）求X变量的各项平方和与自由度1、总平方和及自由度222ccV、、•2x..222263.15263.152TOC\o"1-5"\h\zSSnx）-x2（1.5021.852-...-1.102）84.83251.75485（x）jkn48dfT(x)=kn-1=4X12-1=472、处理间平方和与自由度63.15二0.83dft(x)=k48t(x)k2Il2X・・I2222SS(x)=-瓦Xi2.———=一(18.252+15.402+15.652+13.852)nykn12-1=4-1=33、处理内平方和与自由度SQx）=SSr（x）-sst（x）=1-75-0-83=0.92dfe(x)=dfT(x)-dft(x)=47-3=44(二)求y变量各项平方和与自由度2550.522550.52=6410.3196.76484822550.502133.802)11.68481、总平方和与自由度22y222SSr(y)-yj1-=(12.4012.00-...11.00)kndfT(y)=kn-1=4X2-仁472、处理间平方和与自由度SSt(y)J'yi2.—y..1(141.802130.802144.802nkn12dft(y)=k-1=4-1=33、处理内平方和与自由度SSe(y)=SSr(y)-SS(y)=96.76-11.68=85.08dfe(y)=dfT(y)-dft(y)=47-3=44（三）求x和y两变量的各项离均差乘积和与自由度1、总乘积和与自由度SFT二'「Xjyj-字1j4kn=1.5012.401.8512.00...1.1011.00—63.15550.50412=732.50-63.15550.50=8.25412dfT(x,y)=kn-1=4M2-1=472、处理间乘积和与自由度—1lx..y..SRXiy.-nykn丄(18.25141.8015.40130.1015.65144.8013.85133.80)63.15550.5012412=1.64dft（x,y）=k-1=4-1=33、处理内乘积和与自由度SPe=SFT-SR=8.25-1.64=6.61dfe（x,y）=dfy（x,y）-dft（x,y）=47-3=44平方和、乘积和与自由度的计算结果列于表10—3。变异来源dfSSxSSySRy处理间（t）30.8311.681.64处理内（误差）（e）440.9285.086.61总变异（T）471.7596.768.25表10—3x与y的平方和与乘积和表（四）对x和y各作方差分析（表10—4）变异来源dfX变量y变量F值SSMSFSSMSF处理间30.830.2813.33**11.683.892.02F0.05=2.82处理内（误差）440.920.02185.081.93总变异471.7596.76F0.01=4.26表10—4初生重与50日龄重的方差分析表分析结果表明，4种处理的供试仔猪平均初生重间存在着极显著的差异，其50日龄平均重差异不显著。须进行协方差分析，以消除初生重不同对试验结果的影响，减小试验误差，揭示出可能被掩盖的处理间差异的显著性。（五）协方差分析1、误差项回归关系的分析误差项回归关系分析的意义是要从剔除处理间差异的影响的误差变异中找出50日龄重(y)与初生重(x)之间是否存在线性回归关系。计算出误差项的回归系数并对线性回归关系进行显著性检验，若显著则说明两者间存在回归关系。这时就可应用线性回归关系来校正y值(50日龄重)以消去仔猪初生重(x)不同对它的影响。然后根据校正后的y值(校正50日龄重)来进行方差分析。如线性回归关系不显著，则无需继续进行分析。回归分析的步骤如下：离回归平方和与相应的自由度(1)计算误差项回归系数，回归平方和，从误差项的平方和与乘积和求误差项回归系数byx(e)_SFeSSe(x)卷7・1848(10-10)误差项回归平方和与自由度22(10-11)SFe26.612SSR(e)=SSe(x)0.9247.49dfR(e)=1误差项离回归平方和与自由度ssr(e)=sse(y)-ss^=85.08-47.49=37.59(10-12)dfr(e)=dfe(y)-dfR(e)=44-1=43(2)检验回归关系的显著性(表10—5)表10—5哺乳仔猪50日龄重与初生重的回归关系显著性检验表变异来源SSdfMSFF0.01误差回归47.49147.4954.32**7.255误差离回归37.59430.8742误差总和85.0844F检验表明，误差项回归关系极显著，表明哺乳仔猪50日龄重与初生重间存在极显著的线性回归关系。因此，可以利用线性回归关系来校正y,并对校正后的y进行方差分析。2、对校正后的50日龄重作方差分析(1)求校正后的50日龄重的各项平方和及自由度利用线性回归关系对50日龄重作校正,并由校正后的50日龄重计算各项平方和是相当麻烦的，统计学已证明，误差平方和及自由度等于其相应变异项的离回归平方和及自由度，由度可直接由下述公式计算。①校正50日龄重的总平方和与自由度，即总离回归平方和与自由度'SFT2SSr—SSp(y)-SSR(y)—SSr(y)96.76SST(x)校正后的总平方和、因此，其各项平方和及自8.25257.851.75(10-13)dfT=dfT(y)-dfR(y)=47-1=46②校正50日龄重的误差项平方和与自由度，即误差离回归平方和与自由度26.61285.0837.59SSe(x)0.92SSb-SSe(y)_SS^(e)—SS^(y)SFe(10-14)dfe=dfe(y)-dfe(R)=44-仁43上述回归自由度均为1，因仅有一个自变量X。③校正50日龄重的处理间平方和与自由度SS=SSt_SSe=57.87-37.59=20.28(10-15)dft=dfr-dfe=k-1=4-1=3(2)列出协方差分析表，对校正后的50日龄重进行方差分析(表10—6)查F表：Fo.o1(3,43)=4.275(由线性内插法计算)，由于F=7.63>F0.01Q43),P<0.01，表明对于校正后的50日龄重不同食欲添加剂配方间存在极显著的差异。故须进一步检验不同处理间的差异显著性，即进行多重比较。表10—6表10-2资料的协方差分析表变异来源dfS3ssSRyb校正50日龄重的方差分析Fdf'SS'MS处理间(t)30.8311.681.64机误(e)440.9285.086.617.18484337.590.8742总和(T)471.7596.768.254657.87校正处理间320.286.767.63**3、根据线性回归关系计算各处理的校正50日龄平均重误差项的回归系数byx(e)表示初生重对50日龄重影响的性质和程度，且不包含处理间差异的影响，于是可用byx(e)根据平均初生重的不同来校正每一处理的50日龄平均重。校正50日龄平均重计算公式如下：yi.=yi.-byx(e)(xi・一'乂..)(10-16)公式中：yi.为第i处理校正50日龄平均重；yi.为第i处理实际50日龄平均重(见表10—2);xi.为第i处理实际平均初生重(见表10—2);x6315x..为全试验的平均数，x..二仝二^6315=1.3156kn48byx(e)为误差回归系数，byx(e)=7.1848将所需要的各数值代入(10—佝式中，即可计算出各处理的校正50日龄平均重(见表10—7)。表10—7各处理的校正50日龄平均重计算表处理乂).-X..byx(e)(xi.-x..)实际50日龄平均重校正50日龄平均重yi.-byx(e)(xi.-X..)对照1.52-1.3156=0.20447.1848>0.2044=1.468611.8211.82-1.1686=10.3514配方11.28-1.3156=-0.03567.1848X(-0.0356)=-0.258810.8410.84+0.2558=12.0758配方21.30-1.3156=-0.01567.1848X(-0.0156)=-0.112112.0712.07+0.1121=12.1821配方31.15-1.3156=-0.16567.1848X(-0.1656)=-1.189811.1511.15+1.1898=12.33984、各处理校正50日龄平均重间的多重比较各处理校正50日龄平均重间的多重比较，即各种食欲添加剂的效果比较。(1)t检验检验两个处理校正平均数间的差异显著性，可应用t检验法：—F—F引.—yj.t:-(10-17)Sy；%’2(x「一Xj.)IS「y=MS；-+j—(10-18)2j.¥nSSe(x)式中，yi.-y,为两个处理校正平均数间的差异；S・■为两个处理校正平均数差数标准误；y.亠.MSe为误差离回归均方；n为各处理的重复数；Xi.为处理i的x变量的平均数；Xj.为处理j的X变量的平均数；SS(x)为X变量的误差平方和例如，检验食欲添加剂配方1与对照校正50日龄平均重间的差异显著性:yi.—丫2.=10.3514-12.0758=-1.7244MSe=37.59/43=0.8742n=12=1.52,x2.=1.28,SSe(x)=0.92将上面各数值代入(10—18)式得:SZ「0.87422+(^2811=0.4477120.92于是10.3514-12.07580.4477二-3.85),|t|查t值表，当自由度为43时(见表10—6误差自由度)，t°.01(43)=2.70(利用线性内插法计算>t0.01(43),PV0.01,表明对照与食欲添加剂1号配方校正50日龄平均重间存在着极显著的差异，这里表现为1号配方的校正50日龄平均重极显著高于对照。其余的每两处理间的比较都须另行算出Sy_-,再进行t检验。yi.-¥j.(2)最小显著差数法利用t检验法进行多重比较，每一次比较都要算出各自的S_,yi.-yj.'比较麻烦。当误差项自由度在20以上，x变量的变异不甚大(即x变量各处理平均数间差异不显著)，为简便起见，可计算一个平均的S_，一，采用最小显著差数法进行多重比较。S_*_*yi.jj.yi.—yj.Syi.-yj.SSt(x)SSe(x)(k-1)的计算公式如下：(10-19)公式中SSt(x)为x变量的处理间平方和。然后按误差自由度查临界t值，计算出最小显著差数：LSW兀二，(10-20)本例x变量处理平均数间差异极显著，不满足“x变量的变异不甚大”这一条件，不应采用此处所介绍的最小显著差数法进行多重比较。为了便于读者熟悉该方法，仍以本例的数据说明之。此时SyiT2^0.8742120.8310.92(4-1)=0.4354由dfe=43,查临界t值得：t0.05(43)=2.017,t°.01(43)=2.70于是LSD°.05=2.017@4353=0.878LSD0.01=2.700.4353=1.175不同食欲添加剂配方与对照校正50日龄平均重比较结果见表10—8。表10—8不同食欲添加剂配方与对照间的效果比较表食欲添加剂配方校正50日龄平均重对照校正50日龄平均重差数112.075810.35141.7244**212.182110.35141.8307**312.339810.35141.9884**多重比较结果表明：食欲添加剂配方1、2、3号与对照比较，其校正50日龄平均重间均存在极显著的差异，这里表现为配方1、2、3号的校正50日龄平均重均极显著高于对照。（3）最小显著极差法当误差自由度在20以上，x变量的变异不甚大，还可以计算出平均的平均数校正标准误Sy，利用LSR法进行多重比较。Sy的计算公式如下：Sy=然后由误差自由度dfe和秩次距k查SSR表（或q表）LSR-.=SSR.SySSt(x)1SSe（x）（k-1）_，计算最小显著极差：(10-21)(10-22)时MSe=0.8742，n=12,Sy■0.8742-*0.83121|L0.92(4_1)"0.3078对于【例10.1】资料，由于不满足“x变量的变异不甚大”这一条件，不应采用此处所介绍的LSR法进行多重比较。为了便于读者熟悉该方法，仍以【例10.1】的数据说明之。此SS（x）=0.83，S£（x）=0.92，k=4，代入（10—21）式可计算得:SSR值与LSR值见表10—9。表10—9SSR值与LSR值表秩次距k234SSR0.052.863.013.10SSR0.013.823.994.10LSR0.050.8830.9290.957LSR0.011.1791.2321.266各处理校正50日龄平均重多重比较结果见表10—10。表10—10各处理校正50日龄平均重多重比较表（SSR法）处理y'-10.3514yi'-12.0758-12.1821配方312.33981.9884**0.26400.1577配方212.18211.8307**0.1063配方112.0758**1.7244对照10.3514多重比较结果表明：食欲添加剂配方3、2、1号的哺乳仔猪校正50日龄平均重极显著高于对照，不同食欲添加剂配方间哺乳仔猪校正50日龄平均重差异不显著。习 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 1、何为试验控制？如何对试验进行统计控制？2、什么是均积、协方差？均积与协方差有何关系？、对试验进行统计控制的协方差分析的步骤有哪些？、一饲养试验，设有两种中草药饲料添加剂和对照三处理，重复9次，共有27头猪参与试验，两个月增重资料如下。由于各个处理供试猪只初始体重差异较大，试对资料进行协方差分析。中草药饲料添加剂对猪增重试验结果表（单位：kg）处理2号添加剂1号添加剂对照组观测指标初重x增重y初重x增重y初重x增重y30.535.527.529.528.526.524.525.021.519.522.518.523.021.520.018.532.028.5观20.520.522.524.519.018.0测21.025.524.527.516.516.0值28.531.526.028.535.030.522.522.518.519.022.520.518.520.528.531.515.516.021.524.520.518.517.016.0（b=0.9832，线性回归关系极显著）。5、四种配合饲料的比较试验，每种饲料各有供试猪10头，供试猪的初始重（kg）及试验后的日增重（kg）列于下表，试对试验结果进行协方差分析。处理I号料H号料山号料W号料观测指标始重x增重y始重x增重y始重x增重y始重x增重y360.89280.64280.55320.52300.80270.81220.62270.58260.74270.73260.58250.64观230.80240.67220.58230.62测260.85250.77230.66270.54值300.68230.67200.55280.54200.73200.64220.60200.55190.68180.65230.71240.44200.80170.59180.55190.51160.58200.57170.48170.51（b=0.0073，线性回归关系极显著）