經(jīng)濟學第七章方差分析與正交試驗設計初步課件

1、第七章方差分析與正交試驗設計初步7/27/20221版權(quán)所有 BY 統(tǒng)計學課程組本章重點與難點重點: 是了解和掌握方差分析與正交試驗設計的思想和方法。難點: 方差分析與正交試驗設計的方法。7/27/20222版權(quán)所有 BY 統(tǒng)計學課程組學習目標 通過本章的學習，正確理解方差分析與正交試驗設計的思想和有關概念，了解各種方法的應用條件和范圍，能夠應用軟件進行有關的計算，能說明計算結(jié)果的含義。7/27/20223版權(quán)所有 BY 統(tǒng)計學課程組案 例 染整工藝對布的縮水率有影響嗎 ? 某公司為考察染整工藝是否對布的縮水率有顯著影響進行了試驗。試驗中采用三種染整工藝，分別對四種布樣進行了處理，測得縮水率的

2、百分比資料如7-1表所示:7/27/20224版權(quán)所有 BY 統(tǒng)計學課程組問題：(1) 如何判斷染整工藝類型對縮水率是否有顯著影響？若影響顯著，應采用哪種工藝？(2) 如果還有一個因素B（試驗配方，共三種）也對縮水率產(chǎn)生作用，應該如何選擇最優(yōu)方案？(3) 如果還有兩個或者兩個以上的因素也對縮水率產(chǎn)生作用，應該如何選擇最優(yōu)方案？ 為了回答上述問題，本章將介紹單因素方差分析、雙因素方差分析和正交試驗設計法的有關概念、分析思想、數(shù)學模型、方差分析表和正交試驗設計的基本方法等。7/27/20225版權(quán)所有 BY 統(tǒng)計學課程組第一節(jié) 方差分析的基本思想 第二節(jié) 單因素方差分析 第四節(jié) 雙因素方差分析 第

3、四節(jié) 正交試驗設計初步 第七章方差分析與正交試驗設計初步7/27/20226版權(quán)所有 BY 統(tǒng)計學課程組第一節(jié) 方差分析的基本思想 一、方差分析的有關概念和基本思想 二、顯著性檢驗7/27/20227版權(quán)所有 BY 統(tǒng)計學課程組一、方差分析的有關概念和基本思想 （一）簡單平均數(shù) 方差分析（Analysis of Variance，ANOVA）是一種檢驗多個總體均值是否相等的統(tǒng)計方法。 一般將方差分析研究的對象稱為因素，而因素中的內(nèi)容稱為水平。若方差分析同時針對兩個因素進行，則稱為雙因素方差分析。 7/27/20228版權(quán)所有 BY 統(tǒng)計學課程組一、方差分析的有關概念和基本思想 在本章案例中，縮

4、水率就是試驗指標，染整工藝是所要檢驗的因素（又稱因子），三種不同的工藝可看成是該因素的三種水平，故這是一個單因素三水平的試驗。 從表7.1可知，12個數(shù)據(jù)各不相同。一方面，同一種工藝對不同種布樣的縮水率是不同的，其差異可以看成是由于隨機因素造成的；另一方面，不同工藝對各布樣的縮水率也是不同的，這既可能是由于染整工藝類型不同造成的，也有可能是由于隨機因素造成的。 7/27/20229版權(quán)所有 BY 統(tǒng)計學課程組一、方差分析的有關概念和基本思想 要判斷隨機因素和工藝差別哪個是造成縮水率不同的主要原因，可假設三種不同的工藝為三個不同的總體，將此問題歸結(jié)為判斷三個總體是否具有相同的分布。又由于經(jīng)常假定

5、遇到的是正態(tài)總體，且在進行試驗時，除了要檢驗的因素，其他條件是盡可能保持一致的，于是可以認為每個總體的方差是相同的。這樣一來，推斷幾個總體是否具有相同分布，就可以化為檢驗幾個具有相同方差的正態(tài)總體均值是否相等的問題。這里還要假定各水平觀察值為來自正態(tài)總體的隨機樣本，各總體相互獨立且方差相同。這些假定在實際中一般難以嚴格滿足，但應對數(shù)據(jù)進行處理，使其近似地滿足正態(tài)分布。 7/27/202210版權(quán)所有 BY 統(tǒng)計學課程組二、顯著性檢驗 設檢驗的因素有m個水平，分別記為 ，在每個水平下做k次試驗，觀察值 表示第i個水平下的第j個試驗值。又設 是m個相互獨立且方差相等的正態(tài)總體， 的總體均值 ，則方

