經(jīng)濟(jì)學(xué)第七章方差分析與正交試驗(yàn)設(shè)計(jì)初步課件

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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