在ＥＸＣＥＬ中實現(xiàn)多總體方差的Ｂａｒｔｌｅｔｔ齊性檢驗

摘 要：在總體服從正態(tài)分布的前提下，方差齊性與否決定了單因素多水平的模型能否采用方差分析來檢驗均數(shù)間的差異。主要介紹在EXCEL中，Bartlett齊性檢驗的自動計算工具的設(shè)計與實現(xiàn)。工具完成后可以重復(fù)使用，大大減少以后齊性檢驗計算的工作量。

關(guān)鍵詞： 方差齊性；總體；檢驗；EXCEL

中圖分類號： G8032 文章編號：1009-783X(2007)01-0114-03 文獻(xiàn)標(biāo)識碼： A

在體育教學(xué)和運動訓(xùn)練等的科學(xué)實驗中，對影響體育教學(xué)成績及運動競賽的成績的原因的探究，一直是當(dāng)代體育科研中研究的主線。例如，在運動訓(xùn)練中，為更加有效地提高運動成績，通常需要考察不同的運動強度、不同的運動量和不同的運動持續(xù)時間等因素對不同的專項運動成績的影響，目的是為了找出適合不同專項的運動強度、運動量、運動持續(xù)時間的較佳組合。又如，我們從運動系體操專業(yè)的學(xué)生中隨機抽取條件相似的20名學(xué)生隨機分成4組，每組5人，由4位教師施以不同的教學(xué)方法，教20個具有相當(dāng)難度的體操動作，并規(guī)定每個動作的計分標(biāo)準(zhǔn)，試教一學(xué)期后舉行測試，測得各組得分，見下表。現(xiàn)假定每組的得分服從正態(tài)分布，則這4種教學(xué)方法的效果間是否有顯著性差異的問題就是我們迫切需要了解的。

如果僅僅從上例每組的總分上看，顯然四種不同的教法帶來了四種不同的學(xué)生得分，分值上肯定有差異，但這種差異主要是由隨機誤差引起的，還是主要是由于教學(xué)方法的不同而引起的，即是否有顯著性差異的統(tǒng)計結(jié)論，還須經(jīng)統(tǒng)計檢驗后才能得出。若用兩個樣本間均數(shù)差異的顯著性檢驗方法來處理本類問題的話，需要做6次檢驗。若這樣的試驗安排共有N組，則需要做N（N-1）/2次兩兩比較，這一方面，顯然太麻煩了，另一方面，當(dāng)設(shè)定兩兩比較時，犯第一類錯誤的概率α=0.05，則N個獨立樣本兩兩比較時，每次比較不犯第一類錯誤的概率為0.95^N(N-1)／2，相應(yīng)犯第一類錯誤的概率為1-0.95^N(N-1)／2，遠(yuǎn)遠(yuǎn)大于事先設(shè)定的0.05。因此，多個均數(shù)比較時不宜采用我們熟知的t檢驗作兩兩比較，應(yīng)采用一種新的統(tǒng)計處理方法來實現(xiàn)。

解決這一類問題的方法是方差分析。它最早由英國統(tǒng)計學(xué)家費舍（R.A.Fisher）在1923年提出，最初用于生物學(xué)和農(nóng)業(yè)試驗方面，后于1946年由斯內(nèi)德克（G.W.Snedecor）進(jìn)一步加以完善。為紀(jì)念費舍的杰出貢獻(xiàn)，又把它稱為F檢驗。現(xiàn)在它在體育領(lǐng)域中也得到了廣泛的應(yīng)用。

方差分析是在總體服從正態(tài)分布且方差齊性的假設(shè)下展開的，在滿足總體正態(tài)性但方差不齊時，此法不可用，而只能改用方差不齊時兩均數(shù)差異的顯著性檢驗的方法來進(jìn)行兩兩均數(shù)間的比較。因此，這里很有必要來考慮方差的齊性檢驗的問題。本文主要介紹在EXCEL中如何來實現(xiàn)多總體方差的Bartlett齊性檢驗的自動計算。