6、差分析實際上就是要檢驗假設：7/27/202211版權(quán)所有 BY 統(tǒng)計學課程組令n=mk通常稱 為組平均數(shù)，稱 為總平均數(shù)。7/27/202212版權(quán)所有 BY 統(tǒng)計學課程組平方和分解公式： 7/27/202213版權(quán)所有 BY 統(tǒng)計學課程組 稱為總離差平方和，它是描述所有數(shù)值離散程度的數(shù)量指標。 稱為組內(nèi)平方和或誤差平方和，是觀察值與組內(nèi)平均數(shù)之差的平方和，它反映了組內(nèi)（即在同一水平之下）樣本的隨機波動。 的自由度 ，其組內(nèi)方差為 。 稱為組間平方和，是組內(nèi)平均數(shù)與總平均數(shù)之差的平方和，它反映了因素水平的不同及隨機因素引起的差異。 的自由度 ，其組間方差為 。7/27/202214版權(quán)所有

7、BY 統(tǒng)計學課程組 由平方和分解公式可知， 與 的比值反映兩種差異的大小，比值越大，說明由于因素水平不同引起的差異越顯著。根據(jù)統(tǒng)計推斷的有關定理和推論，統(tǒng)計量因此，F(xiàn)作為檢驗 是否成立的檢驗統(tǒng)計量。 7/27/202215版權(quán)所有 BY 統(tǒng)計學課程組 對于給定的顯著性水平，由F分布表可查出相應自由度的臨界值 。如果， 則拒絕原假設 ，此時說明因素A對試驗指標起顯著影響；如果 ，則接受原假設 ，此時說明因素A的不同水平對試驗指標的影響不顯著。7/27/202216版權(quán)所有 BY 統(tǒng)計學課程組第二節(jié) 單因素方差分析一、各水平試驗次數(shù)相等的方差分析 試驗次數(shù)相等的單因素方差分析的具體步驟：1建立假設

8、7/27/202217版權(quán)所有 BY 統(tǒng)計學課程組2計算有關均值及平方和7/27/202218版權(quán)所有 BY 統(tǒng)計學課程組在表7.1中增加若干計算欄，計算有關均值，如表7.3所示。7/27/202219版權(quán)所有 BY 統(tǒng)計學課程組7/27/202220版權(quán)所有 BY 統(tǒng)計學課程組3列方差分析表7/27/202221版權(quán)所有 BY 統(tǒng)計學課程組對于本章案例，方差分析表如表7.5所示：7/27/202222版權(quán)所有 BY 統(tǒng)計學課程組4統(tǒng)計決策 對于顯著性水平 ，檢驗統(tǒng)計量 ，說明不同工藝方法的差異顯著。又由于 ,故第一種工藝方法（即 ）對布料縮水率的影響顯著小于其他方法，應予采用。7/27/20

9、2223版權(quán)所有 BY 統(tǒng)計學課程組二、試驗次數(shù)不等的方差分析試驗中，有時各水平下的試驗次數(shù)不相等，如表7.6所示：7/27/202224版權(quán)所有 BY 統(tǒng)計學課程組試驗次數(shù)不等的單因素方差分析的計算步驟與試驗次數(shù)相等的完全一樣，只是將K改為 即可。 例7.1 為了對幾個行業(yè)的服務質(zhì)量進行評價，某市消費者協(xié)會對該地的旅游業(yè)、居民服務業(yè)、公路客運業(yè)和保險業(yè)分別抽取了不同數(shù)量的企業(yè)。每個行業(yè)中的這些企業(yè)在服務內(nèi)容、服務對象、企業(yè)規(guī)模等方面基本相同。經(jīng)統(tǒng)計，最近一年消費者對這23家企業(yè)投訴的次數(shù)資料如表7-7所示，消費者協(xié)會想知道:這幾個行業(yè)之間的服務質(zhì)量是否有顯著差異？如果有，究竟是在哪些行業(yè)之間

