正交試驗方差分析通俗易懂.doc

《正交試驗方差分析通俗易懂.doc》由會員分享，可在線閱讀，更多相關(guān)《正交試驗方差分析通俗易懂.doc（12頁珍藏版）》請在裝配圖網(wǎng)上搜索。

1、 第十一章 正交設(shè)計試驗資料的方差分析 在實(shí)際工作中 ，常常需要同時考察 3個或3個以上的試驗因素 ，若進(jìn)行全面試驗，則試驗的規(guī)模將很大 ，往往因試驗條件的限制而難于實(shí)施 。 正交設(shè)計是安排多因素試驗 、尋求最優(yōu)水平組合的一種 高效率試驗設(shè)計方法。 第一節(jié)、正交設(shè)計原理和方法 (一) 正交設(shè)計的基本概念 正 交 設(shè) 計 是利用正交表來安排多因素試驗、分析試驗結(jié)果的一種設(shè)計方法。它從多因素試驗的全部水平組合中挑選部分有代表性的水平組合進(jìn)行試驗，通過對這部分試驗結(jié)果的分析了解全面試驗的情況，找出最優(yōu)水平組合。 例如， 研究氮、磷、鉀肥施用量對某小麥品種

2、產(chǎn)量的影響： A因素是氮肥施用量，設(shè)A1、A2、A3 3個水平 ； B因素是磷肥施用量，設(shè)B1、B2、B3 3個水平 ； C因素是鉀肥施用量，設(shè)C1、C2、C3 3個水平。 這是一個3因素每個因素3水平的試驗 ，各因素的水平之間全部可能的組合有27種。 如果進(jìn)行全面試驗 ，可以分析各因素的效應(yīng) ，交互作用，也可選出最優(yōu)水平組合。 但全面試驗包含的水平組合數(shù)較多，工作量大 ，由于受試驗場地、經(jīng)費(fèi)等限制而難于實(shí)施 。 如果試驗的主要目的是尋求最優(yōu)水平組合，則可利用正交設(shè)計來安排試驗。 正交設(shè)計的基本特點(diǎn)是：用

3、部分試驗來代替全面試驗，通過對部分試驗結(jié)果的分析，了解全面試驗的情況。 正交試驗是用部分試驗來代替全面試驗，它不可能像全面試驗?zāi)菢訉Ω饕蛩匦?yīng)、交互作用一一分析；當(dāng)交互作用存在時，有可能出現(xiàn)交互作用的混雜。 如對于上述3因素每個因素3水平試驗，若不考慮交互作用，可利用正交表L9(34)安排，試驗方案僅包含9個水平組合，就能反映試驗方案包含27個水平組合的全面試驗的情況，找出最佳的生產(chǎn)條件。 一、正交設(shè)計的基本原理 表11-1 33試驗的全面試驗方案 正交設(shè)計就是從全面試驗點(diǎn)（水平組合）中挑選出有代表性的部分試驗點(diǎn)（水平組合）來進(jìn)行試驗。圖1中標(biāo)有‘9 ’

4、個試驗點(diǎn)，就是利用正交表L9(34)從27個試驗點(diǎn)中挑選出來的9個試驗點(diǎn)。即： (1)A1B1C1 (2)A1B2C2 (3)A1B3C3 (4)A2B1C2 (5)A2B2C3 (6)A2B3C1 (7)A3B1C3 (8)A3B2C1 (9)A3B3C2 上述選擇 ，保證了A因素的每個水平與B因素 、 C 因 素的各個水平在試驗中各搭配一次。 從圖1中可以看到，9個試驗點(diǎn)分布是均衡的 ，在立方體的每個平面上 有且僅有3個試驗點(diǎn)；每兩個平面的交線上有且僅有1個試驗點(diǎn)。 9個試驗點(diǎn)均衡地分布于整個立方體內(nèi) ，

5、有很強(qiáng)的代表性，能夠比較全面地反映全面試驗的基本情況。 二、正交表及其特性 (一) 正交表 表 11-2 是L8(27)正交表，其中 “L”代表正交表；L 右下角的數(shù)字“8”表示有8行，用這張正交表安排試驗包含8個處理 (水平組合) ；括號內(nèi)的底數(shù)“2” 表示因素的水平數(shù)，括號內(nèi) 2的指數(shù)“7”表示有7列，用這張正交表最多可以安排7個2水平因素。 表11-2 L8(27)正交表 2水平正交表還有L4(23)、L16(215)等； 3水平正交表有L9(34)、L27(313) 、…、 等。 (二) 正交表的特性

