生物统计学课件方差分析㈡.ppt





《生物统计学课件方差分析㈡.ppt》由会员分享,可在线阅读,更多相关《生物统计学课件方差分析㈡.ppt(20页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第六章第六章 方差分析(二)方差分析(二)(针对复因素试验结果进行的分析)(针对复因素试验结果进行的分析)n第一节第一节 正交试验数据正交试验数据n第二节第二节 完全随机试验数据完全随机试验数据n第三节第三节 系统分组试验数据系统分组试验数据n第四节第四节 裂区试验数据裂区试验数据*第六章要点提示第六章要点提示 本章讲授复因素试验数据的方差分析,是本课程的难点,它与正交表本章讲授复因素试验数据的方差分析,是本课程的难点,它与正交表的联系非同小可。学习时的联系非同小可。学习时要求掌握正交设计和系统分组设计的试验数据要求掌握正交设计和系统分组设计的试验数据的平方和、自由度分解原理和随机单位组设计(
2、交叉分组)的联系和区别;的平方和、自由度分解原理和随机单位组设计(交叉分组)的联系和区别;学会在有表可套的情况下利用正交表实现处理平方和学会在有表可套的情况下利用正交表实现处理平方和SSt、自由度、自由度DFt再再分解结果的分解结果的“自动检校自动检校”;牢记各因素不同水平之间以及各处理之间的平牢记各因素不同水平之间以及各处理之间的平均数进行多重比较时标准误均数进行多重比较时标准误SE的计算方法。的计算方法。涉及教材内容:涉及教材内容:第六章第三节,第六章第三节,第十二章第八节。第十二章第八节。作业布置:作业布置:教材教材P152-153 P152-153 T17、T18、T19;P269 P
3、269 T9。第一节第一节 正交试验数据正交试验数据6.1 在进行矿物质对架子猪的补饲试验中,考察补饲配方在进行矿物质对架子猪的补饲试验中,考察补饲配方(A)、用量、用量(B)、食盐、食盐(C)三个因素,每个因素各有三个因素,每个因素各有3个水平(个水平(a=3,b=3,c=3);按;按L 9(3 4)进进行正交试验,将行正交试验,将A、B、C依次安排到该正交表的第依次安排到该正交表的第、列构造正交组列构造正交组合合(k=9);试验实施时试验实施时,重复重复2次次(n=2),随机单位组设计,随机单位组设计。所得试验结果所得试验结果套用该正交表整理后列表如下,试予分析。套用该正交表整理后列表如下
4、,试予分析。一、数据整理一、数据整理 本例有因素水平及试验方案安排表,见教材本例有因素水平及试验方案安排表,见教材P251、P252(数据单位:(数据单位:kg)C=T 2/nk=1347.4 2/18=100860.3756 dfT=nk 1=2 9 1=17 SST=Y 2 C=63.4 2+68.9 2+92.82 100860.3756=1978.5444 SST=区组区组 SS r 处理处理 SS t 剩余剩余 SS 分解层次分解层次 SS SS SS SS 再分解层次再分解层次 dfT=区组区组 df r 处理处理 df t 剩余剩余 df 分解层次分解层次 df df df df
5、 再分解层次再分解层次 第一节第一节 正交试验数据正交试验数据ABC空空Tt1111163.467.4130.82122268.987.2156.13133364.966.3131.24212364.386.3150.65223170.288.5158.76231265.866.6132.47313271.489.0160.48321369.591.2160.79332173.792.8166.5T1 T2 T3418.1441.7487.6441.8475.5430.1423.9473.2450.3456.0448.9442.5Tr:612.1 735.3 T1347.4SS416.3344
6、185.2077202.881115.2011(校对校对=819.6243)819.6244df2222(校对校对=8)8第一节第一节 正交试验数据正交试验数据二、二、SST的分解及的分解及SS t的再分解的再分解 SS r Tr 2/k C 843.2355(612.1 2+735.3 2)/9 100860.3756 SS t Tt 2/n C 819.6244(130.8 2+156.12 +166.5 2)/2 100860.3756 SS T2/(nk/a)C 416.3344 SS T2/(nk/b)C 185.2077 SS T2/(nk/c)C 202.8811 SS T2/(
7、nk/3)C 15.2011 剩余剩余 SS SST SS r SS t 315.6845 用空列的用空列的SS实现实现“自动检校自动检校”:SS 15.2011 SS t SS SS SS 15.2022三、列三、列ANOVA表,进行表,进行F-testSOV DF SS MS F F 0.05区组区组 1 843.2355 处理处理 8 819.6244 102.45 3.10ns 3.97 2 416.3344 208.17 6.29*4.10 2 185.2077 92.60 2.80ns 4.10 2 202.8811 101.44 3.07ns 4.10 2 15.2011 模型误
8、差模型误差剩余剩余 8 315.6845 33.09 试验误差试验误差 总总 17 843.2355注:注:F 空列空列S 2/剩余剩余S 2 0.1926 (15.20112)/(315.68438)F 0.05,2,8 4.46 于是再分解得到的于是再分解得到的空列空列SS和分解层次的剩和分解层次的剩余余SS一样,都要视为误差一样,都要视为误差SSe而予以合并。而予以合并。Se 2 (15.2011315.6843)/(28)33.09第一节第一节 正交试验数据正交试验数据四四、多重比较多重比较、A因素三个水平间的比较因素三个水平间的比较 SE=MSe/6=33.096 =2.35按按df
9、e =10查查SSR临界值后比较如下:临界值后比较如下:k LSR 顺序顺序 A 0.05 0.01 0.05 0.01 A3 81.26 a A2 7.40 10.53 A2 73.62 b AB3 7.76 11.12 A1 69.68 b B、F-test已知,已知,B因素三因素三水平无显著差异水平无显著差异;(SE=MSe/6=33.096 =2.35)、F-test已知,已知,C因素三因素三水平差异不显著水平差异不显著;(SE=MSe/6=33.096 =2.35)、F-test已知,各正交组合间已知,各正交组合间差异不显著差异不显著。(SE=MSe/2=33.092 =4.07)三
10、、列三、列ANOVA表,进行表,进行F-testSOV DF SS MS F F 0.05区组区组 1 843.2355 处理处理 8 819.6244 102.45 3.10ns 3.97 2 416.3344 208.17 6.29*4.10 2 185.2077 92.60 2.80ns 4.10 2 202.8811 101.44 3.07ns 4.10 2 15.2011 模型误差模型误差剩余剩余 8 315.6845 33.09 试验误差试验误差 总总 17 843.2355注:注:F 空列空列S 2/剩余剩余S 2 0.1926 (15.20112)/(315.68438)F 0
11、.05,2,8 4.46 于是再分解得到的于是再分解得到的空列空列SS和分解层次的剩和分解层次的剩余余SS一样,都要视为误差一样,都要视为误差SSe而予以合并。而予以合并。Se 2 (15.2011315.6843)/(28)33.09第二节第二节 完全随机试验数据完全随机试验数据 例例6.2 为了从三种不同原料为了从三种不同原料(a=3)和三种不同温度和三种不同温度(b=3)中选择使酒精中选择使酒精产量最高的水平组合,设计了两因产量最高的水平组合,设计了两因素完全随机试验,每一水平组合重素完全随机试验,每一水平组合重复复4次次(n=4),其产量结果套其产量结果套L9(34)进进行整理后列表如
12、右,试予分析。行整理后列表如右,试予分析。一、数据整理一、数据整理 本例数据源自教材本例数据源自教材P133习题习题17C=T 2/nab=1170 2/36=38025 dfT=nab 1=4 9 1=35 SST=Y 2 C=7989 =412+49 2+19 2 38025SST=SS t SS e (分解分解)SSA SS B SS AB (再分解再分解)SS SS ABABABTt1111141492325138212221112252472313336222611654212347595040196522314338333615062312822181462731324835535
13、919583213553847441849332130332619108T1T2T3275408487529406235384376410396329445C=38025 SST=7989SSe=1824.5T1170SS1913.173633.552.67565.17(=6164.51)6164.50df2222(校对校对=8)8第二节第二节 完全随机试验数据完全随机试验数据二、二、SST的分解及的分解及SS t的再分解的再分解SS t Tt 2/n C 6164.5(138 2+722 +108 2)/4 C SS A TA2/nb C 1913.17(2752+408 2+4872)/1
14、2 38025 SS B TB2/na C 3633.5(5292+4062+235 2)/12 38025 SS T2/12 C 52.67(3842+376 2+410 2)/12 38025 SS T2/12 C 565.17(3962+329 2+4452)/12 38025SSe SST SS t 1824.5 用空列的用空列的SSAB实现实现“自动检校自动检校”:SS SS 52.67+565.17 617.84 SS t SSA SS B 617.83ABABABTt111114149232513821222111225247231333622261165421234759504
15、0196522314338333615062312822181462731324835535919583213553847441849332130332619108T1T2T3275408487529406235384376410396329445C=38025 SST=7989SSe=1824.5T1170SS1913.173633.552.67565.17(=6164.51)6164.50df2222(校对校对=8)8第二节第二节 完全随机试验数据完全随机试验数据 正交表在经典试验设计的复因素试验结果的方差分析过程中的作用:正交表在经典试验设计的复因素试验结果的方差分析过程中的作用:只要复
16、因素试验方案只要复因素试验方案“有表可套有表可套”,则不论因素多少,所有数据的整理都可用相应的正交表,则不论因素多少,所有数据的整理都可用相应的正交表一一表完成表完成,无需另外单列两向表计算主效应,无需另外单列两向表计算主效应SS和互作效应和互作效应SS;用正交表算得的空列用正交表算得的空列SS、DF或其累加值就是相应的互作效应或其累加值就是相应的互作效应SS、DF值,这一优点使得按值,这一优点使得按SSt 再分解原理算得的剩余再分解原理算得的剩余SS、DF值可与之值可与之“相互佐证相互佐证”,类似会计记帐时的,类似会计记帐时的“逼角逼角”,是正交表在是正交表在经典试验统计中所能起到的十分独特
17、的功用;经典试验统计中所能起到的十分独特的功用;用正交表参与复因素试验结果的方差分析过程体现了正交表与复因素试验的内在联系,认识用正交表参与复因素试验结果的方差分析过程体现了正交表与复因素试验的内在联系,认识这种必然性对正确理解正交表在建立多元统计模型中的作用非常重要,即:正交试验设计这种必然性对正确理解正交表在建立多元统计模型中的作用非常重要,即:正交试验设计回归正交设计回归正交设计回归旋转设计回归旋转设计(旋转组合设计旋转组合设计);有正交表参与方差分析过程时,处理平方和有正交表参与方差分析过程时,处理平方和SSt 的再分解步骤异常简洁、直观,本来烦琐的的再分解步骤异常简洁、直观,本来烦琐
18、的复因素试验结果的平方和再分解运算变成用几乎是同一个算式循环操作,其程序无须专业人复因素试验结果的平方和再分解运算变成用几乎是同一个算式循环操作,其程序无须专业人员就可编制,员就可编制,个人电脑应用个人电脑应用起来十分方便。起来十分方便。后续章节内容将告诉我们,不仅是完全随机试验,其它后续章节内容将告诉我们,不仅是完全随机试验,其它经典试验设计经典试验设计如随机单位组设计、如随机单位组设计、裂区设计、系统分组设计等得到的数据结果,只要在方案设计时稍加考虑,使得裂区设计、系统分组设计等得到的数据结果,只要在方案设计时稍加考虑,使得SSt 的再分的再分解解“有表可套有表可套”,同样可发挥正交表的功
19、用。这里再将一完全随机试验结果(教材,同样可发挥正交表的功用。这里再将一完全随机试验结果(教材P109)套)套正交表整理如下。正交表整理如下。效应效应ABABAB AB Tt11111122.0 26.5 24.472.921222230.0 27.5 26.083.531333332.4 26.5 27.085.941444430.5 27.0 25.182.652123423.5 25.8 27.076.362214333.2 28.5 30.191.872341238.0 35.5 33.0106.582432126.5 24.0 25.075.593134230.5 26.8 25.5
20、82.8103243136.5 34.0 33.5104.0113312428.0 30.5 24.683.1123421320.5 22.5 19.562.5134142334.5 31.4 29.395.2144231429.0 27.5 28.084.5154324127.5 26.3 28.582.3164413218.5 20.0 19.057.5T1324.9327.2305.3326.4334.7C=36680.4919SST=982.3181SSe=147.4133T:1326.9T2350.1363.8304.6337.3330.3T3332.4357.8328.7323.7
21、335.4SSt:834.9048T4319.5278.1388.3339.5326.5SS44.5106383.7356386.994015.35734.3073(校对(校对=834.9048)第二节第二节 完全随机试验数据完全随机试验数据二、二、SST的分解及的分解及SS t的再分解的再分解SS t Tt 2/n C 6164.5(138 2+722 +108 2)/4 CSS A TA2/nb C 1913.17(2752+408 2+4872)/12 38025 SS B TB2/na C 3633.5(5292+4062+235 2)/12 38025 SS T2/9 C 52.67
22、(3842+376 2+410 2)/12 38025SS T2/9 C 565.17(3962+329 2+4452)/12 38025SSe SST SS t 1824.5 用空列的用空列的SSAB实现实现“自动检校自动检校”:SS SS 52.67+565.17 617.84 SS t SSA SS B 617.83ABABABTt1111141492325138212221112252472313336222611654212347595040196522314338333615062312822181462731324835535919583213553847441849332130
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 生物 统计学 课件 方差分析

限制150内