第六講方差分析與實驗設(shè)計

第六講方差分析與實驗設(shè)計第一頁，共六十四頁，編輯于2023年，星期一本講主題方差分析－－多個總體平均值相等的檢驗Testingfortheequalityofkpopulationmeans多重比較過程Multiplecomparisonprocedures完全隨機設(shè)計Completelyrandomizeddesigns隨機區(qū)組設(shè)計Randomizedblockdesign多因素實驗Factorialexperiments第二頁，共六十四頁，編輯于2023年，星期一方差分析－多個等方差正態(tài)總體均值的假設(shè)檢驗引例：為了了解工人對綜合質(zhì)量管理了解多少，從三個工廠各隨機抽取6名工人進行測驗，得分如下觀測個體工廠1考分工廠2考分工廠3考分185715927575643827362476746957169756858267樣本均值797466樣本方差342032第三頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）現(xiàn)在的問題是檢驗假設(shè)不是所有總體的均值都相等在之下，三個樣本均來自同一總體，由樣本均值的抽樣分布可知，三個樣本均值的均值也是總體均值的最好估計值樣本均值的方差的估計值為

第四頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）

由此得總體方差的樣本間估計值為這個值遠大于總體方差用各個樣本方差作估計的估計值34、20、32，也遠大于各個樣本方差的平均值（或稱為的樣本內(nèi)估計值）

第五頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）

如果前面的假設(shè)正確，基于樣本均值本身之間的差異的的估計值與基于每個樣本內(nèi)數(shù)據(jù)的的估計值應(yīng)該不會有太大的差異，即二者之間的比值應(yīng)接近1才對。如果差異太大，即二者比值遠大于1，說明三個樣本并非來自于一個總體，即前面假設(shè)就不真，那么究竟怎樣才算差異大呢，也就是二者的比例什么情況下才算遠大于1呢？第六頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）考慮m個等方差的正態(tài)總體要檢驗不是所有總體的方差都相等

為此，從第i（）個總體中隨機抽取容量為的樣本，得觀測值如下表第七頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）

觀測個體總體1…i…總體m1……2…………………………………………………………樣本均值樣本方差第八頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）令為所有觀測值的平均值，其中。當時第九頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）兩種變差隨機變差：同一工廠里的工人的成績是不同的。系統(tǒng)變差：不同的工廠里的工人的成績是不同的。

要判別不同的工廠是否是造成工人成績差異的主要因素，歸結(jié)為判斷三個總體是否具有相同分布的問題。第十頁，共六十四頁，編輯于2023年，星期一則總體誤差平方和SST被分解為兩個部分：組間誤差平方SSB和及組內(nèi)誤差平方和SSW多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）第十一頁，共六十四頁，編輯于2023年，星期一記則在之下多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）第十二頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）則對任意給定的，可找到，使利用樣本觀測值求出組間方差估計值與組內(nèi)方差估計之比，若，則拒絕接受原假設(shè)，否則接受原假設(shè)。

第十三頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）

方差分析表方差來源平方和自由度均方和F－值組間SSBM-1MSBMSB/MSW組內(nèi)SSWN-mMSW總計SSTN-1MST第十四頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－方差分析（等方差情況）

方差分析表方差來源平方和自由度均方和F－值組間51622589組內(nèi)4301528.67F0.05(2,15)=3.68總計94617第十五頁，共六十四頁，編輯于2023年，星期一

對任意，的關(guān)于的否定區(qū)域為多個正態(tài)總體均值的假設(shè)檢驗

－多重比較（Fisher’sLSD）第十六頁，共六十四頁，編輯于2023年，星期一在前例中

因此沒有理由拒絕假設(shè)，而有理由拒絕假設(shè)和多個正態(tài)總體均值的假設(shè)檢驗

－多重比較（Fisher’sLSD）第十七頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗

－多重比較（Fisher’sLSD）一、方差具有齊性時均值的比較方法有多種，如：1、LSD：用T檢驗完成各組均值間的配對比較，對多重比較誤差不進行調(diào)整；2、LSDMOD：用T檢驗完成各組間均值的配對比較，但通過設(shè)置每個檢驗的誤差率來控制整個誤差率二、方差不具齊性時均值的比較方法也有多種，如1、Tamhane’sT2：檢驗進行配對比較2：Dunnett’sT3：正態(tài)分布下的配對比較第十八頁，共六十四頁，編輯于2023年，星期一多個正態(tài)總體均值的假設(shè)檢驗SS----SumofSquare（平方和）df----degreeoffreedom(自由度）MS----MeanSquare,Variance(方差）F----由樣本算出的F值P-value----拒絕H0犯錯的概率F-crit----在所要求的顯著水平下的F的臨界值第十九頁，共六十四頁，編輯于2023年，星期一多個等方差正態(tài)總體均值的假設(shè)檢驗

－方差分析方差分析的SPSS實現(xiàn)1、打開SPSS數(shù)據(jù)編輯器（SPSSDataEditor）2、將數(shù)據(jù)按因變量及自變量形成具有兩個變量的數(shù)據(jù)文件，其中自變量即是因素的各不同水平，均用數(shù)字代表，因變量即是各不同水平對應(yīng)的觀測值3、打開主菜單Analyze中的“CompareMean”選項，打開“One-WayANOVA”對話窗，將因變量移入“DependentList”表框，將自變量移入“Factor”表框，點擊“OK”，基本的方差分析便可完成。第二十頁，共六十四頁，編輯于2023年，星期一多個等方差正態(tài)總體均值的假設(shè)檢驗

－方差分析方差分析的SPSS實現(xiàn)1、Polynomial選框：選擇此項，后面的“degree”變?yōu)榭捎?，系統(tǒng)默認為“Linear”。即自、因變量間的線性關(guān)系檢驗。2、Coefficients窗口：有時候需要比較自變量不同取值對應(yīng)的分組數(shù)據(jù)的均值，例如某一組的均值是否與另一組或另幾組的均值相等，或某組均值的多少倍是否與另一組或幾組的均值的多少倍相等，等等。3、Previous、Next：閱讀、修改、輸入前后面的系數(shù)第二十一頁，共六十四頁，編輯于2023年，星期一1）方差分析是多個總體進行的假設(shè)檢驗，必須在前述兩個基本假定前提下進行，即各總體方差必須相等。2）利用F檢驗時，若組內(nèi)方差的自由度(r-1)n很小，則F檢驗靈敏度很低。此時，必須增加試驗次數(shù)或放寬，以達到實踐中檢驗的可靠性。3）在總體不遵從正態(tài)分布時，可采用非參數(shù)分析法或?qū)⒆兞孔鬟m當變換，使變換后的隨機變量遵從正態(tài)分布且方差相等，然后再進行方差分析。第二十二頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－基本術(shù)語實驗：收集樣本數(shù)據(jù)的過程實驗設(shè)計：收集樣本數(shù)據(jù)的計劃響應(yīng)變量：實驗中測量的變量實驗單元：測量響應(yīng)變量的對象因素：凡與響應(yīng)變量有關(guān)系的變量稱為因素因素水平：因素的不同等級處理：實驗中若干個因素的若干水平的特定組合

重復：同一水平或處理的多次測量第二十三頁，共六十四頁，編輯于2023年，星期一例如：一項市場營銷研究的內(nèi)容是考察品牌和貨架位置對咖啡周銷量的影響。分別記錄放在三個貨架位置（底部、頂部和中部）的兩種品牌（品牌A和品牌B）的咖啡銷售量。品牌與貨架的不同組合共有6種，在18周的時間內(nèi)，每周變換一種組合。這種設(shè)計的安排如下表.試識別本實驗單元。貨架位置底中頂品牌A

B第(1,9,14)周第(2,7,16)周第(4,12,17)周第(5,10,13)周第(3,8,18)周第(6,11,14)周實驗設(shè)計－基本術(shù)語第二十四頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－完全隨機化設(shè)計完全隨機化設(shè)計是一種將K種處理隨機指派給各個實驗單元的設(shè)計，或者說是一種從K個總體分別抽取獨立隨機樣本的設(shè)計。例：一公司考慮用三種方法為銷售人員支付酬金：只給傭金，給固定薪金，較低的固定薪金加酬金。為比較哪種方法能調(diào)動銷售人員的極積性，分別從三類銷售人員中隨機抽取7名，檢查其月銷售額，列入下表第二十五頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－完全隨機化設(shè)計傭金固定薪金傭金加固定金165120140981151561309022021012611219510713418715523524080待查175113.29166.17第二十六頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－完全化隨機設(shè)計觀測個體處理1…處理i…處理k個體1……個體2……………………………………………個體……處理均值第二十七頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－完全化隨機設(shè)計在完全化隨機實驗設(shè)計中，總體誤差平方和SST被分解為兩個部分：處理間誤差平方SSTR和及隨機誤差平方和SSE第二十八頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－完全化隨機設(shè)計記則在之下第二十九頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－完全隨機化設(shè)計