10、？如果能找出哪些行業(yè)的服務質(zhì)量最差，就可以建議對消費者權(quán)益保護法中該行業(yè)的某些條款作出修正。7/27/202225版權(quán)所有 BY 統(tǒng)計學課程組7/27/202226版權(quán)所有 BY 統(tǒng)計學課程組解(1) 建立假設 不全相等 (2) 計算有關均值及平方和 7/27/202227版權(quán)所有 BY 統(tǒng)計學課程組(3) 列方差分析表(表7-8)7/27/202228版權(quán)所有 BY 統(tǒng)計學課程組(4) 統(tǒng)計決策對于顯著性水平 由于檢驗統(tǒng)計量 ，所以拒絕原假設 ，即有95%的把握認為不同的行業(yè)之間投訴的差異顯著。7/27/202229版權(quán)所有 BY 統(tǒng)計學課程組Excel中方差分析的計算步驟(1) 點擊“工具

11、”欄中的數(shù)據(jù)分析項；(2) 在分析工具框中連擊“單因素方差分析”；(3) 在對話框的“數(shù)據(jù)區(qū)域”框中鍵入A3：C6；在框中保持0.05不變（也可根據(jù)需要變?yōu)?.01）；在“輸出選項”中鍵入D3；選擇“確定”，輸出結(jié)果如表7-9所示：7/27/202230版權(quán)所有 BY 統(tǒng)計學課程組由于P=0.0207150.05，故拒絕原假設，即有95%的把握認為三種染整工藝縮水率的差異是顯著的。 7/27/202231版權(quán)所有 BY 統(tǒng)計學課程組三、方差分析中的多重比較最小顯著差異法 的基本步驟為：第一步：提出原假設： 第二步：計算各檢驗統(tǒng)計量的 值；第三步：計算LSD，其公式為：第四步：根據(jù)顯著性水平進行

12、決策：如果 的值，則拒絕 ；否則，則接受 7/27/202232版權(quán)所有 BY 統(tǒng)計學課程組對于例7.1，對四個行業(yè)的均值進行多重比較（=0.05） 由題意及計算可知， 第一步：提出假設7/27/202233版權(quán)所有 BY 統(tǒng)計學課程組第二步：計算檢驗統(tǒng)計量的值 第三步：計算LSD。 7/27/202234版權(quán)所有 BY 統(tǒng)計學課程組7/27/202235版權(quán)所有 BY 統(tǒng)計學課程組第四步：進行決策7/27/202236版權(quán)所有 BY 統(tǒng)計學課程組第三節(jié) 雙因素方差分析一、無交互作用的雙因素方差分析 若記一因素為因素A，另一因素為因素B，對A與B同時進行分析，就屬于雙因素方差分，即判斷是否有某

13、一個或兩個因素對試驗指標有顯著影響，兩個因素結(jié)合后是否有新效應。在統(tǒng)計學中將各個因素的不同水平的搭配所產(chǎn)生的新的影響稱為交互作用。我們先討論無交互作用的雙因素方差分析問題，對于有交互作用的雙因素方差分析問題稍后再討論。7/27/202237版權(quán)所有 BY 統(tǒng)計學課程組 假定因素A有r個水平： ；因素B有s個水平： 。在A的r個水平與B的s個水平的每種組合下作一次試驗，可得無交互作用的雙因素方差分析的數(shù)據(jù)結(jié)構(gòu)如表7-10所示：7/27/202238版權(quán)所有 BY 統(tǒng)計學課程組表7-10 雙因素方差分析的數(shù)據(jù)結(jié)構(gòu)7/27/202239版權(quán)所有 BY 統(tǒng)計學課程組在表7-10中：7/27/20224

14、0版權(quán)所有 BY 統(tǒng)計學課程組判斷因素A的影響是否顯著等價于檢驗假設：判斷因素B的影響是否顯著等價于檢驗假設：其中， 表示A的第i個水平所構(gòu)成的總體均值， 表示的B第j個水平所構(gòu)成的總體均值。7/27/202241版權(quán)所有 BY 統(tǒng)計學課程組 對離差總平方和進行分解。與單因素情況類似，能夠證明下列公式成立： 7/27/202242版權(quán)所有 BY 統(tǒng)計學課程組由數(shù)理統(tǒng)計可以證明：當 時 7/27/202243版權(quán)所有 BY 統(tǒng)計學課程組7/27/202244版權(quán)所有 BY 統(tǒng)計學課程組7/27/202245版權(quán)所有 BY 統(tǒng)計學課程組表7-11 雙因素（無交互作用）方差分析表 7/27/2022

