單因素方差分析與多重比較

單因素方差分析與多重比較第一頁，共五十五頁，編輯于2023年，星期五§3.1方差分析引論方差分析及其有關(guān)術(shù)語方差分析的基本思想和原理方差分析的基本假定問題的一般提法第二頁，共五十五頁，編輯于2023年，星期五什么是方差分析(ANOVA)?

(analysisofvariance)檢驗多個總體均值是否相等通過分析觀察數(shù)據(jù)的誤差判斷各總體均值是否相等研究分類型自變量對數(shù)值型因變量的影響一個或多個分類尺度的自變量2個或多個(k個)處理水平或分類一個間隔或比率尺度的因變量有單因素方差分析和雙因素方差分析單因素方差分析：涉及一個分類的自變量雙因素方差分析：涉及兩個分類的自變量第三頁，共五十五頁，編輯于2023年，星期五什么是方差分析?

(例題分析)消費者對四個行業(yè)的投訴次數(shù)

行業(yè)觀測值零售業(yè)旅游業(yè)航空公司家電制造業(yè)12345675766494034534468392945565131492134404451657758【例】為了對幾個行業(yè)的服務(wù)質(zhì)量進行評價，消費者協(xié)會在四個行業(yè)分別抽取了不同的企業(yè)作為樣本。最近一年中消費者對總共23家企業(yè)投訴的次數(shù)如下表第四頁，共五十五頁，編輯于2023年，星期五什么是方差分析?

(例題分析)分析四個行業(yè)之間的服務(wù)質(zhì)量是否有顯著差異，也就是要判斷“行業(yè)”對“投訴次數(shù)”是否有顯著影響作出這種判斷最終被歸結(jié)為檢驗這四個行業(yè)被投訴次數(shù)的均值是否相等如果它們的均值相等，就意味著“行業(yè)”對投訴次數(shù)是沒有影響的，即它們之間的服務(wù)質(zhì)量沒有顯著差異；如果均值不全相等，則意味著“行業(yè)”對投訴次數(shù)是有影響的，它們之間的服務(wù)質(zhì)量有顯著差異第五頁，共五十五頁，編輯于2023年，星期五方差分析中的有關(guān)術(shù)語因素或因子(factor)所要檢驗的對象要分析行業(yè)對投訴次數(shù)是否有影響，行業(yè)是要檢驗的因素或因子水平或處理(treatment)因子的不同表現(xiàn)零售業(yè)、旅游業(yè)、航空公司、家電制造業(yè)就是因子的水平觀察值在每個因素水平下得到的樣本值每個行業(yè)被投訴的次數(shù)就是觀察值第六頁，共五十五頁，編輯于2023年，星期五方差分析中的有關(guān)術(shù)語試驗這里只涉及一個因素，因此稱為單因素四水平的試驗總體因素的每一個水平可以看作是一個總體比如零售業(yè)、旅游業(yè)、航空公司、家電制造業(yè)可以看作是四個總體樣本數(shù)據(jù)被投訴次數(shù)可以看作是從這四個總體中抽取的樣本數(shù)據(jù)第七頁，共五十五頁，編輯于2023年，星期五方差分析的基本思想和原理

(圖形分析)零售業(yè)旅游業(yè)航空公司家電制造第八頁，共五十五頁，編輯于2023年，星期五從散點圖上可以看出不同行業(yè)被投訴的次數(shù)是有明顯差異的即使是在同一個行業(yè)，不同企業(yè)被投訴的次數(shù)也明顯不同家電制造也被投訴的次數(shù)較高，航空公司被投訴的次數(shù)較低行業(yè)與被投訴次數(shù)之間有一定的關(guān)系如果行業(yè)與被投訴次數(shù)之間沒有關(guān)系，那么它們被投訴的次數(shù)應(yīng)該差不多相同，在散點圖上所呈現(xiàn)的模式也就應(yīng)該很接近方差分析的基本思想和原理