方差分析表方差來源平方和自由度均方和F－值組間SSTRM-1MSTRMST/MSE組內(nèi)SSEN-mMSE總計SSTN-1MST第三十頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－隨機化區(qū)組設(shè)計隨機化區(qū)組設(shè)計是一種在M個區(qū)組的每一區(qū)組內(nèi)對K種處理進行比較的設(shè)計。每個區(qū)組含K個匹配的實驗單元和K種隨機分派的處理，并要求將每種處理分派給每個區(qū)組內(nèi)的每一個單元。例：某超市A有一天將3個競爭超市（超市B、超市C和超市D）的49種雜貨項目的售價與本超市售價作了比較。下表列出隨機選出的7個項目的價格。試問：怎樣才能確定這四個超市雜貨項目的平均價格之間是否有差別？第三十一頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－隨機化區(qū)組設(shè)計

雜貨項目超市A超市B超市C超市D麥片1.1美元1.18美元1.39美元1.18美元果凍0.240.240.310.26肥皂0.520.600.630.55克利斯科油1.261.702.271.29潔唇紙0.670.700.790.70金槍魚肉0.630.660.790.63德蒙特豌豆0.430.470.650.47第三十二頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－隨機區(qū)組設(shè)計

觀測個體處理1…處理i…處理k區(qū)組均值區(qū)組1……區(qū)組2……………………………………………………區(qū)組b……處理均值第三十三頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－隨機區(qū)組設(shè)計在隨機區(qū)組實驗設(shè)計中，總體差異平方和被分解為三部分：處理平方和、區(qū)組平方和及誤差平方和：

第三十四頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－隨機區(qū)組設(shè)計令則在之下

第三十五頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－隨機化區(qū)組設(shè)計

方差分析表方差來源平方和自由度均方和F－值處理間SSTRk-1MSTRMST/MSE區(qū)組間SSBLb-1MSBL誤差SSE(k-1)(b-1)MSE總計SSTN-1MST第三十六頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－隨機區(qū)組設(shè)計在前例中

第三十七頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－隨機區(qū)組設(shè)計令則在之下

第三十八頁，共六十四頁，編輯于2023年，星期一實驗設(shè)計－隨機化區(qū)組設(shè)計

方差分析表方差來源平方和自由度均方和F－值處理間SSTRk-1MSTRMST/MSE區(qū)組間SSBLb-1MSBL誤差SSE(k-1)(b-1)MSE總計SSTN-1MST第三十九頁，共六十四頁，編輯于2023年，星期一多因素實驗設(shè)計影響試驗結(jié)果（因變量的取值）的因素同時有多個。主效應(yīng)：因素水平的改變所造成的試驗結(jié)果的改變；交互效應(yīng)：當某一因素的效應(yīng)隨另一因素的水平不同而不同，則稱職這兩個因素之間存在交互作用。由于交互作用引起的試驗結(jié)果的改變稱為交互效應(yīng)。因素間是否存在交互作用，直觀上可通過交互作用圖判斷，理論上可通過專門的數(shù)學方法進行檢驗。第四十頁，共六十四頁，編輯于2023年，星期一因素B因素AA1A2…ArB1x11x12…x1rX1.B2x21x22…x2rX2.………………Bsxs1xs2…xsrXs.X.1X.2…X.r特點：每一處理只有一個觀測值雙因素實驗設(shè)計—不考慮交互作用第四十一頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計－不考慮交互作用例如：考慮不同推進器及不同燃料對射程的影響，試驗結(jié)果如下：推進器1推進器2推進器3燃料158.256.265.3燃料249.154.151.6燃料360.170.939.2燃料475.858.248.7第四十二頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計—不考慮交互作用不同推進器及不同燃料對射程的影響圖第四十三頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計—不考慮交互作用方差分析模型：觀察值=未知常數(shù)+因素A效應(yīng)+因素B效應(yīng)+殘差因素A與因素B間不存在相互作用某一特定的組合（單元）的聯(lián)合效應(yīng)不會大于各因素的主效應(yīng)之和。第四十四頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計—不考慮交互作用離差平方和分解其中：