1 Bartlett方差齊性檢驗的方法

Bartlett法是一種可在各水平重復(fù)測定次數(shù)不等時用來檢驗方差齊性的方法，雖然，當(dāng)各水平重復(fù)測定次數(shù)相等時，可用Cochran提供的檢驗方法，但Bartlett法同樣適用。

2 在EXCEL中進(jìn)行Bartlett方差齊性檢驗的方法

2.1 工作表的安排

在用Bartlett法進(jìn)行方差齊性檢驗時，為使計算相對自動化，而且能在必要時可以擴展到更多的水平數(shù)，需要在同一個工作簿中動用兩張工作表。

其中，工作表1用來存放實驗結(jié)果的原始數(shù)據(jù)。工作表2用來存放中間和最后的運算結(jié)果。

2.2 原始數(shù)據(jù)清單的格式要求

在第一張工作表中，第一行從A到Z的各列中用來建立單因素相應(yīng)各水平的名稱，即字段名。字段名應(yīng)由字母或漢字開頭加EXCEL中的其它合法字符組成。考慮到后續(xù)運作的需要，字段名的長度不要過長，最好在10個字符以內(nèi)。

A列中用來存放因素的第一個水平的字段名和測試值，B列中用來存放因素的第二個水平的字段名和測試值，余類推。建立的數(shù)據(jù)清單，如圖1所示。需要注意的是在工作表1中除輸入上述必須的數(shù)據(jù)外，不能在空白單元格中輸入與實驗數(shù)據(jù)無關(guān)的任何內(nèi)容。

2.3 在工作表2中用Excel的公式和函數(shù)進(jìn)行計算

根據(jù)Bartlett方差齊性檢驗的計算公式，用Excel公式和函數(shù)計算和存放中間統(tǒng)計量和最后結(jié)果。從而使整個計算過程實現(xiàn)自動化。工作表2中的樣式見圖2。

需要注意的是，在工作表2第12行以上的整個區(qū)域中的許多單元格，是被用來定義存放不同的計算結(jié)果的，因此，與此無關(guān)的內(nèi)容，不能在本區(qū)域中出現(xiàn)。

2.3.1 單因素水平數(shù)（即分組數(shù)）的確定

