第四章回归与相关分析ppt课件

《第四章回归与相关分析ppt课件》由会员分享，可在线阅读，更多相关《第四章回归与相关分析ppt课件（88页珍藏版）》请在装配图网上搜索。

1、 实验设计与分析 第四章第四章 回归与相关回归与相关 主讲教师主讲教师 李得孝李得孝 4-1 直线回归与相关Linear Regression and Correlation变量之间的关系n统计相关关系 n1)因果关系回归分析 n2)平行关系相关分析354045500.20.40.60.811.201r 4-1-1直线回归ixiy7067555251525160641616.31610.91440.01440.71423.31471.31421.81547.11533.0表表4-1-3 播种至齐穗天数播种至齐穗天数x与总积温与总积温y的关系的关系图图4-1-4 y关于关于x的回归散点图的回归散

2、点图(Scatter plot)5055606570140014501500155016001650回归关系的模拟步骤：回归关系的模拟步骤：一、绘制散点图：一、绘制散点图：根据初等函数及其图像知识确定近根据初等函数及其图像知识确定近似的函数似的函数y=f(x)y=0+x0,1,2,iiiyxiny=0+x找到一条直线，即，确定0和的值，使偏向平方和最小。二二.参数估计参数估计利用最小二乘法估计0和：niiiniiexyQ12012)(niiiieniiiexbxbybQbxbybQ101000)(20)(2niiiebxbyQbb12000min)(,xbynybnxbllxxyyxxnxxn

3、yxyxbiiiixxxyiiiiiniiiiniiiiiii0221221)()()()(异同。不同教科书所用符号的的计算公式；，提问xyyyxxlll,:ixiy7067555251525160641616.31610.91440.01440.71423.31471.31421.81547.11533.0 522ix0000.58x4.13504iy4889.1500y 307202ix22.203122702iy4.787733iiyx0000.4449522307202xxl0689.4906894.1350422.203122702yyl2000.447894.135045224.7

4、87733xyl实例计算：实例计算：0860.100000.4442000.4478xxxyllb4988.9150 xbyb回归方程(Regression equation)为：xbxby0860.104988.9150Intercept回归截距Slope斜率确定直线回归方程的方法有最小二乘法、高斯牛顿法、麦夸特法、确定直线回归方程的方法有最小二乘法、高斯牛顿法、麦夸特法、牛顿法、梯度法、正割法等等。牛顿法、梯度法、正割法等等。提问：回归方程能否真是存在？三三.显著性检验显著性检验1)检验回归方程的显著性检验回归方程的显著性 222)()()(yyQyyUyyleyy平方和的分解：50556

5、06570140014501500155016001650y yy),(yx剩余平方和，残差平方和，离回归平方和残差回归平方和与剩余平方和的计算回归平方和与剩余平方和的计算eyyxxxyyyxxxyyyyyexxyyQUllllllUlxxyybUlyyyyyyyyyyyyyyQSSSPSSbxxbxbbbxbyyUyyl所以，/2)(2)(2)()()()()(/)()()(?)(2222222222200222/2nQSexy回归直线的显著性检验回归直线的显著性检验)2,1()2(nFnQUFe)()(),cov(yVxVyxyyxxxylllr 方法1:方差分析法方法2:相关系数法后讲2

6、5.12*0532.81)29(7823.39002866.45167)2(7823.39002866.451670689.490682866.451670000.444/2000.4478/01.022FnQUFUlQllUeyyexxxy798.0*9594.00689.490680000.4442000.447801.0rlllryyxxxy2)b的显著性检验的显著性检验1203.10000.444*)29(7823.3900/)2(/xxyxxebSSSlnQS499.3*0030.91203.10860.10)2()2(01.02tsbtntlnQblbsbtbxxexxb3)b0

7、3)b0 的假设检验的假设检验(过原点过原点998.2*99.1345.65499.91545.65)4445891(794.3900)1(2)7(01.00220tSbtSSxnnQSbxeb0:;0:000BHBHA在实践运用中，有时要检验回归直线过定点c,d)的问题，此时只需把察看点(x,y)变为(x-c,y-d)就成为过原点的问题了 四四.预测与控制预测与控制 当0bbxy显著时，可用于预测、控制等但必需留意，运用时x的取值范围只能在拟合回归方程时所用样本的范围内，不能外推。002020()12()112xxyxxxxQSnlnxxQSnlny0000,(2)(2)tnStnSyyyy