6、1、任一列中，不同數(shù)字出現(xiàn)的次數(shù)相同 例如L8(27)中不同數(shù)字只有1和2，它們各出現(xiàn)4次；L9(34)中不同數(shù)字有1、2和3，它們各出現(xiàn)3次 。 2、任兩列中，同一橫行所組成的數(shù)字對出現(xiàn)的次數(shù)相同 例如 L8(27)的任兩列中(1, 1), (1, 2), (2, 1), (2, 2)各出現(xiàn)兩次；L9(34)任兩列中 (1, 1), (1, 2), (1, 3), (2, 1), (2, 2), (2, 3), (3, 1), (3, 2), (3, 3)各出現(xiàn)1次。即每個因素的一個水平與另一因素的各個水平互碰次數(shù)相等，表明任意兩列各個數(shù)字之間的搭配是均勻的。

7、 用正交表安排的試驗，具有均衡分散和整齊可比的特點(diǎn)。 均衡分散，是指用正交表挑選出來的各因素 水 平 組合在全部水平組合中的分布是均衡的 。 由 圖11-1可以看出，在立方體中 ，任一平面內(nèi)都包含 3 個 試驗點(diǎn)， 任兩平面的交線上都包含1個試驗點(diǎn)。 整齊可比是指每一個因素的各水平間具有可比性。 因為正交表中每一因素的任一水平下都均衡地包含著另外因素的各個水平，當(dāng)比較某因素不同水平時，其它因素的效應(yīng)都彼此抵消。如在A、B、C 3個因素中，A因素的 3 個水平 A1、A2、A3 條件下各有 B、C 的 3 個不同水平，即： 在這9個水平組合中，A因素各水平下包括

8、了B、C因素的3個水平，雖然搭配方式不同，但B、C皆處于同等地位，當(dāng)比較A因素不同水平時，B因素不同水平的效應(yīng)相互抵消，C因素不同水平的效應(yīng)也相互抵消。所以A因素3個水平間具有可比性。同樣，B、C因素3個水平間亦具有可比性。 (三) 正交表的類別 1、相同水平正交表 各列中出現(xiàn)的最大數(shù)字相同的正交表稱為相同水平正交表。 L4(23)、L8(27)、L12(211)等各列中最大數(shù)字為2，稱為兩水平正交表； L9(34)、L27(313)等各列中最大數(shù)字為3，稱為3水平正交表。 2、 混合水平正交表 各列中出現(xiàn)的最大數(shù)字不完全相同的正交表稱為 混合水

9、平正交表。 L8(4124)表中有一列最大數(shù)字為4，有4列最大數(shù)字為2。 也就是說該表可以安排1個4水平因素和4個2水平因素。 L16(4423)，L16(4212)等都混合水平正交表。 三、正交設(shè)計方法 【例111】 某水稻栽培試驗選擇了3個水稻優(yōu)良品種(A)：二九矮、高二矮、窄葉青 ， 3種密度(B)： 15、20、25（萬苗/666.7m2）；3種施氮量(C)： 3、5、8（kg/666.7m2），試采用正交設(shè)計安排一個試驗方案。 (一) 確定試驗因素及其水平, 列出因素水平表 表11-3 因素水平表

10、 (二) 選用合適的正交表 根據(jù)因素、水平及需要考察的交互作用的多少來選擇合適的正交表。 選用正交表的原則是：既要能安排下試驗的全部因素(包括需要考查的交互作用)，又要使部分水平組合數(shù)（處理數(shù)）盡可能地少。 一般情況下，試驗因素的水平數(shù)應(yīng)恰好等于正交表記號中括號內(nèi)的底數(shù)；因素的個數(shù)（包括需要考查交互作用）應(yīng)不大于正交表記號中括號內(nèi)的指數(shù)；各因素及交互作用的自由度之和要小于所選 正交表 的 總 自由度，以便估計試驗誤差。 若各因素及交互作用的自由度之和等于所選正交表總自由度，則可采用有重復(fù)正交試驗來估計試驗誤差。 此例有3個3水平因素，若不考察交互作