15、46版權(quán)所有 BY 統(tǒng)計學課程組 例7.2 為提高某種產(chǎn)品的合格率，考察原料用量和來源地對其是否有影響。原料來源地有三個：甲、乙、丙；原料用量有三種：現(xiàn)有量、增加5%、增加8%。每個水平組合各作一次試驗，得到的數(shù)據(jù)如表7-12所示。試分析原料用量和來源地對產(chǎn)品合格率的影響是否顯著？7/27/202247版權(quán)所有 BY 統(tǒng)計學課程組表7-12 產(chǎn)品合格率數(shù)據(jù)7/27/202248版權(quán)所有 BY 統(tǒng)計學課程組解:(1) 建立假設(2) 計算相應的均值和平方和：7/27/202249版權(quán)所有 BY 統(tǒng)計學課程組(3) 列方差分析表7/27/202250版權(quán)所有 BY 統(tǒng)計學課程組(4) 統(tǒng)計決策對于

16、顯著性水平 =0.05，查表得臨界值因為 即根據(jù)現(xiàn)有數(shù)據(jù)，有95%的把握可以推斷原料來源地對產(chǎn)品合格率的影響不大，而原料用量對合格率有顯著影響。由于 為最優(yōu)水平。既然原料來源地對產(chǎn)品合格率的影響不顯著，在保證質(zhì)量的前提下，可以選擇運費最省的地方作為原料來源地選擇時的首選。如果丙地的運費最省，則最優(yōu)方案為 。 7/27/202251版權(quán)所有 BY 統(tǒng)計學課程組Excel應用以例7.2為例，設數(shù)據(jù)輸入到Excel工作表的區(qū)域為A2:C4，再按以下步驟進行：(1) 點擊“工具”欄中的數(shù)據(jù)分析項。(2) 分析工具框中連擊“方差分析：無重復的雙因素方差分析”。(3) 在對話框的“數(shù)據(jù)區(qū)域”框中鍵入B3：

17、D5；在框中保持0.05不變（也可根據(jù)需要變?yōu)?.01）；在“輸出選項”中鍵入E2；選擇“確定”，輸出結(jié)果如表7-14所示： 7/27/202252版權(quán)所有 BY 統(tǒng)計學課程組表7-14 無重復雙因素分析方差分析結(jié)果7/27/202253版權(quán)所有 BY 統(tǒng)計學課程組二、有交互作用的雙因素方差分析 所謂交互作用，簡單來說就是不同因素對試驗指標的復合作用，因素A和B的綜合效應不是二因素效應的簡單相加。為了能分辨出兩個因素的交互作用，一般每組試驗至少作兩次。7/27/202254版權(quán)所有 BY 統(tǒng)計學課程組表7-15 有交互作用的雙因素方差分析數(shù)據(jù)結(jié)構(gòu)7/27/202255版權(quán)所有 BY 統(tǒng)計學課程

18、組2建立假設7/27/202256版權(quán)所有 BY 統(tǒng)計學課程組7/27/202257版權(quán)所有 BY 統(tǒng)計學課程組對這一模型可設如下三個假設：7/27/202258版權(quán)所有 BY 統(tǒng)計學課程組3方差分析與單因素方差分析的平方和分解類似，有 7/27/202259版權(quán)所有 BY 統(tǒng)計學課程組7/27/202260版權(quán)所有 BY 統(tǒng)計學課程組7/27/202261版權(quán)所有 BY 統(tǒng)計學課程組表7.16 雙因素（有交互作用）方差分析表7/27/202262版權(quán)所有 BY 統(tǒng)計學課程組 例7.3 某公司想將橡膠、塑料和軟木的板材沖壓成密封墊片出售。市場上有兩種不同型號的沖壓機可供選擇。為了能對沖壓機每小