(圖形分析)第九頁，共五十五頁，編輯于2023年，星期五僅從散點圖上觀察還不能提供充分的證據(jù)證明不同行業(yè)被投訴的次數(shù)之間有顯著差異這種差異也可能是由于抽樣的隨機性所造成的需要有更準(zhǔn)確的方法來檢驗這種差異是否顯著，也就是進行方差分析所以叫方差分析，因為雖然我們感興趣的是均值，但在判斷均值之間是否有差異時則需要借助于方差這個名字也表示：它是通過對數(shù)據(jù)誤差來源的分析判斷不同總體的均值是否相等。因此，進行方差分析時，需要考察數(shù)據(jù)誤差的來源。方差分析的基本思想和原理第十頁，共五十五頁，編輯于2023年，星期五1. 比較兩類誤差，以檢驗均值是否相等2. 比較的基礎(chǔ)是方差比3. 如果系統(tǒng)(處理)誤差顯著地不同于隨機誤差，則均值就是不相等的；反之，均值就是相等的4. 誤差是由各部分的誤差占總誤差的比例來測度的方差分析的基本思想和原理第十一頁，共五十五頁，編輯于2023年，星期五方差分析的基本思想和原理

(兩類誤差)隨機誤差因素的同一水平(總體)下，樣本各觀察值之間的差異比如，同一行業(yè)下不同企業(yè)被投訴次數(shù)是不同的這種差異可以看成是隨機因素的影響，稱為隨機誤差

系統(tǒng)誤差因素的不同水平(不同總體)下，各觀察值之間的差異比如，不同行業(yè)之間的被投訴次數(shù)之間的差異這種差異可能是由于抽樣的隨機性所造成的，也可能是由于行業(yè)本身所造成的，后者所形成的誤差是由系統(tǒng)性因素造成的，稱為系統(tǒng)誤差第十二頁，共五十五頁，編輯于2023年，星期五方差分析的基本思想和原理

(兩類方差)數(shù)據(jù)的誤差用平方和(sumofsquares)表示，稱為方差組內(nèi)方差(withingroups)因素的同一水平(同一個總體)下樣本數(shù)據(jù)的方差比如，零售業(yè)被投訴次數(shù)的方差組內(nèi)方差只包含隨機誤差組間方差(betweengroups)因素的不同水平(不同總體)下各樣本之間的方差比如，四個行業(yè)被投訴次數(shù)之間的方差組間方差既包括隨機誤差，也包括系統(tǒng)誤差第十三頁，共五十五頁，編輯于2023年，星期五方差分析的基本思想和原理

(方差的比較)若不同不同行業(yè)對投訴次數(shù)沒有影響，則組間誤差中只包含隨機誤差，沒有系統(tǒng)誤差。這時，組間誤差與組內(nèi)誤差經(jīng)過平均后的數(shù)值就應(yīng)該很接近，它們的比值就會接近1若不同行業(yè)對投訴次數(shù)有影響，在組間誤差中除了包含隨機誤差外，還會包含有系統(tǒng)誤差，這時組間誤差平均后的數(shù)值就會大于組內(nèi)誤差平均后的數(shù)值，它們之間的比值就會大于1當(dāng)這個比值大到某種程度時，就可以說不同水平之間存在著顯著差異，也就是自變量對因變量有影響判斷行業(yè)對投訴次數(shù)是否有顯著影響，實際上也就是檢驗被投訴次數(shù)的差異主要是由于什么原因所引起的。如果這種差異主要是系統(tǒng)誤差，說明不同行業(yè)對投訴次數(shù)有顯著影響第十四頁，共五十五頁，編輯于2023年，星期五方差分析的基本假定每個總體都應(yīng)服從正態(tài)分布對于因素的每一個水平，其觀察值是來自服從正態(tài)分布總體的簡單隨機樣本比如，每個行業(yè)被投訴的次數(shù)必需服從正態(tài)分布各個總體的方差必須相同各組觀察數(shù)據(jù)是從具有相同方差的總體中抽取的比如，四個行業(yè)被投訴次數(shù)的方差都相等觀察值是獨立的比如，每個行業(yè)被投訴的次數(shù)與其他行業(yè)被投訴的次數(shù)獨立第十五頁，共五十五頁，編輯于2023年，星期五方差分析中的基本假定在上述假定條件下，判斷行業(yè)對投訴次數(shù)是否有顯著影響，實際上也就是檢驗具有同方差的四個正態(tài)總體的均值是否相等如果四個總體的均值相等，可以期望四個樣本的均值也會很接近四個樣本的均值越接近，推斷四個總體均值相等的證據(jù)也就越充分樣本均值越不同，推斷總體均值不同的證據(jù)就越充分第十六頁，共五十五頁，編輯于2023年，星期五方差分析中基本假定如果原假設(shè)成立，即H0:m1=m2=m3=m4四個行業(yè)被投訴次數(shù)的均值都相等意味著每個樣本都來自均值為、方差為2的同一正態(tài)總體