11、用，則各因素自由度之和為因素個數(shù) (水平數(shù)-1) = 3 (3-1) =6，小于L9(34)總自由度 9-1=8，故可以選用L9(34)； 若要考察交互作用，則應(yīng)選用L27(313)，此時所安排的試驗方案實(shí)際上是全面試驗方案。 (三) 表頭設(shè)計 表頭設(shè)計就是把挑選出的因素和要考察的交互作用分別排入正交表的表頭適當(dāng)?shù)牧猩稀? 在不考察交互作用時，各因素可隨機(jī)安排在各列上；若考察交互作用，就應(yīng)按該正交表的交互作用列表安排 各 因 素與交互作用。 此例不考察交互作用，可將品種(A)、密度(B)和施氮量 (C)依次安排在L9(34)的第1、2、3列上，第4 列

12、 為空列，見表2-4。 表11-4 表頭設(shè)計 L9（34）表頭設(shè)計 L8(27) 表頭設(shè)計 (四) 列出試驗方案 把正交表中安排因素的各列(不包含欲考察的交互作用列)中的每個數(shù)字依次換成該因素的實(shí)際水平，就得到一個正交試驗方案。 表11-5 正交試驗方案 第二節(jié) 正交試驗資料的方差分析 若各號試驗處理都只有一個觀測值，則稱之為單個觀測值正交試驗； 若各號試驗處理都有兩個或兩個以上觀測值，則稱之為有重復(fù)觀測值正交試驗。 一、 單個觀測值正交試驗資料的方差分析 對【例11-1】用L9(34)安排試驗方案后，各號試驗只進(jìn)行一次，試

13、驗結(jié)果列于表2-6。試對其進(jìn)行方差分析。 表11-6 正交試驗結(jié)果計算表 Ti為各因素同一水平試驗指標(biāo)之和 ，T為9個試驗號的試驗指標(biāo)之和； 為各因素同一水平試驗指標(biāo)的平均數(shù)。 該試驗的9個觀測值總變異由A因素、B因素、C因素及誤差變異4部分組成，因而進(jìn)行方差分析時平方和與自由度的分解式為： SST = SSA + SSB + SSC+SSe dfT = dfA + dfB + dfC + dfe 用n表示試驗(處理)數(shù)；a、b、c表示A、B、C因素的水平數(shù)；ka、kb、kc表示A、B

14、、C因素的各水平重復(fù)數(shù)。本例，n=9、a=b=c=3、 ka=kb=kc=3。 1、計算各項平方和與自由度 矯正數(shù) C = T2/n = 37112/9 = 1530169.00 總平方和 SST =Σx2-C =（340.02+422.52+…+462.52） -1530169.00 =21238.00 A因素平方和 SSA=Σ/ka-C =(1201.52+1291.52+

15、1218.02)/3 -1530169.00 =1530.50 B因素平方和 SSB = Σ/kb-C =(1092.02+1278.52+1340.52)/3 -1530169.00 =11153.17 C因素平方和 SSC=Σ/kc-C =(1142.52+1245.02+1323.52)/3 -1530169.00 =5492.17 誤差平方和

16、 SSe=SST-SSA-SSB-SSC =21238.00-1530.5-11153.17 -5492.17 =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 = dfT-dfA-dfB-dfC = 8-2-2-2 = 2 2、列出方差分

17、析表，進(jìn)行F檢驗 表11-7 方差分析表 F 檢驗結(jié)果表明，三個因素對產(chǎn)量的影響都不顯著。究其原因可能是本例試驗誤差大且誤差自由度小(僅為2)，使檢驗的靈敏度低，從而掩蓋了考察因素的顯著性。 由于各因素對增重影響都不顯著，不必再進(jìn)行各因素水平間的多重比較。此時，可從表11-6中選擇平均數(shù)大的水平A2、B3、C3組合成最優(yōu)水平組合 A2B3C3。 若F檢驗結(jié)果3個因素對試驗指標(biāo)的影響顯著或極顯著，進(jìn)行各因素水平間多重比較常采用SSR法。 本例是選用相同水平正交表 L9(34)安排的試驗，A、B、C因素各水平重復(fù)數(shù)相同，即ka=kb=kc=3，它們的標(biāo)準(zhǔn)誤相同，即