——總誤差平方和

——

因素A的誤差平方和

——因素B的誤差平方和

——

隨機誤差平方和第四十五頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計—不考慮交互作用假設(shè)檢驗：判斷因素A的影響是否顯著等價于檢驗假設(shè)H01：.1=.2=…=.r(A因素各水平均值相等)；H11:

至少有一個等式不成立統(tǒng)計量判斷因素B的影響是否顯著等價于檢驗假設(shè)H02：1.=2.=…=s.(B因素各水平均值相等)；H12:

至少有一個等式不成立統(tǒng)計量如果F>F，則拒絕H0。第四十六頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計—不考慮交互作用SPSS實現(xiàn)（1）

1、在SPSSDataEditor中將數(shù)據(jù)依左圖形式編輯；2、打開主菜單Analyze中的“GeneralLineralModel”選項，打開“Univariate”對話窗，將因變量移入“DependentList”表框，將自變量移入“FixedFactor”表框；3、打開“Model”對話窗第四十七頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計—不考慮交互作用SPSS實現(xiàn)（2）1、點黑“Custom”2、在“BuildTerm”欄，選擇“MainEffects”；3、將“Factor”欄內(nèi)的自變量移入“model”欄；4、點擊“Continue”5、基本設(shè)置完畢，點擊“Ok”。第四十八頁，共六十四頁，編輯于2023年，星期一SPSS輸出結(jié)果結(jié)果表明：采用不同的推進器及不同的燃料對火箭的射程均沒有顯著影響。雙因素實驗設(shè)計—不考慮交互作用第四十九頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計—考慮交互作用

因素B總和B1B2Bb

因素AA1x111,…,x11rx211,…,x21r…xb11,…,xb1r

A2x121,…,x12rx221,…,x22r…xb21,…,xb2r

……………

Aax1a1,…,x1arx2a1,…,x2ar…xba1,…,xbar

總和

特點：各處理有重復值第五十頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計—考慮交互作用例如：考虛4種不同燃料與3種不同型號的推進器的搭配對火箭射程（單位：海里）的影響，每種搭配發(fā)射了2次火箭，共做了24次試驗，得數(shù)據(jù)如下：

要求分析燃料和推進器及搭配的不同是否對火箭的射程有顯著影響。零假設(shè)是沒有影響。推進器1推進器2推進器3燃料158.2,52.656.2,41.265.3,60.8燃料249.1,42.854.1,50.551.6,48.4燃料360.1,58.370.9,73.239.2,40.7燃料475.8,71.558.2,51.048.7,41.4第五十一頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計—考慮交互作用方差分析模型

i=1,2,…,a;j=1,2,…,a;k=1,2,…,r觀察值=未知常數(shù)+因素A效應(yīng)+因素B效應(yīng)+互作效應(yīng)+殘差第五十二頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計－考慮交互作用平方各分解：在考慮交互作用存在的雙因素（因素A和因素B）實驗設(shè)計中，總體差異平方和SST被分解為四部分：A因素平方和SSA、B因素平方和SSB、互作平方和SSAB及誤差平方和SSE：

第五十三頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計－考慮交互作用其中第五十四頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計－考慮交互作用令

第五十五頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計－考慮交互作用假設(shè)檢驗：判斷因素A的影響是否顯著等價于檢驗

H01：α1=α2=…=αr=0

(A因素各水平均值相等)；或：H11:

至少有一個等式不成立統(tǒng)計量如果F>F，則拒絕H0。第五十六頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計－考慮交互作用假設(shè)檢驗：判斷因素B的影響是否顯著等價于檢驗

H02：β1=β2=…=βs=0（B因素各水平均值相等)；或：H12:

至少有一個等式不成立統(tǒng)計量

如果F>F，則拒絕H0。第五十七頁，共六十四頁，編輯于2023年，星期一雙因素實驗設(shè)計－考慮交互作用假設(shè)檢驗判斷因素A、B的交互影響是否顯著等價于檢驗

H03：（αβ）ij=0

(互作因素各水平均值