Xf(X)1

2

3

4

第十七頁，共五十五頁，編輯于2023年，星期五方差分析中基本假定若備擇假設(shè)成立，即H1:mi(i=1，2，3，4)不全相等至少有一個總體的均值是不同的四個樣本分別來自均值不同的四個正態(tài)總體Xf(X)3

1

2

4

第十八頁，共五十五頁，編輯于2023年，星期五問題的一般提法設(shè)因素有k個水平，每個水平的均值分別用1、2、、k

表示要檢驗k個水平(總體)的均值是否相等，需要提出如下假設(shè)：H0:12…k

H1:1,2,，k不全相等設(shè)1為零售業(yè)被投訴次數(shù)的均值，2為旅游業(yè)被投訴次數(shù)的均值，3為航空公司被投訴次數(shù)的均值，4為家電制造業(yè)被投訴次數(shù)的均值，提出的假設(shè)為H0:1234

H1:1,2,3,4

不全相等第十九頁，共五十五頁，編輯于2023年，星期五§3.2單因素方差分析數(shù)據(jù)結(jié)構(gòu)分析步驟關(guān)系強度的測量用Excel進行方差分析第二十頁，共五十五頁，編輯于2023年，星期五單因素方差分析的數(shù)據(jù)結(jié)構(gòu)

(one-wayanalysisofvariance)

觀察值(j)因素(A)i

水平A1水平A2

…水平Ak12::n

x11x21…xk1x12x22…xk2::::::::x1n

x2n…xkn第二十一頁，共五十五頁，編輯于2023年，星期五分析步驟提出假設(shè)構(gòu)造檢驗統(tǒng)計量統(tǒng)計決策第二十二頁，共五十五頁，編輯于2023年，星期五提出假設(shè)一般提法H0:m1=m2=…=mk自變量對因變量沒有顯著影響H1:m1

，m2

，…

，mk不全相等自變量對因變量有顯著影響注意：拒絕原假設(shè)，只表明至少有兩個總體的均值不相等，并不意味著所有的均值都不相等第二十三頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量構(gòu)造統(tǒng)計量需要計算水平的均值全部觀察值的總均值誤差平方和均方(MS)

第二十四頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(計算水平的均值)假定從第i個總體中抽取一個容量為ni的簡單隨機樣本，第i個總體的樣本均值為該樣本的全部觀察值總和除以觀察值的個數(shù)計算公式為式中：ni為第i個總體的樣本觀察值個數(shù)

xij為第i個總體的第j個觀察值

第二十五頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(計算全部觀察值的總均值)全部觀察值的總和除以觀察值的總個數(shù)計算公式為第二十六頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(例題分析)第二十七頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(計算總誤差平方和SST)全部觀察值與總平均值的離差平方和反映全部觀察值的離散狀況其計算公式為前例的計算結(jié)果：

SST=(57-47.869565)2+…+(58-47.869565)2=115.9295第二十八頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(計算水平項平方和SSA)各組平均值與總平均值的離差平方和反映各總體的樣本均值之間的差異程度，又稱組間平方和該平方和既包括隨機誤差，也包括系統(tǒng)誤差計算公式為前例的計算結(jié)果：SSA=1456.608696第二十九頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(計算誤差項平方和SSE)每個水平或組的各樣本數(shù)據(jù)與其組平均值的離差平方和反映每個樣本各觀察值的離散狀況，又稱組內(nèi)平方和該平方和反映的是隨機誤差的大小計算公式為前例的計算結(jié)果：SSE=2708第三十頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(三個平方和的關(guān)系)總離差平方和(SST)、誤差項離差平方和(SSE)、水平項離差平方和(SSA)之間的關(guān)系SST=SSA+SSE前例的計算結(jié)果：

4164.608696=1456.608696+2708第三十一頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(三個平方和的作用)