18、 單個觀測值正交試驗資料的方差分析，其誤差是由“空列”來估計的。然而“空列”并不空，實(shí)際上是被未考察的交互作用所占據(jù)。 這種誤差既包含試驗誤差，也包含交互作用，稱為模型誤差。 若交互作用不存在，用模型誤差估計試驗誤差是可行的；若因素間存在交互作用，則模型誤差會夸大試驗誤差，有可能掩蓋考察因素的顯著性。 試驗誤差應(yīng)通過重復(fù)試驗值來估計。所以，進(jìn)行正交試驗最好能有二次以上的重復(fù)。正交試驗的重復(fù)，可采用完全隨機(jī)或隨機(jī)區(qū)組設(shè)計。 二、 有重復(fù)觀測值正交試驗資料的方差分析 【例114】 為了探討花生銹病藥劑防治效果的好壞，進(jìn)行了藥劑種類（A）、濃度（B）、劑

19、量（C）3因素試驗，各有3個水平，選用正交表L9(34)安排試驗。 試驗重復(fù)2次，隨機(jī)區(qū)組設(shè)計。正交試驗方案及試驗結(jié)果(產(chǎn)量 kg/小區(qū)，小區(qū)面積133.3m2)見表11—10，對試驗結(jié)果進(jìn)行方差分析。 用r表示試驗處理的重復(fù)數(shù)(區(qū)組數(shù))； n，a、b、c，ka、kb、kc的意義同上。 此例 r=2； n=9， a=b=c=3， ka=kb=kc=3。 表11-10 防治花生銹病藥劑種類、濃度、劑量正交試驗方案及結(jié)果計算表 Ti為各因素同一水平試驗指標(biāo)之和 ，T為9個試驗號的試驗指標(biāo)之和； 為各因素同一水平試驗指標(biāo)的平均數(shù)。 對于有重復(fù)、且重復(fù)采

20、用隨機(jī)區(qū)組設(shè)計的正交試驗，總變異可以劃分為處理間、區(qū)組間和誤差變異三部分，而處理間變異可進(jìn)一步劃分為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+SSr+SSe1+ SSe2

21、 dfT = dfA + dfB + dfC + dfr + dfe1 + dfe2 其中：SSr為區(qū)組間平方和；SSe1為模型誤差平方和；SSe2為試驗誤差平方和；SSt為處理間平方和； dfr 、 dfe1 、dfe2 、dft 為相應(yīng)自由度。 注意 ，對于重復(fù)采用完全隨機(jī)設(shè)計的正交試驗，在平方和與自由度劃分式中無 SSr、dfr項。 1、計算各項平方和與自由度 矯正數(shù) C =T2/ r n =549.02/（29）=16744.50 總平方和 SST=Σx2-C =28.02+35.02+…

22、+30.02-16744.50 =246.62 區(qū)組間平方和 SSr=ΣT2r /n-C =(273.52+275.52)/9- 16744.50 =0.22 處理間平方和 SSt = ΣT2t / r - C =(56.52+69.82+…+59.42)/2-16744.50 =245.96 A因素平方和 SSA = ΣT2A / kar - C = (

23、191.02+184.42+173.62)/(32) - 16744.50 =25.72 B因素平方和 SSB =ΣT2B / kbr - C =(191.42+169.72+187.92)/(32) - 16744.50 =45.24 C因素平方和 SSC = ΣT2C / kcr - C =(165.82+195.42+187.82)/(32) -16744.50 =78

24、.77 模型誤差平方和 SSe1 = SSt – SSA – SSB - SSC =245.96- 25.72- 45.24.- 78.77 = 96.23 試驗誤差平方和 SSe2 =SST – SSr - SSt =246.62- 0.22- 245.96 = 0.44 總自由度 dfT=rn-1=29-1=17 區(qū)組自由度 dfr=r-1=2-1=1 處理自由度 df

25、t=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 試驗誤差自由度 dfe2=dfT-dfr-dft =17-1-8 = 8 2、列出方差分析表，進(jìn)行F檢驗 表11-10 有重復(fù)觀測值正交試驗資料的方差分析表 首先