8、000,0(2)(2)yytnStnSyy 例3的99%预测图 45505560657075130014001500160017001800 xy0860.104988.915yStyy)7(01.0 yStyy)7(01.0 yStyy)7(01.0yStyy)7(01.04-2 相关分析及检验：相关分析及检验：)()(/),cov(yVxVyxyyxxxylllr 0:0H21/(2)2rrtr Srt nn)/(22aaatdftr这种检验不能对作出区间估计，亦不能进展H0:=c(c0)的检验！？处理途径：将r转换成Z1)相关系数等于非零常数的检验相关系数等于非零常数的检验.0,11ln

9、21;0,11ln21时当时当rrrrrrz那么z近似服从正态分布 31.0,11ln21;0,11ln21nzz时当时当),(2zzN)1,0(Nzuzzu检验:，【例4-2-1】在第一节【例4-1-1】中，12n*9394.0r，试检验与91.0的差别显著性，并给出95%的置信区间7329.100.32ln219394.019394.01ln2111ln21rrz5275.122.21ln2191.0191.01ln2111ln21z3333.0312131nz6163.03333.05275.17329.1u0796.13333.096.17329.113862.23333.096.17

10、329.127931.016642.816642.81e1e1e1e2.15922.159222111r9832.012144.11812144.1181e1e1e1e4.77254.772522222r2)两个相关系数的比较两个相关系数的比较【例4-2-2】在研讨贵阳水牛中，用39头牛的实测数据，得体重与胸围的相关系数为r1=0.8256*,，体重与体斜长的相关系数为r2=0.6678*,试比较两相关系数的的差别显著性 (参考成组数据的假设检验！1742.11744.08256.1ln218256.018256.01ln211z8068.03322.06678.1ln216678.01667

11、8.01ln212z2357.03391339131312121nnzz56.12357.08068.01742.12121zzzzu结论4-3 y关于x的非线性回归分析常用的一些简单非线性关系，都可经过变量代换变成直线回归分析，这种方法称为拟线性化法。此种方法获得的非线性方程普通不是最好的回归方程，只需用非线性最小二乘回归才干得到较为理想的结果【例4.1.4】用表4-1-4数据拟合Logistic曲线100.1Rxzae0 1020304050020406080100图图4-1-8Z与与x 的的Logistic关系曲线关系曲线 转换?zz100zzy)100(lnxz53.527.57143

12、03.316780040106.414.6520002.6827323931514.65.8493151.7663245722031.42.1847130.7814846422545.61.1929820.1764564373060.40.6556290.4221599873575.20.3297871.10.3075784090.20.1086472.2196470414595.40.0482183.0320222755097.50.0256413.663516646表表4-1-4 越冬棉红铃虫的化蛹进度越冬棉红铃虫的化蛹进度计算直线回归方程的SAS程序为：Data new;Input x

13、y;Cards;165 1.25 117180 1.85 160175 5.92 80120 9.42 259.98 28 11.01 40;Run;利用利用reg的程序的程序*/proc reg;model y=x;plot y*x /conf95;/*指定给出回归曲线图指定给出回归曲线图*/run;利用利用glm的程序：的程序：proc glm;model y=x;run;计算相关系数的SAS程序为：Data new;Input x y;Cards;70 1616.3 67 1610.9 55 1440.052 1440.7 51 1423.3 52 1471.351 1421.8 61

14、1547.1 64 1533.0;Proc corr;/用于计算相关系数的过程用于计算相关系数的过程Var x y;Run;4-3 对直线回归分析的进一步讨论n回归直线的失拟问题分析回归直线的失拟问题分析n回归直线间的比较回归直线间的比较n 1两条回归直线的比较两条回归直线的比较n 2k 条回归直线的比较条回归直线的比较4-3-1 回归直线的失拟问题分析回归直线的失拟问题分析 表表4-1-1 4-1-1 吸附率吸附率y y关于浓度关于浓度x x的直线回归的有反复的数据的直线回归的有反复的数据149.016.616.733.316.6516.67249.316.816.833.616.8016.

15、78349.516.816.933.716.8516.84449.816.917.033.916.9516.95550.017.017.134.117.0517.02650.217.017.134.117.0517.08iix1 iy2iy.iyiy0.3403xy在K个点上均作了m次反复察看，那么真正的察看误差平方和 122222111111111()kmkmkkrkeijiijiijiijijiijiQyyymyyym)1(1mkfe22111()()kmneijiiiiji mkQyyyy22211112()kmmnijiijiiiijji mkyymy yyy2221112()kmni

16、jiiiiiiji mkymymy yyy122211(2)()kneiiiiiiii mkQmymymyyyy12211()()kneiiiiii mkQmyyyy 12211()()knLfeeiiiiiiQQQmyyyy失拟平方和 12(1)Lfeefffnk m122222111111111()kmkmkkrkeijiijiijiijijiijiQyyymyyym失拟性检验：16,2,12,10,6(2 1)6,4eeLfkmnfff1262211112beijiijiQyy025.0185.342421.3424)1.346.333.33(2121.342422210.03160.0

17、250.0066LfeeQQQ110.0066 40.001650.3960.025 60.00417LfLfLfeeQfFQf10.05(,),(4,6)4.53LfLfeFF ffF故回归直线是不失拟的。进展失拟检验后，再进展回归方程的显著性检验 参见系统分组4-3-2 回归直线间的比较回归直线间的比较1.两条回归直线的比较两条回归直线的比较【例4-2-1】测定玉米地方种类“石榴子和“七叶白的叶片长宽乘积x，cm2和叶面积y，cm2的关系，试比较二者的异同经计算有,1070822;18,94.223,83.3191111xxlnyx;69447.0,82818.1,743652,51686

18、3*10111bbllxyyy11111.828180.69447;419.0,16.eeyxQf七叶白：,1351824;22;63.284,54.4082222xxlnyx;69719.0,20000.0;942483,658513*20222bbllxyyy22220.200000.69719;1423.3,20.eeyxQf 石榴子：1)平行性检验平行性检验 210:H21:AH12121212212(4)11bbxxxxbbbbtt nnSll122121244eeeQQQnnnn175.51363.1842422180.4193.142320095.0175.51135182411