在計算過程中，多次用到水平數(shù)這個值，由于水平數(shù)會隨著實際問題的變化而變化，為達(dá)到自動計算的目的，必須用函數(shù)來對它進(jìn)行實際的測定，具體可通過函數(shù)COUNTIF(Sheet1!A:Z，\">\"\"\")來計算該問題的實際水平數(shù)。其中A:Z表示，由A到Z共26列和EXCEL默認(rèn)的最大行組成的區(qū)域，也即測定的最大水平數(shù)定義為26。本文開發(fā)的計算工具中，定義的最大水平數(shù)為26。如果實際的水平數(shù)比這還大，可通過修改Z值使之更大即可。

2.3.2 總觀測值數(shù)（各樣本含量的總數(shù)）的確定

同樣用函數(shù)COUNT(Sheet1!A:Z)可計算得到該問題中總觀測值數(shù)。同上，如果實際的水平數(shù)比這還大，可通過修改Z值使之更大即可統(tǒng)計出各樣本含量的總數(shù)。

在本計算工具的開發(fā)過程中，上述分組數(shù)和總觀測值數(shù)的計算函數(shù)將被嵌套在其它函數(shù)中用來實現(xiàn)具體的功能。

2.3.3 各水平樣本含量的確定

工作表1中，第一列即A列的樣本含量的確定， 可用COUNT(Sheet1!A:A) 計算得到該列的觀測值數(shù)。由于工作表1中的水平至少3個或3個以上，因此，須在工作表2中用多個處在同一行的連續(xù)的單元格來存放計算結(jié)果，這就帶來不確定性。當(dāng)預(yù)設(shè)有公式的單元格數(shù)超過實際的因素的水平數(shù)時，如不作特殊處理將在后續(xù)的計算中會帶來不必要的麻煩。為了使工作表2的相應(yīng)列中列出各水平的樣本含量數(shù)，而對超出對應(yīng)水平數(shù)的列中不顯示任何信息，從而不影響后續(xù)的計算，因此，可使用IF（）函數(shù)來實現(xiàn)上述想法，即當(dāng)某列中的觀測值數(shù)為0時，在工作表2的相應(yīng)單元格中，不顯示任何信息。考慮到其它各列的樣本含量的統(tǒng)計方法同其是一樣的，可通過使用填充句柄復(fù)制的方法實現(xiàn)公式的粘貼，為能得到圖2所示的結(jié)果，故在工作表2的B2單元格中輸入“=IF(COUNT(Sheet1!A:A)=0，\"\"，COUNT(Sheet1!A:A))”，其它各列與第2行對應(yīng)單元格中的計算公式可通過向右拖曳填充句柄的方式來實現(xiàn)公式的錄入。同前一致，在工作表2中，從B2單元格起已將計算各樣本含量數(shù)的公式設(shè)置到AA列的AA2元格。這樣最多可同時處理26個水平。如水平數(shù)更多時，用上法繼續(xù)向右復(fù)制粘貼公式即可。

2.3.4 各樣本自由度的確定

樣本的自由度等于樣本含量減1。工作表2中的第三行［l4］被設(shè)計用來存放計算得到的各樣本自由度的結(jié)果。主要是為了用它來計算 的結(jié)果。它可通過上一行的對應(yīng)數(shù)據(jù)減1來求得。考慮到上述的其它原因，因此在單元格B3中輸入“=IF(B2=\"\"，\"\"，B2-1) ”，即可計算得到該列樣本對應(yīng)的自由度。其它各列與第3行對應(yīng)單元格中的計算公式可通過向右拖曳填充句柄的方式來實現(xiàn)公式的錄入，一直到AA3單元格。

2.3.5 各樣本方差的計算

在EXCEL中，樣本的方差可通過調(diào)用VAR()函數(shù)來計算。區(qū)域A:A表示整個A列，函數(shù)VAR(A:A)可用來求A列樣本的方差。因此，在B4單元格中輸入“=IF(COUNT(Sheet1!A:A)=0，\"\"，VAR(Sheet1!A:A))”，其它各列與第4行對應(yīng)單元格中的計算公式可通過向右拖曳填充句柄的方式一直到AA4單元格結(jié)束來實現(xiàn)公式的錄入。

2.3.6 離均差平方和的計算

離均差平方和這一行的設(shè)置是專門用來計算合并方差的。在EXCEL中，調(diào)用DEVSQ()函數(shù)可計算各樣本的離均差平方。同樣本方差的計算方法一樣，只需在B5單元格［16］中輸入“=IF(COUNT(Sheet1!A:A)=0，\"\"，DEVSQ(Sheet1!A:A))”，其它各列與第5行對應(yīng)單元格中的計算公式可通過向右拖曳填充句柄的方式一直到AA5單元格結(jié)束來實現(xiàn)公式的錄入。

2.3.7 方差的對數(shù)與自由度乘積的計算

方差的對數(shù)與自由度乘積的計算是為計算q值做準(zhǔn)備的。在EXCEL中，以10為底的方差的對數(shù)的計算，可直接調(diào)用LOG10( )函數(shù)來完成。因此，只需在B6單元格中輸入“=IF(B2=\"\"，\"\"，LOG10(B4)*B3)”，即可得到第一個樣本的方差的對數(shù)與其自由度乘積的計算結(jié)果，其它各列與第6行對應(yīng)單元格中的計算公式的錄入，可通過向右拖曳填充句柄的方式一直到AA6單元格結(jié)束來實現(xiàn)。

2.3.8 合并方差的計算

合并方差的計算是為計算最終的χ²值做準(zhǔn)備的。根據(jù)組合樣本方差的計算公式，在B7單元格中輸入“=SUM(B5:I5)/(COUNT(Sheet1!A:Z)-COUNTIF(Sheet1!A:Z，\">\"\"\"))”，即可求得合并方差。

2.3.9 q值的計算

總觀測值數(shù)與水平數(shù)之差可通過“=COUNT(Sheet1!A:Z)-COUNTIF(Sheet1!A:Z，\">\"\"\")”來求得，方差的對數(shù)與自由度乘積的總和可用“=SUM(B6:AA6)”來求得，故在B8單元格中輸入“=(COUNT(Sheet1!A:Z)-COUNTIF(Sheet1!A:Z，\">\"\"\"))*LOG10(B7)-SUM(B6:AA6)”，即可求得q值。

2.3.12 P值的計算

調(diào)用卡方分布函數(shù)CHIDIST()，可計算得到給定的χ^{2值在給定的自由度下的概率值。在B11單元格中輸入“=CHIDIST(B10，COUNTIF(Sheet1!A:Z，\">\"\"\")-1)”，即可求得原假設(shè)成立的概率值。}

2.3.13 得出分析結(jié)論

當(dāng)p<0.05時，拒絕方差齊性的原假設(shè)，當(dāng)p>=0.05時，沒有理由拒絕方差齊性的原假設(shè)。故在B12單元格中輸入“=IF(B11>0.05，“方差齊性”，“方差不齊”)”。

我們將上述公式在“方差齊性檢驗”工作簿的工作表2中對應(yīng)地做了設(shè)置，并用8水平的一組數(shù)據(jù)進(jìn)行測試，運行正常，用表1中的數(shù)據(jù)進(jìn)行了校驗，證明工作表2中預(yù)設(shè)的公式是正確的。表明上述的思路是正確的。

3 實例分析

在田徑的鉛球教學(xué)中，為了研究三種不同的鉛球教學(xué)方法的效果，將某年級三個班中，同齡的各種運動能力基本相同的男生分成三個組，分別按以下三種不同的方法進(jìn)行教學(xué)，方法一（A1）、方法二（A2）、方法三（A3）。在三個月后，經(jīng)過多次課的教學(xué)，以同樣的標(biāo)準(zhǔn)測得各組成績，見圖3。試分析三種方案的教學(xué)效果的方差間是否齊性？

在EXCEL中的解題步驟：

1打開“方差齊性檢驗”的工作簿。

2在工作表1中按圖三所示輸入字段名和原始數(shù)據(jù)。

3單擊工作表2的標(biāo)簽，在工作表2中出現(xiàn)如表4所示的結(jié)果，不拒絕方差齊性的原假設(shè)。

4自動計算工具制作的注意事項。

1) 原始數(shù)據(jù)的錄入應(yīng)嚴(yán)格按文中“原始數(shù)據(jù)清單的格式要求”中的有關(guān)說明進(jìn)行，方可在工作表2中得到正確的運算結(jié)果。

2）工作表2的許多單元格中的IF()函數(shù)不能缺少，它不僅是用來美化工作表2的，更是用來確保計算的正確性的。

3) 注意函數(shù)與公式輸入的正確性，這是自動檢驗工具完成的關(guān)鍵。尤其應(yīng)該注意公式和函數(shù)中的標(biāo)點符號必須是英文半角符號。

4 小結(jié)

1）本文通過對Bartlett齊性檢驗中計算公式的分解，合理利用EXCEL中的函數(shù)，實現(xiàn)了多總體方差的Bartlett齊性檢驗的自動計算。

2）方差齊性自動檢驗工具一旦制作完成，即可重復(fù)使用。使用簡單方便，幾乎一勞永逸。大大提高了運算速度。

3）可以根據(jù)教學(xué)對象，引導(dǎo)學(xué)生按此文敘述的步驟自己實現(xiàn)自動計算工具的制作，教師也可以將自己制作的檢驗工具拷貝給學(xué)生，供學(xué)生使用。

本文中所涉及到的圖表、注解、公式等內(nèi)容請以PDF格式閱讀原文。