19、時所生產(chǎn)的墊片數(shù)進行比較，并確定哪種機器使用何種材料生產(chǎn)墊片的能力更強，該公司使用每臺機器對每一種材料分別運行三段時間，得到的試驗數(shù)據(jù)（每小時生產(chǎn)的墊片數(shù)）如表7-16所示，試運用方差分析確定最優(yōu)方案。7/27/202263版權(quán)所有 BY 統(tǒng)計學課程組7/27/202264版權(quán)所有 BY 統(tǒng)計學課程組解(1) 建立假設： (2) 計算相應的均值和平方和：7/27/202265版權(quán)所有 BY 統(tǒng)計學課程組7/27/202266版權(quán)所有 BY 統(tǒng)計學課程組(3) 列方差分析表（7-19）7/27/202267版權(quán)所有 BY 統(tǒng)計學課程組(4) 統(tǒng)計決策由于 ,說明不僅沖壓機的型號和墊片材料對墊片數(shù)

20、量有顯著影響，而且其交互作用也是顯著的。由結(jié)構(gòu)均值表可知，在沖壓機中,第一種的均值較大；墊片材料中,木的均值較大，故最優(yōu)方案是 。 7/27/202268版權(quán)所有 BY 統(tǒng)計學課程組有交互作用的雙因素方差分析中的Excel應用 以例7.3為例，將數(shù)據(jù)輸入到Excel工作表的區(qū)域為B2:D7，沖壓機的二個水平A1、A2在第一列，墊片材料的水平B1、B2、B3在第一行。再按以下步驟進行：(1) 點擊“工具”欄中的數(shù)據(jù)分析項。(2) 分析工具框中連擊“方差分析：可重復雙因素方差分析”。(3) 在對話框的“數(shù)據(jù)區(qū)域”框中鍵入A1:D7；在“每一樣本的行數(shù)” 框中鍵入3。在框中保持0.05不變（也可根據(jù)

21、需要變?yōu)?.01）；在“輸出選項”中鍵入空白地方，如F2；選擇“確定”，輸出結(jié)果如表7-20所示： 7/27/202269版權(quán)所有 BY 統(tǒng)計學課程組7/27/202270版權(quán)所有 BY 統(tǒng)計學課程組7/27/202271版權(quán)所有 BY 統(tǒng)計學課程組第四節(jié) 正交試驗設計初步一、試驗設計的概念和設計原則試驗設計有三個基本原則：重復性、隨機化和區(qū)組化。(1) 重復性原則：重復性是指對一項試驗要在相同的條件下重復進行若干次。只有進行多次的試驗，才能掌握其規(guī)律性。(2) 隨機化原則：隨機化是指試驗材料和試驗地點都要隨機地確定。這樣進行試驗得出的結(jié)論才具有客觀性和普遍性，且每次進行的試驗都可認為是相互獨

22、立的。(3) 區(qū)組化原則：一組試驗，試驗者總希望在相同或近似相同的條件下進行，以便在相互比較中得出正確的結(jié)論。 7/27/202272版權(quán)所有 BY 統(tǒng)計學課程組在安排試驗時需要注意三點： 第一、盡量減少試驗誤差。在試驗時，要盡可能使對試驗產(chǎn)生影響的其他因素達到理論中要求的精確程度，這樣我們在對目的指標進行測算時才能達到較為客觀的結(jié)果。在一項試驗中誤差是必然存在的，我們的目的就是要盡量減少該誤差。第二、盡量減少試驗次數(shù)。試驗的次數(shù)越多，進行試驗所耗費的人力、物力也就會相應地增多，最佳的試驗就是用盡可能少的費用來獲得最有效的試驗結(jié)果。第三、所設計的試驗要便于對指標值進行統(tǒng)計分析。影響一個事件的指

23、標是多樣的，在進行設計時要選擇比較容易測算且對事件影響較大的指標進行計算。7/27/202273版權(quán)所有 BY 統(tǒng)計學課程組 在現(xiàn)實問題中，影響試驗指標的因素通常有很多個，要考察它們就涉及多因子的試驗設計問題。多因子試驗中一個很困難的問題就是因子數(shù)較多，而又需要進行多次的試驗。例如，有10個因子對某一指標有顯著影響，而每個因子取兩個水平進行比較，那么就有210=1024個不同的水平組合，即每個水平組合作一次試驗，就需要進行1024次試驗，耗費的人力、財力可想而知，這在實際中是不可行的。因此，需要按照一定的方法從中選擇一部分進行試驗。常用的試驗設計方法有正交試驗設計法、參數(shù)設計法、回歸設計法、均