19、070822121bbS2941.00095.000272.00095.069719.069447.02121bbSbbt030.2)36()4(05.02105.0tnnt0bbxy2)共截距检验共截距检验 02010:H0201:AH01020102010212222121212(4)11bbxxxxbbbbtt nnSxxnnll0467.4175.51135182454.408107082283.319221181220201bbS030.2)36(5019.00467.420000.082818.105.0tt3)在平行情况下的共截距检验在平行情况下的共截距检验(重合性检验重合性检验

20、 027.2,2073.00927.52826.03383.105.002010201tSbbtbb6960.024226461686135135182410708229424837436522121xxSSSSSPSPb3383.183.3196960.094.22301b2862.054.4086960.063.28402b11.33830.6960 xy 20.28620.6960 xy)(21221xxxxyyyyellbllQ平3212nnQe8518.48)13518241070822(6960.0658513516863371220927.58518.48135182410708

21、22)54.40883.319(221181)(112221221210201xxxxbbllxxnnS4在平行性被回绝下的共截距检验在平行性被回绝下的共截距检验 b0可解释为加权平均回归直线的截距(权重：ni/(n1+n2)2222212121222222111121212221110)()(xxxxnnxyxxxyxxyynnxbnxbnyb200()()ijiijijiijijeijijQybb xybb xy共222221112222222121120)()(xlnxlnbxyxxyxyxxxxijij2122221122221222221112)(0 xlnxlnxxxlnxlnSx

22、xxxxxxxb2.k条回归直线的比较条回归直线的比较,21NnnnkieeQQ 2iieeffNk 第一步：重合性分析0002010:kH，k21 假设它成立，那么把k个样本合并,计算出：bxby0 xxyyelblQ2重 2 Nfe重 eeQQ重)1(2)2()2(kkNNffee重)2,1()2/()1/()(kNkFkNQkQQFeee重重重合：合并回归直线不重合：平行性分析第二步：平行性分析kH210:i0不全相等 xxkxxxxxykxyxyllllllb2121 0iiibybx ki,2,1 2011()inkijiijeijQybbx平kinjiijiijixxbyy112)

23、(iixxiyyilbl2 由k条回归直线变为 k条平行回归直线所引起的剩余平方和变化为 eeQQ平(1)(2)1eefNkNkk平)2,1()2/()1/()(kNkFkNQkQQFeee平平第三步：共截距分析,:0002010kHi不全相等 1 1 12220021,2,kkkijijijjjiijjnb xn b xn b xbyNx ybxbikxijijiijexbbyQ20)(共iijxxiiijiijjijijijxlnbxyxy220222()/(1)(1,2)/(2)eeeeQQkFF kNkQNk共k条回归直线的比较条回归直线的比较:实例实例，【例4-2-2】为研讨某4个作

24、物种类从播种至出穗的总积温单位：d与天数x的关系，在同一管理条件下，各察看了7个点，if)8(eQ(1)fi(2)SSxi(3)SSyi(4)SPi(5)(4)/(2)bi(6)b0i（7）(1)-1 f ei（8）(3)-(4)2/(2)Q ei(9)品种1673897939.48364.011.33701.2353147.426394.950109.74291.410.87865.3353474.836928.0151457.711809.012.73755.9651185.646830.0185460.412325.514.85773.8552426.710234.5品种内（平行性）24

25、2890.9484967.236789.912.73618.86759.80756.23890.672316775。0（Qe（平）品种间3696.75486388.359225.9总变异（重合性）273587.7971355.646015.812.83750.1226381157。5（Qe(重）先检验重合性87.3)20,6(,8.12020/5.102346/)5.102345.381157(01.0*FF重87.3)20,6(,8.12020/5.102346/)5.102345.381157(01.0*FF重重合性被回绝，检验平行性10.3)20,3(,26.420/5.102343/)

26、5.102340.16775(05.0FF平重合性与平行性均被回绝，但是仔细斟酌，4条回归直线中的第4条与其他3条差别甚大，先把它剔除，仅比较前3条回归线的异同前3条回归直线的种类内平行性计算：*75.4415/9.78074/)9.78079.100970(重F24.115/9.78072/)9.78074.9096(平F检验阐明前3条回归直线是平行的，其方程为：1669.37 11.87xy2808.4711.87xy3814.44 11.87xy进一步研讨第4条与第3条的关系由于二者的回归截距很相近：228.2)10(,334.2908.073.1285.1405.0*tt平773.85

27、755.960.43740.93t4-4 多元线性回归与相关年 月 日4-4 多元回归与相关提纲n多元回归方程的建立n多元回归方程的检验与偏回归系数的检验n通径分析与偏相关分析x1x2xjxpy1x11x12x1jx1py12x21x22x2jx2py2.ixi1xi2xijxipyi.nxn1xn2xnjxnpyn和T.1T.2T.jT.pT.y表表4-4-1 多元回归数据多元回归数据iippiiixxxy2211001 122ppbb xb xb xy)(221100ppxbxbxbyb111222()()()pppyb xxb xxbxxyninippipiijexxbxxbyyQ112

28、1112)()(nijijpippiijexxxxbxxbyybQ1111)()()(2pj,2,1,0)()(2111 nipknijijkikkjijixxxxbxxyy4-4-1 多元回归方程的建立：pypppppyppypplblblbllbbbbbllbbblbl22112222212111212111111212122212ppxxpppplllllllllL12yyxypylllLT12(,)pb bbbxxxyL bL1xxxybL Lppxbxbxbyb221101()xxkjp pcL2T1200(,),1,2,1,jjjjxxbNcjpbNnx L x1,)(1,)(12

29、21121222pnfULQpflblblbyyUnfnTyyyLeyyeniUpypyyiniyyiyy)1(22RLQLURyyyy复相关系数复相关系数R (采用大写以与简单相关系数采用大写以与简单相关系数r相区别相区别 多元相关或称复相关多元相关或称复相关Multiple Correlation:是在是在m1个个变量中，变量中，m个变量的综合和一个变量的相关。个变量的综合和一个变量的相关。4-4-2 多元线性回归的统计检验 1.多元线性回归方程的显著性检验0:210pH:AHj不完全为0)1,()1(pnpFpnQpUFe)1,()1()1(22pnpFpnRpRF计算复相关系数后，直接

30、查计算复相关系数后，直接查Ra表表df=n-p-1,变量个数变量个数 p+1)1(pnFpFRaaa2 偏回归系数偏回归系数 bj 的检验的检验 0:0jH0:jAH)1(1pntpnQcbtejjjj)1,1(12pnFpnQcbFejjjj或 1pnQcSejjbj例4-4-1】支崇远对贵阳成年水牛39头实测了体重(y/kg)、胸围(x1/cm)、体斜长(x2/cm)和体高(x3/cm)，资料列于表4-4-2中，试建立三元线性回归方程 编号12345678体重443.5507.5462.5514.0471.5545540.5536胸围194200194211205204315207体长14

31、6150150153153153154142体高122.1123.5126.5134.5129.5125.5133.0128.5编号910111213141516体重468550.5492.0583442.5439.5477.5450胸围201199200210194190203194体长153160149160140147148135体高128.7127.5123.5138.6124121129118编号1718192021222324体重466480422413.5471.0414.5410428.5胸围190190185183193188179193体长135138.5140130145

32、133140140体高122124119.5114123.5119119116编号2526272829303132体重468517.5578620481702420491胸围190195.5207.5211203220197194体长155150160150137165142149体高120.5129.5128.5132.5130142.2124122编号33343536373839体重515483505465460404496胸围198200197192185187194体长150135153144154151152体高131.5128124119.5119.5123.0120.0表表4-4

33、-2 39头成年水牛实测体重、体尺数据表头成年水牛实测体重、体尺数据表25.847574,5.5741;5.1520646,7693222211xxxx233124885.1,613312.86;1133959xxx x1323965451.85,720207.7;19034.5x xx xy2129434166.25;3772275.25,2814701yx yx y7.23954923yx由表4-4-2数据计算一级统计数据：由一级统计数据计算二级统计数据：06.488,26.125,22.147,26.197321yxxx9359.315239)(212111xxl3974.261539)

34、(222222xxl5344.138439)(232333xxl212121123205.160139)(lxxxxl313131138103.181339)(lxxxxl323232232987.113839)(lxxxxl0897.14411039)(22yylyy1090.1759839)(111yxyxly4551.1292039)(222yxyxly9526.1119239)(333yxyxly平方和lxx交叉积lxy3152.93591601.32051813.81031601.32052615.39741138.29871813.81031138.29871384.5344xxL

35、17598.109012920.455111192.9526xyL计算 并估计各偏回归系数得出回归方程：1xxL10.00129958.51792230.00163248.51792230.00060100.00038250.0016340.00038250.00317522xx LT1T123(,)(3.49698,1.984559,1.871436)xxxyb b bL L3371.7283322110 xbxbxbyb123728.3371 3.496981.9845591.871436xxx y计算U,Q,U1,U2,U3并检验 5351.108128332211yyylblblbU5

36、546.35981ULQyye4418.94100012995.049698.3211211cbU4201.65530006010.0984559.1222222cbU0016.110300317522.0871436.1233233cbUUUUUUj8635.17066321表表4-4-3 4-4-3 三元线性回归方程的方差分析三元线性回归方程的方差分析12.4)35,1(05.0F145.7)35,1(01.0F40.4)35,3(01.0F变异来源 df（自由度）SS（平方和）MS（均方）FFaX119410.4418(U1)9410.44189.154*X216553.4201(U2)

37、6553.42016.3746*X311103.0016(U3)1103.00161回归3（p）108128.5546(U)36042.84535.06*剩余35（n-p-1）35981.5546(Q)1028.0444总变异38（n-1）144110.0897(Lyy)275%,yyRU l*0.8662R 01109.4871436.100317522.00016324.049698.33333111bccbb10502.2871436.100317522.00003825.0984559.13333222bccbb5335.1070250016.11035351.1081283*UUU5

38、562.370840016.11035546.359813*UQQ,36,2*eUff*9478.51362QUF25.5)36,2(01.0F0687.6132*21*1*0 xbxbyb12613.06874.011092.10502xx y%3.740897.1441105335.1070252R剔除x3后:111222()()()pppyb xxbxxbxxyT()ybxxT120000T120000001,()()1,1()()xxxxNynyNynyyxxLxxxxLxx002T1002T1001()()(1)11()()(1)xxeyxxeSQnpnSQnpnyxxLxxxxLx

39、x其中 10.00129958.51792230.00163248.51792230.00060100.00038250.0016340.00038250.00317522xx L200-7693/39=2.7436150-5741.5/39=2.7821123.5-4885.1/39=-1.7590(2.7436,2.7821,1.7590)例如:4-4-3 用线性回归处理多项式非线性回归问题用线性回归处理多项式非线性回归问题 ppxxxy2210,221ppxxxxxxppxxxy2211031622521421322110zzzzzzzy,2142132211zzxzxzxzx33221

40、10 xxxy如:如:令:令:表表4-4-4 4-4-4 废品率废品率y y与化学成分与化学成分x x的记录的记录xx 122xx xx 122xx y(%)(0.01%)y(%)(0.01%)1.303411560.444016001.003612960.564116810.733713690.304217640.903814440.424318490.813915210.354318490.703915210.404520250.603915210.414722090.504016000.60482304354045500.20.40.60.811.2图4-4-1 废品率y与化学成分x的散

41、点图及回归曲线图,02.10y 6263.01602.10y,6511x 6875.40166511x,267092x 3125.166916267092x44.221)(1212111xnxl1513685)(1222222xnxl18283)(121212112xxnxxll649.11)(1111yxnyxly05.923)(1222yxnyxly4982.1)(122ynylyy正那么方程组为05.923151368518283649.111828344.2212121bbbb,009301.0,8205.021bb484.1822110 xbxbyb218.4840.82050.00

42、9301xxy2102bbx，20210244b bbby极值:11.44009301.02)8205.0(0 x204 18.484 0.009301(0.8205)0.394 0.009301 y当某化学成分含量在0.44%左右时，平均废品率最小，约为0.39%4-4-4 通径分析与偏相关分析通径分析与偏相关分析 偏相关Partial Correlation：是在其他M2个变量皆固定时，指定的两个变量间的相关。偏回归系数bi不能反映自变量的相对重要性，由于：bi 带有详细的单位，单位不同无法比较；即使单位一样，Xi的变异程度不同，也不能比较。但可以采用规范化的偏回归系数，也称通径系数Pat

43、h coefficient：yxiixyiiSSSSbnSSnSSbP)1/(/1)1/(/1即对分子和分母分别除以即对分子和分母分别除以Y和和Xi的规范差，就可以消除单位和变异度的影响，的规范差，就可以消除单位和变异度的影响，其统计学意义是假设添加一个规范差单位，其统计学意义是假设添加一个规范差单位，Y将添加或减少将添加或减少Pi个规范差单位。个规范差单位。通径分析本质上是规范化的多元线性回归分析*11221212PPPyPyxxxxxxbbbSSSSypjSSbbnlSnlSyjjjjjjyyy,2,1,1,1*1b*2b*3b*eb12r23r1x2x3x13ry图图4-5-1 y关于各

44、关于各x的通径图的通径图 yyyrbbrbrrbrbbrrbrbrb3*3*232*1312*323*2*1211*1313*1212*1*T12(,),pb bbb)1,()1/()1(/)1/(/22pnpFpnRpRpnQpUF)1(1)1(2*pntpnRcbtjjjj2*T1 12 2yyp pyxyRUb rb rb r b R*T*2*112ppxxjjjkkjjj kbb r bb R b2*2(=2()jjjjkkjjjkkjkRbxyRRb r bxxy对 的直接决定系数）,相关路对 的相关决定系数jkjjyjjkjjbrbRRR2*2)(2决策系数【例4-5-2】关于小麦

45、产量y与其构成要素x1百粒重、x2(每株穗数、x3(每穗粒数x4(每穗粒重的通径分析实验为随机区组设计，参试种类10个，反复3次误差自在度为18，遗传型的自在度为9，遗传相关系数的显著性临界值为 x1X2X3X4yX110.2740.7060.5250.050X210.3000.4740.477X30.6650.256X410.440735.0)9(,602.0)9(01.005.0rr050.0525.0706.0274.0*4*3*2*1bbbb477.0474.0300.0274.0*4*3*2*1bbbb256.0665.0300.0706.0*4*3*2*1bbbb440.0665.

46、0474.0525.0*4*3*2*1bbbb313.0321.0*2*1bb908.0180.1*4*3bb14 42.02240.09101.29640.15650.09101.29460.09350.6280()1.29640.09352.62491.10930.15650.62801.10932.1175xxjicR861.0)9(,931.0,867.001.0*2RRrbRjyj325.3)9(,262.2)9(,)9(9)1(01.005.02*tttRcbtjjjj*4*3*2*12659.10,9826.11,5259.4,7136.3tttt%4.82,%2.139,%8.

47、9,%3.1024232221RRRR%2.22,%6.30,%5.53,%5.523141312RRRR%5.142,%9.263424RR通径直接作用间接作用X 对y的总作用决策系数x1对y0.3210.0867.10.8330.4770.050 x2对y0.3130.08820.00.3540.4300.477x3对y1.1800.22779.00.0940.06040.256x4对y0.9080.1690.1482.60.7850.440对y0.3650.36513.3yxx21yxx31yxx41yxx12yxx32yxx41yxx13yxx23yxx43yxx14yxx24yxx3

48、4偏相关分析偏相关分析 112112122(0)112121111pypypppyyyyprrrrrrrrrrrrR1111211111211212222212222112121212pypypypypppppypppppyyyypyyyyypyyllllccccllllccccllllccccllllccccAyyypyypyppppypypcccccccccccccccc2121222221111211yyjjjyjycccrjjiiijijcccr观 察 值x126.731.330.433.934.633.830.427.0 x273.459.065.958.264.664.662.17

49、1.4y10089591051102210971103992945观 察 值和平均X133.330.431.533.134.0410.431.5692x264.564.161.156.059.8824.763.4385y1074102910049951045133241024.9231【例4-5-3】测定13块中籼南京11号高产田的每亩穗数(x1),每穗实粒数(x2)和每亩稻谷产量(y).结果列于表4-5-4，试求偏相关系数 11121212221279.6077110.1046934.7692110.1046295.910838.9385943.769238.938528244.9231yy

50、yyyylllllllll A10.1633250.0615000.0055420.0615000.0265380.0020920.0055420.0020920.000223A*221112127174.09108.2956077.791046.110lllr 1*111943.76920.629479.6077 28244.9231yyyylrl l0135.09231.282449108.29539385.382222yyyylllr*221112129341.0026538.0163325.061500.0cccr*11119183.0000223.0163325.0005542.0y

51、yyycccr*22228600.0000223.0026538.0002092.0yyyycccr简单相关系数偏相关系数4-5 协方差分析年 月 日回想实验条件的控制n为了提高实验的准确度和准确度，必需使一切实验单元或区组内的实验单元的实验条件一致，这叫实验控制 n假设实验目的为y，而影响实验单元不一致的要素X可以量测，那么X叫做方差分析中的协变量，可进展统计控制例4-3-1】为比较三种肥料A1,A2,A3对苹果的增产效果，选了24株同龄苹果树，随机分成三组，并在第一年记下各树的产量x,kg，第二年每组随机施用一种肥料，并记录各树的产量 y,kg，结果见表4-3-1，试进展方差分析 表表4-

52、3-1 4-3-1 施用三种肥料的苹果产量施用三种肥料的苹果产量kg/kg/株株肥料观察值(xijj,yij)和平均A1x1j475853464956544440750.875y1j546663515666615046758.375A2x2j525364585961636647659.500y2j545367626263646949461.750A3x3j444846505957585341551.875y3j525854617064696649461.750129854.083145560.6254-5-1 单要素完全随机设计实验的协方差分析单要素完全随机设计实验的协方差分析表表4-3-2

53、4-3-2 不思索协变量不思索协变量X X的方差分析的方差分析1565.39375.30变异来源DFSSMSF肥料间260.75030.375误 差21830.87539.565总变异23891.625变异来源DFSSMSF肥料间2356.0833178.04176.34(Pr=0.0070)误 差21589.7528.0833总变异23945.8333表表 苹果树的根底消费力苹果树的根底消费力x x的方差分析的方差分析*19819.0r*29777.0r*39581.0r707.0)6(05.0r835.0)6(01.0r表表4-3-3 三种肥料供试树中三种肥料供试树中y关于关于x的直线回归

54、分析的直线回归分析iiiixbyb0ief)1/()3()2()7(2ieQeQ)(平eQ)(重eQ变异来源(1)SSxi(2)SSyi(3)SPi(4)=(3)/(1)bi(5)(6)(8)A1180.875293.875226.3751.25165.2979610.5543A2178.000223.500195.0001.09553.432669.8764A3230.875313.500257.7501.11643.8365625.746646.1773肥料内0.2100（平行589.750830.875675.1251.15156.76712048.8305性）2.0135肥料间356.

55、08360.75086.6252总变异945.833891.625765.7500.809616.839222271.6709（重合性）三条回归直线的重合性检验：*64.4520/8305.482/)8305.486709.271(重F5.85(2,20)0.01F三条回归直线的平行性检验：(48.830545.1773)/20.517146.1773/18F平55.3(2,18)0.05F10.21001.151558.375 1.1515(50.875)xx x y26.7671 1.151561.750 1.1515(59.500)xx x y32.0135 1.151561.750 1

56、.1515(51.875)xx xy阐明树的根底消费力x每添加1kg，那么y平均添加1.1515kg，其规范差为0.0643kg 0643.0750.589208305.48xxeeblfQS平平()ijiijiyb xxy假设把根底消费力影响校正到同一平均根底消费力程度)(xxbij.)(xxbyijij才是Ai的真正察看值()(.)1.1515(54.083)ij xijijijijyyb xxyx回归校正值 表表4-3-4 4-3-4 的校正值的校正值 ijy)(xijy)(xijy.iT.T肥料A162.156 61.490 64.247 60.308 61.853 63.793 61

57、.09661.611496.554A256.399 54.247 55.581 57.490 56.338 55.035 53.732 55.278444.100A363.61165.005 63.308 65.702 64.338 60.641 64.490 67.247514.3421455.00038.8820924145522.arTC75.381247.67490.61156.62222CSST92.3328/)342.514100.444554.496(222CSSA83.48ATeSSSSSS201211)1(rafe平*18.6820/83.482/92.332AF85.5)2

58、0,2(01.0F069.62)083.54875.50(1515.1375.58.)(1515.111)(1xxyyx512.55)083.54500.59(1515.1750.61.)(1515.122)(2xxyyx293.64)083.54875.51(1515.1750.61.)(1515.133)(3xxyyx224.2)(1)(3xxyy557.6)(2)(1xxyy224.20103bb557.60201bb0135.2,7671.6,2100.0030201bbb*64.4520/8305.482/)8305.486709.271(20/2/)(平平重重eeeQQQF三种肥料

59、均值间的比较变成三条平行的回归直线截距间的比较，有协变量的单要素完全随机等反复实验方差分析有以下几点：具有协变量的方差分析过程是a条回归直线比较的过程，它既等价于校正平均数间的显著性检验，又等价于三条平行回归直线截距间的显著性检验，因此这个方差分析过程又称为离回归分析 rjaixyijxijiyij,2,1;,2,1)(ijxijyiijxy)(模型充分阐明了这里的协方差分析就是 关于X的直线回归分析 iijy表表4-3-5 单要素完全随机实验的协方差分析离回归分析单要素完全随机实验的协方差分析离回归分析iefieQieMS)(31eeQQ 变异来源df回归计算离回归分析 SSxSSySPbi

60、F(1)处理+误差(总回归)23945.833891.625765.7500.809622271.6709(2)处理间2356.08360.75086.625(3)误差(处理内)21589.750830.875675.1251.15152048.83052.4415(4)重合性，即校正平均数间或回归截距间2222.8404111.420245.64*00212()11ikeikbbikxxxxxxaeQxxSnnlllf平平958.075.589)5.59875.50(824415.220201bbS784.075.589)875.50875.51(824415.220103bbS923.07

61、5.589)5.59875.51(824415.220203bbS)1)1(0000ratSbbtkibbki*128446.6958.0)7671.6(21.0t*318361.2784.0)21.0(0135.2t*325131.9923.0)7671.6(0135.2t086.2)20(05.0t845.2)20(01.0t4-5-2 单要素随机区组实验的协方单要素随机区组实验的协方差分析差分析表表4-3-6 “南优南优3号在不同施肥条件下的颖花数号在不同施肥条件下的颖花数x和结实率和结实率y)(xiT)(yiT ix iy)(xiy)()(,yjxjTTjjyx,)(xjy 区组处理I

62、II处理总和处理平均校正平均数 xyxy 14.59584.32618.911194.45559.564.76 24.09654.11628.201274.10063.566.03 33.94644.11648.051284.02564.065.95 43.90663.57697.471353.73567.567.22 53.45713.79677.241383.62069.067.84 63.48713.38726.861433.43071.568.87 73.39713.03746.421453.21072.568.18 83.14723.24696.381413.19070.566.02

63、 93.34693.04696.381383.19069.064.53 104.12614.76548.881154.44057.562.64 114.12634.75568.871194.43559.564.60 123.84673.60627.441293.72064.564.10 133.96644.50608.461244.23062.065.53 143.03753.01716.041463.02073.067.22 52.3993753.21910105.601847 3.74266.933.80165.003.771465.9643 66.7065.23表表4-3-7 4-3-7

64、 例例2 2的平方和与偏向积和的平方和与偏向积和)(xxXlSS)(yyylSS)(xylSP变异来源df区 组 间10.024026.030.79处 理 间136.8732694.4666.37误 差130.837282.476.44总 变 异277.7344802.9673.60表表4-3-8 4-3-8 例例2 2中中X X与与y y的方差分析的方差分析变异来源df颖花数（X）结实率（y）SSMSFSSMSF区组间10.02400.0240126.0326.034.10处理间136.87320.52888.20*694.4653.428.43*误差130.83720.064582.476

65、.34表表4-3-9 例例2中有关处置的协方差分析离回归分析中有关处置的协方差分析离回归分析 iefieQieMS582.1744.2342.4)(31eeQQ 变异来源df回归计算离回归分析 SSxSSySPbiF(1)处理+误差(总回归)267.7104776.9372.819.4432589.3786(2)处理间136.8732694.4666.37(3)误差(处理内)130.837282.476.447.6921232.93152.744(4)重合性，即校正平均数间或回归截距间1356.44714.342*6.44/0.8327 82.470.7750r 校正结实率没有显著性差别，表4

66、-3-8的结论是错的)()(3002)(2xkxikiyyxxkiebbSlxxrMSS处理内321(1)exxDxxMSlSral处理间处理内)(3eDftSLSD表表4-3-9 例例2中有关处置的协方差分析离回归分析中有关处置的协方差分析离回归分析 iefieQieMS4187.57443.28708.14)(31eeQQ 变异来源df回归计算离回归分析 SSxSSySPbiF(1)区组+误差(总回归)140.8612108.50-7.23-8.39531347.8023(2)区组间10.024026.03-0.79(3)误差(区组内)130.837282.47-6.44-7.69231232.93152.7443(4)重合性，即校正平均数间或回归截距间114.870814.870875.4)12,1(05.0F33.9)12,1(01.0F