SST反映全部數(shù)據(jù)總的誤差程度；SSE反映隨機誤差的大??；SSA反映隨機誤差和系統(tǒng)誤差的大小如果原假設(shè)成立，則表明沒有系統(tǒng)誤差，組間平方和SSA除以自由度后的均方與組內(nèi)平方和SSE和除以自由度后的均方差異就不會太大；如果組間均方顯著地大于組內(nèi)均方，說明各水平(總體)之間的差異不僅有隨機誤差，還有系統(tǒng)誤差判斷因素的水平是否對其觀察值有影響，實際上就是比較組間方差與組內(nèi)方差之間差異的大小第三十二頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(計算均方MS)各誤差平方和的大小與觀察值的多少有關(guān)，為消除觀察值多少對誤差平方和大小的影響，需要將其平均，這就是均方，也稱為方差計算方法是用誤差平方和除以相應(yīng)的自由度三個平方和對應(yīng)的自由度分別是SST的自由度為n-1，其中n為全部觀察值的個數(shù)SSA的自由度為k-1，其中k為因素水平(總體)的個數(shù)SSE的自由度為n-k第三十三頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(計算均方MS)

組間方差：SSA的均方，記為MSA，計算公式為

組內(nèi)方差：SSE的均方，記為MSE，計算公式為第三十四頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(計算檢驗統(tǒng)計量F)將MSA和MSE進行對比，即得到所需要的檢驗統(tǒng)計量F當(dāng)H0為真時，二者的比值服從分子自由度為k-1、分母自由度為n-k的F分布，即第三十五頁，共五十五頁，編輯于2023年，星期五構(gòu)造檢驗的統(tǒng)計量

(F分布與拒絕域)如果均值相等，F(xiàn)=MSA/MSE1a

F分布F(k-1,n-k)0拒絕H0不拒絕H0F第三十六頁，共五十五頁，編輯于2023年，星期五統(tǒng)計決策

將統(tǒng)計量的值F與給定的顯著性水平的臨界值F進行比較，作出對原假設(shè)H0的決策根據(jù)給定的顯著性水平，在F分布表中查找與第一自由度df1＝k-1、第二自由度df2=n-k相應(yīng)的臨界值F

若F>F

，則拒絕原假設(shè)H0

，表明均值之間的差異是顯著的，所檢驗的因素對觀察值有顯著影響若F<F

，則不拒絕原假設(shè)H0

，不能認(rèn)為所檢驗的因素對觀察值有顯著影響第三十七頁，共五十五頁，編輯于2023年，星期五單因素方差分析表

(基本結(jié)構(gòu))第三十八頁，共五十五頁，編輯于2023年，星期五單因素方差分析

(例題分析)第三十九頁，共五十五頁，編輯于2023年，星期五用Excel進行方差分析第四十頁，共五十五頁，編輯于2023年，星期五用Excel進行方差分析

第1步：選擇“工具”下拉菜單第2步：選擇“數(shù)據(jù)分析”選項第3步：在分析工具中選擇“單因素方差分析”

，然后選擇“確定”第4步：當(dāng)對話框出現(xiàn)時

在“輸入?yún)^(qū)域”方框內(nèi)鍵入數(shù)據(jù)單元格區(qū)域在方框內(nèi)鍵入0.05（可根據(jù)需要確定）在“輸出選項”中選擇輸出區(qū)域用Excel進行方差分析第四十一頁，共五十五頁，編輯于2023年，星期五§3.3方差分析中的多重比較多重比較的意義多重比較的方法第四十二頁，共五十五頁，編輯于2023年，星期五方差分析中的多重比較

(multiplecomparisonprocedures)通過對總體均值之間的配對比較來進一步檢驗到底哪些均值之間存在差異可采用最小顯著差異方法(LSD),鄧肯法（Duncan）第四十三頁，共五十五頁，編輯于2023年，星期五SPSS方差分析第四十四頁，共五十五頁，編輯于2023年，星期五方差分析過程單因素方差分析的菜單圖高級多元方差分析菜單第四十五頁，共五十五頁，編輯于2023年，星期五簡單的一維方差分析返回飼料比較數(shù)據(jù)飼料ABCD133.8151.2193.4225.8125.3149.0185.3224.6143.1162.7182.8220.4128.9143.8188.5212.3135.7153.5198.6第四十六頁，共五十五頁，編輯于2023年，星期五單因素方差分析--主對話框第四十七頁，共五十五頁，編輯于2023年，星期五

對照比較對話框返回第四十八頁，共五十五頁，編輯于2023年，星期五

均值多重比較的對話框

第四十九頁，共五十五頁，編輯于2023年，星期五輸出統(tǒng)計量對話框

第五十頁，共五十五頁，編輯于2023年，星期五分析輸出1返回描述統(tǒng)計量例：分析不同飼料對豬體重的影響data07-01

第五十一頁，共五十五頁，編輯于2023年，星期五多重比較結(jié)果的表示方法(一)列梯形表法(二)