26、檢驗MSe1與MSe2差異的顯著性，若經(jīng)F檢驗不顯著，則可將其平方和與自由度分別合并，計算出合并的誤差均方，進(jìn)行F檢驗與多重比較，以提高分析的精度；若F檢驗顯著，說明存在交互作用 ，二者不能合并 ， 此時只能以MSe2進(jìn)行F檢驗與多重比較。 本例MSe1 / MSe2=802.00** ，模型誤差均方 MSe1 與試驗誤差均方 MSe2 差異極顯著，說明試驗因素間交互作用極顯著，只能以試驗誤差均方 MSe2 進(jìn)行F檢驗與多重比較。 F檢驗結(jié)果表明，藥劑種類（A）、濃度（B）、劑量（C）3 因素對花生產(chǎn)量都有極顯著影響；區(qū)組間差異不顯著 。 3、 多重比較 (1)

27、 若模型誤差顯著，說明試驗因素間存在交互作用，各因素所在列有可能出現(xiàn)交互作用的混雜，此時各試驗因素水平間的差異已不能真正反映因素的主效，因而進(jìn)行各因素水平間的多重比較無多大實(shí)際意義，但應(yīng)進(jìn)行試驗處理間的多重比較，以尋求最處理，即最優(yōu)水平組合。進(jìn)行各試驗處理間多重比較時選用試驗誤差均方MSe2。模型誤差顯著，還應(yīng)進(jìn)一步試驗，以分析因素間的交互作用。 (2) 若模型誤差不顯著 ，說明試驗因素間交互作用不顯著，各因素所在列有可能未出現(xiàn)交互作用的混雜，此時各因素水平間的差異能真正反映因素的主效，因而進(jìn)行各因素水平間的多重比較有實(shí)際意義，并從各因素水平間的多重比較中選出各因素的最優(yōu)水平相組合，得到最

28、優(yōu)水平組合。 進(jìn)行各因素水平間的多重比較時，用合并的誤差均方 MSe=（SSe1+ SSe2）/（dfe1+ dfe2） 此時可不進(jìn)行試驗處理間的多重比較。 本例模型誤差極顯著，說明因素間存在交互作用，不必進(jìn)行各因素水平間的多重比較，應(yīng)進(jìn)行試驗處理間的多重比較 ， 以尋求最處理，即最優(yōu)水平組合。為了讓讀者了解多重比較的方法，下面仍對各因素水平間、各試驗處理間進(jìn)行多重比較。 （1）A、B、C因素各水平平均數(shù)的多重比較 表11-12 A因素各水平平均數(shù)的多重比較表(SSR法) 表11-13 B因素各水平平均數(shù)的多重比較表(SSR法) 表1

29、1-14 C因素各水平平均數(shù)的多重比較表(SSR法) 因為 由dfe=8和k=2, 3, 查得SSR值并計算出LSR值列于表11-15。 表11-15 SSR值與LSR值表 多重比較結(jié)果表明： A因素各水平平均產(chǎn)量間 、B因素各水平平均產(chǎn)量間 、C因素各水平平均產(chǎn)量間差異顯著或極顯著。各因素的最優(yōu)水平為A1、B1、C2。 注意，本例模型誤差顯著，試驗因素間存在交互作用 ，不宜從各因素水平間的多重比較中選出各因素的最優(yōu)水平相組合來得到最優(yōu)水平組合。 （2）各試驗處理平均數(shù)間的多重比較 表11-16 各試驗處理平均數(shù)多重比較表(LSD法) 因為 由df

30、e=8, 查得t0.05(8)=2.306，t0.01(8)=3.355， 計算出LSD值為： LSD0.05= t0.05(8) =2.3060.245=0.565 LSD0.01= t0.01(8) =3.3550.245=0.822 各試驗處理間平均數(shù)多重比較結(jié)果，除第2號試驗處理與第7號試驗處理 、第3號試驗處理與第 6 號試驗處理平均產(chǎn)量差異不顯著外，其余各試驗處理平均產(chǎn)量間差異極顯著或顯著，最優(yōu)水平組合為第 2 號試驗處理A1B2C2（或第7號試驗處理A3B1C3） 本例模型誤差顯著，試驗因素間存在交互作用，應(yīng)以試驗處理間的多重比較尋求的最優(yōu)水平組合， 即第2號試驗處理 A1B2C2 （或第7號試驗處理 A3B1C3 ）為該試驗的最優(yōu)水平組合。