24、勻設計法、混料設計法等。限于篇幅，這里僅介紹正交試驗設計的基礎知識，其他方法請讀者自行閱讀試驗設計的有關書籍。7/27/202274版權(quán)所有 BY 統(tǒng)計學課程組二、多因子試驗問題 在多因子試驗中，各因子又有不同的水平數(shù)，我們的目的是要從這些因子不同的水平組合中，找出一組或幾組組合使所要求的指標達到最優(yōu)。下面以一個二因子的例子來具體認識一下多因子試驗問題。例7.4 為提高合金鋼的強度，同時考慮碳(C)含量(因子)及鈦(Ti)與鋁(AL)的含量(因子)對強度Y的影響，希望找出最佳的含量組合，使強度Y達到最大（表7-21）。7/27/202275版權(quán)所有 BY 統(tǒng)計學課程組表7-21 合金鋼的試驗數(shù)

25、據(jù)7/27/202276版權(quán)所有 BY 統(tǒng)計學課程組7/27/202277版權(quán)所有 BY 統(tǒng)計學課程組三、交互作用 一個因子的水平好壞或好壞的程度受另一個因子水平制約的情況，稱為因子A和B的交互作用，記作 或AB。 因子A和B的交互作用可以用圖形較為直觀地表示。 7/27/202278版權(quán)所有 BY 統(tǒng)計學課程組7/27/202279版權(quán)所有 BY 統(tǒng)計學課程組 因子間的交互作用會隨著因子個數(shù)的增加而增加。如四個因子A，B，C，D間的交互作用有以下幾類：(1) 二級交互作用有6個： AB， AC， AD， BC， BD， CD；(2) 三級交互作用有4個： ABC， ABD， ACD， BCD

26、；(3) 四級交互作用有1個： ABCD 。 交互作用共有11個，比因子個數(shù)還多。實踐經(jīng)驗表明，多數(shù)交互作用是不存在或者很小以至可以忽略不計的，一般我們主要考慮部分二級交互作用，但具體考察哪些二級交互作用還要根據(jù)下面的方法來決定。 7/27/202280版權(quán)所有 BY 統(tǒng)計學課程組四、正交表及其類型 正交試驗設計法是利用正交表科學地挑選試驗條件，合理安排試驗的方法，是研究與處理多因素試驗的一種科學方法。 正交設計有兩個重要的特點：（1）任意一對因子（也稱因素）的任一水平組合必在試驗中出現(xiàn)，且出現(xiàn)的次數(shù)相同；（2）總試驗次數(shù)比全面試驗（所有因子的任一水平組合都進行搭配）要少許多次。例如，對7個二

27、水平因子進行全面試驗要進行27=128次，而用正交表安排試驗只需要作8次。用正交表合理地安排試驗，可以做到省時、省力、省錢，還能得到令人滿意的檢驗效果，因此這種方法在改進產(chǎn)品質(zhì)量、研究采用新工藝、試制新產(chǎn)品、了解設備工藝性能以及改進技術管理等方面都有廣泛的應用。 7/27/202281版權(quán)所有 BY 統(tǒng)計學課程組1正交表及其特性 正交表是正交設計的工具，是運用組合數(shù)學理論在正交拉丁方的基礎上構(gòu)造的一種規(guī)格化表格，符號為：其中，L為正交表符號； n為正交表的行數(shù)(試驗次數(shù)，試驗方案數(shù))； j為正交表中每一列因子的水平個數(shù)； i為正交表的列數(shù)(試驗因子的個數(shù))。 7/27/202282版權(quán)所有 B

28、Y 統(tǒng)計學課程組 例如， 表示一張8行7列，每列有兩個水平的正交表（表7-22），它表示在這個試驗中，要作8次不同條件的試驗，共有7個影響因素，每個因素都有兩個水平，我們稱這張表為兩水平的正交表。7/27/202283版權(quán)所有 BY 統(tǒng)計學課程組7/27/202284版權(quán)所有 BY 統(tǒng)計學課程組 正交表具有正交性，這是指它有如下兩個特征：（1）每列中不同的數(shù)字重復的次數(shù)相同。在表7.22中，每列有兩個不同的數(shù)字：1，2，每一個各出現(xiàn)4次。（2）將任意兩列的同行數(shù)字看成一個數(shù)對，那么任意可能數(shù)對重復的次數(shù)相等。在 中，任意兩列有4種可能數(shù)對：(1，1)， (1，2)，(2，1)，(2，2)，每一

