《《正交试验方差分析》PPT课件.ppt》由会员分享，可在线阅读，更多相关《《正交试验方差分析》PPT课件.ppt（74页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 正交设计试验正交设计试验资料的方差分析资料的方差分析 在实际工作中在实际工作中，常常需要同时考察，常常需要同时考察 3个或个或3个以上的试验因素个以上的试验因素，若进行全面试，若进行全面试验，则试验的规模将很大验，则试验的规模将很大，往往因试验条，往往因试验条件的限制而难于实施件的限制而难于实施。正交设计正交设计是安排多因素试验是安排多因素试验、寻求最寻求最优水平组合优水平组合的一种的一种 高效率试验设计方法。高效率试验设计方法。一、正交设计原理和方法 (一)正交设计的基本概念 正正 交交 设设 计计 是利用是利用正交表正交表来安排来安排多因素试验多因素试验、分分析试验结果析试验结果的一种设
2、计方法的一种设计方法。它从多因素试验的全部。它从多因素试验的全部水平组合中挑选部分有代表性的水平组合进行试验，水平组合中挑选部分有代表性的水平组合进行试验，通过对这部分试验结果的分析了解全面试验的情况通过对这部分试验结果的分析了解全面试验的情况 ，找出最优水平组合。找出最优水平组合。例如，例如，研究氮研究氮、磷、钾肥施用量、磷、钾肥施用量对某小麦品种产对某小麦品种产量的影响：量的影响：A因素是氮因素是氮肥施用量肥施用量，设，设A1、A2、A3 3个水平个水平；B因素是因素是磷肥施用量磷肥施用量，设，设B1、B2、B3 3个水平个水平；C因素是因素是钾肥施用量钾肥施用量，设，设C1、C2、C3
3、3个水平。个水平。这是一个这是一个3因素因素每个因素每个因素3水平的试验水平的试验，各因素的，各因素的水平之间全部可能的组合有水平之间全部可能的组合有27种。种。如果如果进行全面试验进行全面试验，可以分析各因素的，可以分析各因素的效应效应，交互作用，也可选出最优水平组合，交互作用，也可选出最优水平组合。但全面试验包含的水平组合数较多，工但全面试验包含的水平组合数较多，工作量大作量大，由于受试验场地、经费等限制而难，由于受试验场地、经费等限制而难于实施于实施。如果试验的如果试验的主要目的主要目的是寻求是寻求最优水平组最优水平组合合，则可利用正交设，则可利用正交设计来安排试验。来安排试验。正交设计
4、的正交设计的基本特点基本特点是：是：用部分试验来代替全面用部分试验来代替全面试验，通过对部分试验结果的分析，了解全面试验的试验，通过对部分试验结果的分析，了解全面试验的情况。情况。正交试验是用部分试验来代替全面试验，它不可正交试验是用部分试验来代替全面试验，它不可能像全面试验那样对各因素效应、交互作用一一分析；能像全面试验那样对各因素效应、交互作用一一分析；当交互作用存在时，有可能出现交互作用的混杂当交互作用存在时，有可能出现交互作用的混杂。表表1 1 33试验试验的全面的全面试验试验方案方案C1C2C3A1B1A1B1C1A1B1C2A1B1C3B2A1B2C1A1B2C2A1B2C3B3A
5、1B3C1A1B3C2A1B3C3A2B1A2B1C1A2B1C2A2B1C3B2A2B2C1A2B2C2A2B2C3B3A2B3C1A2B3C2A2B3C3A3B1A3B1C1A3B1C2A3B1C3B2A3B2C1A3B2C2A3B2C3B3A3B3C1A3B3C2A3B3C3 图图图图1 31 3因素每个因素因素每个因素因素每个因素因素每个因素3 3水平试验点的均衡分布图水平试验点的均衡分布图水平试验点的均衡分布图水平试验点的均衡分布图 正正交交设设计计就就是是从从全全面面试试验验点点（水水平平组组合合）中中挑挑选选出出有有代代表表性性的的部部分分试试验验点点（水水平平组组合合）来来进进
6、行行试试验验。图图1中中标标有有 9个个试试验验点点，就就是是利利用用正正交交表表L9(34)从从27个试验点中挑选出来的个试验点中挑选出来的9个试验点。即：个试验点。即：(1)A1B1C1 (2)A1B2C2 (3)A1B3C3 (4)A2B1C2 (5)A2B2C3 (6)A2B3C1 (7)A3B1C3 (8)A3B2C1 (9)A3B3C2 上述选择上述选择，保证了，保证了A因素的每个水平与因素的每个水平与B因素因素、C 因因 素的各个水平在试验中各搭配一次。素的各个水平在试验中各搭配一次。从图从图1中可以看到，中可以看到，9个个试验点分布是均衡的试验点分布是均衡的，在，在立方体的每个
7、平面上立方体的每个平面上 有且仅有有且仅有3个试验点；每两个平个试验点；每两个平面的交线上有且仅有面的交线上有且仅有1个试验点。个试验点。9个试验点均衡地分布于整个立方体内个试验点均衡地分布于整个立方体内，有很强，有很强的代表性，能够比较全面地反映全面试验的基本情况。的代表性，能够比较全面地反映全面试验的基本情况。正交表是一种特别的表格，是正交设计的基本工具。正交表是一种特别的表格，是正交设计的基本工具。我们只介绍它的记号、特点和使用方法。我们只介绍它的记号、特点和使用方法。二、正交表及其特性,(一)正交表 记号及含义记号及含义 正交表的列数正交表的列数（最多能安排的因素个数，（最多能安排的因
8、素个数，包括交互作用、误差等）包括交互作用、误差等）正交表的行数正交表的行数（需要做的试验次数）（需要做的试验次数）各因素的水平数各因素的水平数（各因素的水平数相等）各因素的水平数相等）q正交表正交表的代号的代号如如 表示表示 表示各因素的表示各因素的水平数水平数为为2，做做8次试验次试验，最多考虑，最多考虑7个个因素因素（含交互作用）的（含交互作用）的正正交表交表。表表2 L8(27)正交表正交表 2水平正交表还有水平正交表还有L4(23)、L16(215)等；等；3水平正交表有水平正交表有L9(34)、L27(313)、等。等。(二二)正交表的特性正交表的特性 1 1、任一列中，不同数字出
9、现的次数相同、任一列中，不同数字出现的次数相同 例如例如L8(27)中不同数字只有中不同数字只有1和和2，它们各出现，它们各出现4次；次；L9(34)中不同数字有中不同数字有1、2和和3，它们各出现，它们各出现3次次。2 2、任两列中，同一横行所组成的数字对出、任两列中，同一横行所组成的数字对出现的次数相同现的次数相同 例如例如 L8(27)的任两列中的任两列中(1,1),(1,2),(2,1),(2,2)各出现两次；各出现两次；L9(34)任两列中任两列中 (1,1),(1,2),(1,3),(2,1),(2,2),(2,3),(3,1),(3,2),(3,3)各出现各出现1次。即次。即每个
10、因素的一个水平与另一因素的各个水平互碰次数每个因素的一个水平与另一因素的各个水平互碰次数相等，表明任意两列各个数字之间的相等，表明任意两列各个数字之间的搭配是均匀搭配是均匀的。的。用正交表安排的试验，具有用正交表安排的试验，具有均衡分散均衡分散和和整齐可比整齐可比的特点。的特点。均衡分散均衡分散，是指用正交表挑选出来的各因素是指用正交表挑选出来的各因素 水水 平平 组合在全部水平组合中的分布是均衡的组合在全部水平组合中的分布是均衡的。由由 图图1可以看出，在立方体中可以看出，在立方体中，任一平面内都包含，任一平面内都包含 3 个个 试试验点，验点，任两平面的交线上都包含任两平面的交线上都包含1
11、个试验点。个试验点。整齐可比整齐可比是指每一个因素的各水平间具是指每一个因素的各水平间具有可比性。有可比性。因为正交表中每一因素的任一水平下都因为正交表中每一因素的任一水平下都均衡地包含着另外因素的各个水平，当比较均衡地包含着另外因素的各个水平，当比较某因素不同水平时，其它因素的效应都彼此某因素不同水平时，其它因素的效应都彼此抵消。如在抵消。如在A、B、C 3个因素中，个因素中，A因素的因素的 3 个水平个水平 A1、A2、A3 条件下各有条件下各有 B、C 的的 3 个不同水平，即：个不同水平，即：在在这这9个个水水平平组组合合中中，A因因素素各各水水平平下下包包括括了了B、C因因素素的的3
12、个个水水平平，虽虽然然搭搭配配方方式式不不同同，但但B、C皆皆处处于于同同等等地地位位，当当比比较较A因因素素不不同同水水平平时时，B因因素素不不同同水水平平的的效效应应相相互互抵抵消消，C因因素素不不同同水水平平的的效效应应也也相相互互抵抵消消。所所以以A因因素素3个个水水平平间间具具有有可可比比性性。同同样样，B、C因因素素3个水平间亦具有可比性。个水平间亦具有可比性。(三)正交表的类别 1 1、相相同同水水平平正正交交表表 各各列列中中出出现现的的最最大大数数字字相相同同的正交表称为相同水平正交表。的正交表称为相同水平正交表。L4(23)、L8(27)、L12(211)等等各各列列中中最
13、最大大数数字字为为2，称为两水平正交表；称为两水平正交表；L9(34)、L27(313)等等各各列列中中最最大大数数字字为为3，称称为为3水水平正交表。平正交表。2、混合水平正交表混合水平正交表 各列中出现的最各列中出现的最大数字不完全相同的正交表称为大数字不完全相同的正交表称为 混合水平正混合水平正交表。交表。L8(4124)表中有一列最大数字为表中有一列最大数字为4，有，有4列最大数字为列最大数字为2。也就是说该表可以安排也就是说该表可以安排1个个4水平因素和水平因素和4个个2水平因素。水平因素。L16(4423)，L16(4212)等都混合水平正等都混合水平正交表。交表。三、正交设计方法
14、 【例例1】某水稻栽培试验选择了某水稻栽培试验选择了3个水稻优良个水稻优良品种品种(A)：二九矮、高二矮、窄叶青：二九矮、高二矮、窄叶青，3种种密度密度(B)：15、20、25（万苗（万苗2）；）；3种种施氮量施氮量(C)：3、5、8（2），），试采用正交设计安排一个试验方案。试采用正交设计安排一个试验方案。(一一)确定试验因素及其水平确定试验因素及其水平,列出因素水列出因素水平表平表 水水 平平因因 素素品种品种(A)密度密度(B)施氮量施氮量(C)1二九矮二九矮(A1)15(B1)3(C1)2高二矮高二矮(A2)20(B2)5(C2)3窄叶青窄叶青(A3)25(B3)8(C3)表表3 因素
15、水平表因素水平表 (二)选用合适的正交表 根根据据因因素素、水水平平及及需需要要考考察察的的交交互互作作用用的的多多少少来选择合适的正交表。来选择合适的正交表。选选用用正正交交表表的的原原则则是是：既既要要能能安安排排下下试试验验的的全全部部因因素素(包包括括需需要要考考查查的的交交互互作作用用)，又又要要使使部部分分水水平平组合数（处理数）尽可能地少。组合数（处理数）尽可能地少。一般情况下，试验因素的水平数应恰好一般情况下，试验因素的水平数应恰好等于正交表记号中括号内的底数；因素的个等于正交表记号中括号内的底数；因素的个数（包括需要考查交互作用）应不大于正交数（包括需要考查交互作用）应不大于
16、正交表记号中括号内的指数；表记号中括号内的指数；各因素及交互作用各因素及交互作用的自由度之和要小于所选的自由度之和要小于所选 正交表正交表 的的 总总 自由自由度度，以便估计试验误差。，以便估计试验误差。若各因素及交互作用的自由度之和等于若各因素及交互作用的自由度之和等于所选正交表总自由度，则可所选正交表总自由度，则可采用有重复正交采用有重复正交试验试验来估计试验误差。来估计试验误差。此此例例有有3个个3水水平平因因素素，若若不不考考察察交交互互作作用用，则则各各因因素素自自由由度度之之和和为为因因素素个个数数(水水平平数数-1)=3 (3-1)=6，小小于于L9(34)总总自自由由度度 9-
17、1=8，故故可可以以选选用用L9(34)；若若要要考考察察交交互互作作用用，则则应应选选用用L27(313)，此此时时所所安排的试验方案实际上是安排的试验方案实际上是全面试验方案全面试验方案。(三)表头设计 表头设计就是把挑选出的因素和要考察表头设计就是把挑选出的因素和要考察的交互作用分别排入正交表的表头适当的列的交互作用分别排入正交表的表头适当的列上。上。在不考察交互作用时，各因素可随机安在不考察交互作用时，各因素可随机安排在各列上；若考察交互作用，就应按该排在各列上；若考察交互作用，就应按该正正交表的交互作用列表交表的交互作用列表安排安排 各各 因因 素与交互作素与交互作用。用。表表2 L
18、8(27)正交表正交表 L8(27)表表头设计头设计因素数因素数列列 号号12345673ABABCACBC4ABABCDCACBDBCADD4ABCDABCBDACDBCAD5ADEBCDABCECBDACBEDAEBCEAB表表4 表头设计表头设计列列 号号1234因因 素素ABC空空 此例不考察交互作用，可将品种此例不考察交互作用，可将品种(A)、密、密度度(B)和施氮量和施氮量(C)依次安排在依次安排在L9(34)的第的第1、2、3列上，第列上，第4 列列 为空列，见表为空列，见表4。(四)列出试验方案 把正交表中安排因素的各列把正交表中安排因素的各列(不包含欲考察的不包含欲考察的交互
19、作用列交互作用列)中的每个数字依次换成该因素的实际中的每个数字依次换成该因素的实际水平，就得到一个水平，就得到一个正交试验方案正交试验方案。表表5 正交试验方案正交试验方案试验号因 素ABC12311(二九矮)1(15)1(3)21(二九矮)2(20)2(5)31(二九矮)3(25)3(8)42(高二矮)1(15)2(5)52(高二矮)2(20)3(8)62(高二矮)3(25)1(3)73(窄叶青)1(15)3(8)83(窄叶青)2(20)1(3)93(窄叶青)3(25)2(5)第二节 正交试验资料的方差分析 若若各各号号试试验验处处理理都都只只有有一一个个观观测测值值，则则称称之之为为单个观
20、测值正交试验单个观测值正交试验；若若各各号号试试验验处处理理都都有有两两个个或或两两个个以以上上观观测测值值，则称之为则称之为有重复观测值正交试验有重复观测值正交试验。一、单个观测值正交试验资料的方差分析 对对【例例1】用用L9(34)安安排排试试验验方方案案后后，各各号号试试验验只只进行一次，试验结果列于表进行一次，试验结果列于表6。试对其进行方差分析。试对其进行方差分析。表表6 正交试验结果计算表正交试验结果计算表试验试验号号因因 素素产产量量ABC(1)(2)(3)1111340.0(x1)2122422.5(x2)3133439.0(x3)4212360.0(x4)5223492.5(
21、x5)6231439.0(x6)7313392.0(x7)8321363.5(x8)9332462.5(x9)T11201.51092.01142.53711.0(T)T21291.51278.51245.0T31218.01340.51323.5 400.50364.00380.83430.50426.17415.00406.00446.83441.17 Ti为各因素同一水平试验指标之和为各因素同一水平试验指标之和，T为为9个试个试验号的试验指标之和验号的试验指标之和；为为各因素同一水平各因素同一水平试验试验指指标标的平均数的平均数。该该试试验验的的9个个观观测测值值总总变变异异由由A因因素
22、素、B因因素素、C因因素素及及误误差差变变异异4部部分分组组成成，因因而而进进行行方方差差分分析析时时平平方方和和与自由度的分解式为：与自由度的分解式为：SST=SSA+SSB+SSC+SSe dfT =dfA +dfB +dfC +dfe 用用n表表示示试试验验(处处理理)数数；a、b、c表表示示A、B、C因因素素的的水水平平数数；ka、kb、kc表表示示A、B、C因因素素的的各各水水平平重复数。本例，重复数。本例，n=9、a=b=c=3、ka=kb=kc=3。1 1、计算各项平方和与自由度、计算各项平方和与自由度 矫正数矫正数 C=T2/n=37112/9=总平方和总平方和 SST=x2-
23、C =（222）A因素平方和因素平方和 SSA A=/ka-C =222222)/3)/3 BB因素平方和因素平方和因素平方和因素平方和 SSSSB B =/k/kb b-C-C 222222)/3)/3 C因素平方和因素平方和 SSC=T2C/kc-C 222)/3 =5492.17 误差平方和误差平方和 SSe=SST-SSA-SSB-SSC =3062.16 总自由度总自由度 dfT =n-1=9-1=8 A因素自由度因素自由度 dfA=a-1=3-1=2 B因素自由度因素自由度 dfB=b-1=3-1=2 C因素自由度因素自由度 dfC=c-1=3-1=2 误差自由度误差自由度 dfe
24、=dfT-dfA-dfB-dfC =8-2-2-2=2 2 2、列出方差分析表，进行、列出方差分析表，进行F F检验检验表表7 方差分析表方差分析表变变异来源异来源SSSSdfdfMSMSF FF0.05(2,2)品种品种(A)(A)1530.501530.502 2765.25765.2511 19.0019.00密度密度(B)(B)11153.1711153.172 25576.595576.593.643.64nsns施氮量施氮量(C)(C)5492.175492.172 22746.092746.091.791.79 ns ns误误差差3062.163062.162 21531.081
25、531.08总变总变异异21238.0021238.008 8表表6 正交试验结果计算表正交试验结果计算表试验试验号号因因 素素产产量量ABC(1)(2)(3)1111340.0(x1)2122422.5(x2)3133439.0(x3)4212360.0(x4)5223492.5(x5)6231439.0(x6)7313392.0(x7)8321363.5(x8)9332462.5(x9)F 检检验验结结果果表表明明，三三个个因因素素对对产产量量的的影影响响都都不不显显著著。究究其其原原因因可可能能是是本本例例试试验验误误差差大大且且误误差差自自由由度度小小(仅仅为为2)，使使检检验验的的灵
26、灵敏敏度度低低，从从而而掩掩盖盖了了考考察察因因素素的的显著性。显著性。由由于于各各因因素素对对增增重重影影响响都都不不显显著著，不不必必再再进进行行各各因因素素水水平平间间的的多多重重比比较较。此此时时，可可从从表表6中中选选择择平平均均数数大的水平大的水平A2、B3、C3组合成组合成最优水平组合最优水平组合 A2B3C3。若若F检验结检验结果果3个因素个因素对试验对试验指指标标的影响的影响显显著或极著或极显显著，著，进进行各因素水平行各因素水平间间多重比多重比较较常采用常采用SSR法。法。本例是本例是选选用相同水平正交表用相同水平正交表 L9(34)安排安排的的试验试验，A、B、C因素各水
27、平重复数相因素各水平重复数相同，同，即即ka=kb=kc=3 3，它，它们们的的标标准准误误相相同，即同，即 单单个个观观测测值值正正交交试试验验资资料料的的方方差差分分析析，其其误误差差是是由由“空空列列”来来估估计计的的。然然而而“空空列列”并并不不空空，实实际际上上是被未考察的交互作用所占据。是被未考察的交互作用所占据。这这种种误误差差既既包包含含试试验验误误差差，也也包包含含交交互互作作用用，称称为为模型误差模型误差。若若交交互互作作用用不不存存在在，用用模模型型误误差差估估计计试试验验误误差差是是可可行行的的；若若因因素素间间存存在在交交互互作作用用，则则模模型型误误差差会会夸夸大大
28、试验误差，有可能掩盖考察因素的显著性试验误差，有可能掩盖考察因素的显著性。试验误差应通过重复试验值来估试验误差应通过重复试验值来估计。计。所以，进行正交试验最好能有二所以，进行正交试验最好能有二次以上的重复。正交试验的重复，可次以上的重复。正交试验的重复，可采用采用完全随机或随机区组设计完全随机或随机区组设计。二二、有重复观测值正交试验资料的方差分析有重复观测值正交试验资料的方差分析 【例例2】为了探讨花生锈病药剂防治效果的好坏，为了探讨花生锈病药剂防治效果的好坏，进行了进行了药剂种类药剂种类（A）、）、浓度（浓度（B）、剂量）、剂量（C）3因因素试验，各有素试验，各有3个水平，选用正交表个水
29、平，选用正交表L9(34)安排试验。安排试验。试验重复试验重复2次，随机区组设计。正交试验方案及试验结次，随机区组设计。正交试验方案及试验结果果(产量产量 kg/小区，小区面积小区，小区面积2)见表见表10，对试验结果进，对试验结果进行方差分析。行方差分析。用用r表示试验处理的重复数表示试验处理的重复数(区组数区组数)；n，a、b、c，ka、kb、kc的意义同上的意义同上。此例此例 r=2；n=9，a=b=c=3，ka=kb=kc=3。表表10 10 防治花生锈病药剂种类、浓度、剂量正交试验方案及结果计算表防治花生锈病药剂种类、浓度、剂量正交试验方案及结果计算表试验试验号号 因因 素素 产量产
30、量x(kg/小区小区)TtABC区区组组I区区组组II(1)(2)(3)11（百菌清）（百菌清）1（高）（高）1（80）28.028.556.528.2521（百菌清）（百菌清）2（中）（中）2（100）35.034.869.834.9031（百菌清）（百菌清）3（低）（低）3（120）32.232.564.732.3542（敌锈敌锈灵）灵）1（高）（高）2（100）33.033.266.233.1052（敌锈敌锈灵）灵）2（中）（中）3（120）27.427.054.427.2062（敌锈敌锈灵）灵）3（低）（低）1（80）31.832.063.831.9073（波（波尔尔多）多）1（高）（
31、高）3（120）34.234.568.734.3583（波（波尔尔多）多）2（中）（中）1（80）22.523.045.522.7593（波（波尔尔多）多）3（低）（低）2（100）29.430.059.429.70T1191.0191.4165.8 273.5275.5549.0T2184.4169.7195.4 T3173.6187.9187.8 31.8331.9027.6330.7328.2832.5728.9331.3231.30 Ti为各因素同一水平试验指标之和为各因素同一水平试验指标之和，T为为9个试个试验号的试验指标之和验号的试验指标之和；为为各因素同一水平各因素同一水平试验试
32、验指指标标的平均数的平均数。对对于于有有重重复复、且且重重复复采采用用随随机机区区组组设设计计的的正正交交试试验验，总总变变异异可可以以划划分分为为处处理理间间、区区组组间间和和误误差差变变异异三三部部分分，而而处处理理间间变变异异可可进进一一步步划划分分为为A因因素素、B因因素素、C因因素素与与模模型型误误差差变变异异四四部部分分。此此时时，平平方方和和与与自自由由度分解式为：度分解式为：SST=SSt+SSr+SSe2 dfT=dft +dfr+dfe2 而而 SSt=SSA+SSB+SSC+SSe1 dft =dfA+dfB+dfC+dfe1 于是于是 SST=SSA+SSB+SSC+S
33、Sr+SSe1+SSe2 dfT=dfA+dfB+dfC+dfr+dfe1+dfe2 其中：其中：SSr为区组间平方和；为区组间平方和；SSe1为为模型误差平方和；模型误差平方和；SSe2为试验误差平方和；为试验误差平方和；SSt为处理间平方和；为处理间平方和；dfr、dfe1、dfe2 、dft 为相应自由度。为相应自由度。注意注意，对于重复采用完全随机设计的正交试验，对于重复采用完全随机设计的正交试验，在平方和与自由度划分式中无在平方和与自由度划分式中无 SSr、dfr项。项。1 1、计算各项平方和与自由度、计算各项平方和与自由度 矫正数矫正数 C =T2/r n=2/（29）总平方和总平
34、方和 SST=x2-C =222 区组间平方和区组间平方和 SSr=T2r/n-C =22 处理间平方和处理间平方和 SSt =T2t/r-C 222 A因素平方和因素平方和 SSA =T2A/kar-C =222)/(32)=B B因素平方和因素平方和因素平方和因素平方和 SSSSB B=T T2 2B B/k/kb br-Cr-C 222222 =C C因素平方和因素平方和因素平方和因素平方和 SSSSC C=T T2 2C C/k/kc cr-Cr-C =222222 =模型误差平方和模型误差平方和模型误差平方和模型误差平方和 SSSSe1 e1=SSSSt t SS SSA A SS
35、SSB B-SS-SSC C =试验误差平方和试验误差平方和试验误差平方和试验误差平方和 SSSSe2 e2=SSSST T SS SSr r-SS-SSt t =总自由度总自由度 dfT=rn-1=29-1=17 区组自由度区组自由度 dfr=r-1=2-1=1 处理自由度处理自由度 dft=n-1=9-1=8 A因素自由度因素自由度 dfA=a-1=3-1=2 B因素自由度因素自由度 dfB=b-1=3-1=2 C因素自由度因素自由度 dfC=c-1=3-1=2 模型误差自由度模型误差自由度 dfe1=dft-dfA-dfB-dfC =8-2-2-2=2 试验误差自由度试验误差自由度 df
36、e2=dfT-dfr-dft=17-1-8=8 2 2、列出方差分析表，进行、列出方差分析表，进行F F检验检验表表表表10 10 10 10 有重复观测值正交试验资料的方差分析表有重复观测值正交试验资料的方差分析表有重复观测值正交试验资料的方差分析表有重复观测值正交试验资料的方差分析表变变异来源异来源SSdfMSFF0.05F0.01 A25.72212.86214.33*4.107.55 B45.24222.62377.00*C78.77239.39656.50*区区组组0.2210.223.67ns4.9610.01模型模型误误差差(e1)96.23248.12802.00*试验误试验误
37、差差(e2)0.4480.06总总 的的246.6217 首先检验首先检验MSe1与与MSe2差异的显著性，若经差异的显著性，若经F检验检验不显著，则可将其平方和与自由度分别合并，计算出不显著，则可将其平方和与自由度分别合并，计算出合并的误差均方，进行合并的误差均方，进行F检验与多重比较，以提高分析检验与多重比较，以提高分析的精度；若的精度；若F检验显著，说明存在交互作用检验显著，说明存在交互作用 ，二者不，二者不能合并能合并 ，此时只能以此时只能以MSe2进行进行F检验与多重比较。检验与多重比较。本例本例MSe1/MSe2=802.00*，模型误差，模型误差均方均方 MSe1 与试验误差均方
38、与试验误差均方 MSe2 差异极显差异极显著，著，说明试验因素间交互作用极显著说明试验因素间交互作用极显著，只能，只能以试验误差均方以试验误差均方 MSe2 进行进行F检验与多重比检验与多重比较。较。F检验结果表明，检验结果表明，药剂种类药剂种类（A）、）、浓浓度度（B）、）、剂量剂量（C）3 因素对花生产量都因素对花生产量都有极显著影响；区组间差异不显著有极显著影响；区组间差异不显著。3、多重比较多重比较 (1)若若模模型型误误差差显显著著，说说明明试试验验因因素素间间存存在在交交互互作作用用，各各因因素素所所在在列列有有可可能能出出现现交交互互作作用用的的混混杂杂，此此时时各各试试验验因因
39、素素水水平平间间的的差差异异已已不不能能真真正正反反映映因因素素的的主主效效，因因而而进进行行各各因因素素水水平平间间的的多多重重比比较较无无多多大大实实际际意意义义，但但应应进进行行试试验验处处理理间间的的多多重重比比较较，以以寻寻求求最最优优水水平平组组合合。进进行行各各试试验验处处理理间间多多重重比比较较时时选选用用试试验验误误差差均均方方MSe2。模模型型误误差差显显著著，还还应应进进一一步步试试验验，以以分分析析因因素素间的交互作用。间的交互作用。(2)2)若若模型误差不显著模型误差不显著 ，说明试验因说明试验因素间交互作用不显著，各因素所在列有可能素间交互作用不显著，各因素所在列有
40、可能未出现交互作用的混杂，此时各因素水平间未出现交互作用的混杂，此时各因素水平间的差异能真正反映因素的主效，因而进行各的差异能真正反映因素的主效，因而进行各因素水平间的多重比较有实际意义，并从各因素水平间的多重比较有实际意义，并从各因素水平间的多重比较中选出各因素的最优因素水平间的多重比较中选出各因素的最优水平相组合，得到最优水平组合。水平相组合，得到最优水平组合。进行各因素水平间的多重比较时，用合并进行各因素水平间的多重比较时，用合并的误差均方的误差均方 MSe=（SSe1+SSe2）/（dfe1+dfe2）此时可不进行试验处理间的多重比较。此时可不进行试验处理间的多重比较。本例模型误差极显
41、著，说明因素间存在交本例模型误差极显著，说明因素间存在交互作用，不必进行各因素水平间的多重比较，互作用，不必进行各因素水平间的多重比较，应进行试验处理间的多重比较应进行试验处理间的多重比较 ，以寻求最优以寻求最优水平组合。水平组合。为了让读者了解多重比较的方法，为了让读者了解多重比较的方法，下面仍对各因素水平间、各试验处理间进行多下面仍对各因素水平间、各试验处理间进行多重比较。重比较。（1）A、B、C因素各水平平均数的多重比较因素各水平平均数的多重比较 表表12 A12 A因素各水平平均数的多重比因素各水平平均数的多重比较较表表(SSRSSR法法)A因素因素平均数平均数-28.93-30.73
42、A131.832.90*1.10*A230.731.80*A328.93 表表13 B13 B因素各水平平均数的多重比因素各水平平均数的多重比较较表表(SSRSSR法法)B因素因素平均数平均数 -28.28 -31.32B131.903.62*0.58*B331.323.04*B228.28 表表14 C14 C因素各水平平均数的多重比因素各水平平均数的多重比较较表表(SSRSSR法法)C因素因素平均数平均数 -27.63 -31.30 C232.574.94*1.27*C331.303.67*C127.63因因为为 由由dfe=8和和k=2,3,查得查得SSR值并计算出值并计算出LSR值列于
43、表值列于表15。表表15 SSR值与值与LSR值表值表dfekSSR0.05SSR0.01LSR0.05LSR0.01823.264.740.330.4733.395.000.340.50 多重比较结果表明：多重比较结果表明：A因素各水平平均因素各水平平均产量间产量间、B因素各水平平均产量间因素各水平平均产量间、C因素因素各水平平均产量间差异显著或极显著。各因各水平平均产量间差异显著或极显著。各因素的最优水平为素的最优水平为A1、B1、C2。注意，本例模型误差显著，试验因素间注意，本例模型误差显著，试验因素间存在交互作用存在交互作用，不宜从各因素水平间的多重，不宜从各因素水平间的多重比较中选出
44、各因素的最优水平相组合来得到比较中选出各因素的最优水平相组合来得到最优水平组合。最优水平组合。表表11-16 11-16 各各试验处试验处理平均数多重比理平均数多重比较较表表(LSDLSD法法)试验试验号号平均数平均数-22.75-27.20-28.25-29.70-31.90-32.35-33.10-34.35234.9012.15*7.70*6.65*5.20*3.00*2.55*1.80*0.55734.3511.60*7.15*6.10*4.65*2.45*2.00*1.25*433.1010.35*5.90*4.85*3.40*1.20*0.75*332.359.60*5.15*4.
45、10*2.65*0.45631.909.15*4.70*3.65*2.20*929.706.95*2.50*1.45*128.255.50*1.05*527.204.45*822.75（2）各试验处理平均数间的多重比较）各试验处理平均数间的多重比较 因为因为 由由dfe=8,查得查得t0.05(8)=，t0.01(8)，计算出计算出LSD值为：值为：LSD=t0.05(8)=2.3060.245=0.565 LSD=t0.01(8)=3.3550.245=0.822 各试验处理间平均数多重比较结果，除各试验处理间平均数多重比较结果，除第第2号试验处理与第号试验处理与第7号试验处理号试验处理、第
46、、第3号试号试验处理与第验处理与第 6 号试验处理平均产量差异不显号试验处理平均产量差异不显著外，其余各试验处理平均产量间差异极显著外，其余各试验处理平均产量间差异极显著或显著，著或显著，最优水平组合为第最优水平组合为第 2 号试验处理号试验处理A1B2C2（或（或第第7号试验处理号试验处理A3B1C3）本例模型误差显著，试验因素间本例模型误差显著，试验因素间存在交互作用，应以试验处理间的多存在交互作用，应以试验处理间的多重比较寻求的最优水平组合，重比较寻求的最优水平组合，即第即第2 2号试验处理号试验处理 A A1 1B B2 2C C2 2（或第（或第7 7号试验处号试验处理理 A A3 3B B1 1C C3 3）为该试验的最优水平组）为该试验的最优水平组合。合。