29、對各出現(xiàn)2次。 7/27/202285版權(quán)所有 BY 統(tǒng)計學課程組2正交表的分類7/27/202286版權(quán)所有 BY 統(tǒng)計學課程組五、無交互作用情況下的設計正交表安排試驗的步驟如下：（1） 明確試驗目的，確定要考察的試驗指標。（2） 確定要考察的因子和因子的水平。（3） 選用合格的正交表，進行表頭設計。（4） 根據(jù)試驗號的安排進行試驗，并記錄試驗指標的具體數(shù)據(jù)。（5） 數(shù)據(jù)分析。對一個正交表形式的試驗設計通常有三種分析方法：一是用極差分析各因子對指標影響程度的大小，這是一種較為簡單的直觀分析方法；二是用方差分析進行數(shù)據(jù)分析；三是貢獻率分析法。 7/27/202287版權(quán)所有 BY 統(tǒng)計學課程組

30、例7.5 某化工廠生產(chǎn)的一種產(chǎn)品的收率較低，為此希望通過試驗提高收率。在試驗中考察如下三個因子三個水平（表7-23）：表7-23 因子水平表 7/27/202288版權(quán)所有 BY 統(tǒng)計學課程組解 由于所考察的因子是三水平的，因此選用三水平正交表，又現(xiàn)在只考察三個因子，故選擇 安排試驗。 選定了正交表后把因子放在正交表的列上去，稱為表頭設計。在例7.5中將三個因子置于前三列，將它寫成如下的表頭設計形式： 7/27/202289版權(quán)所有 BY 統(tǒng)計學課程組九次試驗的結(jié)果收率(%)分別是：51 61 58 72 69 59 87 85 84(一) 用極差分析各因子對指標影響程度的大小1用直觀分析表對

31、該試驗進行分析為方便起見，把試驗結(jié)果寫在正交表的右邊一列上，并分別用 表示，所有計算可以在表上進行。7/27/202290版權(quán)所有 BY 統(tǒng)計學課程組 首先來看第一列，該列中的1，2，3分別表示因子的三個水平，按水平號將數(shù)據(jù)分為三組：“1”對應 ，“2”對應 ，“3”對應 。 在第二列，該列中的1，2，3分別表示因子B的三個水平，按水平號將數(shù)據(jù)分為三組：“1”對應 ，“2”對應 ，“3”對應 。同理我們可以從圖中得到第三列的情況。7/27/202291版權(quán)所有 BY 統(tǒng)計學課程組 分別對每一列每個水平的數(shù)據(jù)進行相加，得到這個水平下該因子對應指標值的和，例如，7/27/202292版權(quán)所有 BY

32、 統(tǒng)計學課程組 同理我們可以得到對應C水平的T1，T2，T3的值，用每個水平得出的Ti的值去除以3就得到其均值 。 由以上計算可知，每個因子下間的差異反映了這個因子三個水平間的差異；分別對三個因子的各個水平進行數(shù)值大小比較：因子A的三個水平均值差異較大，其第三個水平的均值最大，故因子A的三水平最好。對第二、三列進行類似地分析，可知因子B的第二個水平好，因子C的第二個水平好。 7/27/202293版權(quán)所有 BY 統(tǒng)計學課程組對第四列也可以進行上述分析，按其中的1，2，3分別將數(shù)據(jù)分為三組，但三組的水平組合相同，因此該列僅反映誤差。綜上可知，使指標達到最佳的水平組合是，即反應溫度為90度、加堿量

33、為48公斤、選取乙類催化劑可以使轉(zhuǎn)化率達到最大。7/27/202294版權(quán)所有 BY 統(tǒng)計學課程組2用極差分析各因子對指標影響程度的大小 一個因子的極差是該因子各水平均值的最大值與最小值之差，如果該值大，則說明改變這一因子的水平會對指標造成較大的變化，所以該因子對指標的影響大，反之，則影響小。極差最大的列所對應的因子是最主要因子。當要求指標愈高愈好時，可選R行中最大者，所對應的水平為優(yōu)水平；當要求指標愈低愈好時，可選R行中最小者其所對應的水平為優(yōu)水平。 7/27/202295版權(quán)所有 BY 統(tǒng)計學課程組在例7.5中各個因子的極差分別是：它們被置于表7-24的最下方一行。從三個因子的極差可知因子

34、A的影響最大，其次是C，而因子B的影響則最小，通常記為 7/27/202296版權(quán)所有 BY 統(tǒng)計學課程組表7-24 例7.5的直觀分析計算表7/27/202297版權(quán)所有 BY 統(tǒng)計學課程組(二) 用方差分析法進行數(shù)據(jù)分析 用方差分析方法來說明影響指標的因子的分析步驟如下：1建立統(tǒng)計模型要對試驗進行方差分析，首先要作幾點假設：(1) 在同一水平組合下全部試驗結(jié)果構(gòu)成一個總體，且服從正態(tài)分布；(2) 各正態(tài)總體的方差相同，均為 ；(3) 各正態(tài)均值與水平組合有關；(4) 不同水平組合下的試驗是相互獨立的。 7/27/202298版權(quán)所有 BY 統(tǒng)計學課程組7/27/202299版權(quán)所有 BY

35、統(tǒng)計學課程組2進行平方和分解總平方和為：7/27/2022100版權(quán)所有 BY 統(tǒng)計學課程組3計算各因子的平方和（1）因子的平方和是此因子各水平均值與總均值間的離差平方和；（2）一個置于正交表第k列上的因子，其平方和 ，其自由度 ；（3） 因子平方和中除了誤差之外只反映此因子的效應間的差異；（4）誤差平方和 等于諸空白列的平方和之和，其自由度等于諸空白列的自由度之和。 7/27/2022101版權(quán)所有 BY 統(tǒng)計學課程組4方差分析表7/27/2022102版權(quán)所有 BY 統(tǒng)計學課程組7/27/2022103版權(quán)所有 BY 統(tǒng)計學課程組5最佳水平選擇(1)對顯著因子應選擇其最好的水平；(2)對不

36、顯著的因子可以任意選擇水平，?？筛鶕?jù)降低成本、操作方便等來選擇其水平。7/27/2022104版權(quán)所有 BY 統(tǒng)計學課程組6最佳水平組合均值估計7/27/2022105版權(quán)所有 BY 統(tǒng)計學課程組7/27/2022106版權(quán)所有 BY 統(tǒng)計學課程組7驗證試驗 在實際問題中,分析所得的最佳水平組合無論是否在試驗中出現(xiàn)，通常都需要作驗證試驗，對最佳水平組合進行驗證。 對于例7.5，根據(jù)介紹的計算公式我們得到如表7-26所示的方差分析表：7/27/2022107版權(quán)所有 BY 統(tǒng)計學課程組7/27/2022108版權(quán)所有 BY 統(tǒng)計學課程組(三) 貢獻率分析 當試驗指標不服從正態(tài)分布時，進行方差分析

37、的依據(jù)就不充足，此時可以通過比較各因子的“貢獻率”來衡量因子作用的大小。7/27/2022109版權(quán)所有 BY 統(tǒng)計學課程組表7-27 貢獻率分析表 7/27/2022110版權(quán)所有 BY 統(tǒng)計學課程組 貢獻率最大的幾個因子是重要因子，而貢獻率與誤差貢獻率差不多的因子通常被認為是不重要的。 將例7.5中的數(shù)據(jù)作貢獻率分析,計算結(jié)果如表7-28所示：7/27/2022111版權(quán)所有 BY 統(tǒng)計學課程組六、有交互作用設計 在多因子試驗中，經(jīng)常會出現(xiàn)兩個或多個因子間有交互作用的情況，此時，要用關于有交互作用的方法進行分析。因子間有交互作用與無交互作用的分析方法是有區(qū)別的，這主要體現(xiàn)在正交表的設計上。7/27/2022112版權(quán)所有 BY 統(tǒng)計學課程組(一) 表頭設計1自由度7/27/2022113版權(quán)所有 BY 統(tǒng)計學課程組2在對有交互作用的試驗進行分析時應遵循的原則7/27/2022114版權(quán)所有 BY 統(tǒng)計學課程組(二) 數(shù)據(jù)的方差分析 以下的分析僅以4個因子及交互作用的二水平正交設計為例，描述統(tǒng)計模型和