昨天的文章，我們對方差分析的整體邏輯進行了初步的介紹，今天將以單向（one-way）方差分析為例，具體梳理方差分析的整個過程。

單向（one-way）方差分析，就是大家很熟悉的單因素方差分析（教科書上叫單向）, 一般也稱完全隨機設(shè)計(completely randomized design)的方差分析，是指將研究對象通過完全隨機化方法，分配至多個不同的處理組，比較多組的效應(yīng)指標是否存在差別。

先看如下案例：

為了解大骨節(jié)病與糧食中微量元素硒含量之間的關(guān)系，某研究團隊調(diào)查了A（渭源縣）、B（青州市）兩個大骨節(jié)病區(qū)和C（泰山區(qū)）、D（長清區(qū)）兩個非大骨節(jié)病區(qū)。

每個病區(qū)隨機抽取20戶農(nóng)戶并采集面粉，檢測面粉中硒元素含量(μg/kg)，試分析這4個地區(qū)面粉中硒含量是否存在差異。具體的數(shù)據(jù)情況如下表1。

表1 四地區(qū)面粉硒元素含量樣本數(shù)據(jù)表

我們將上述數(shù)據(jù)繪制成圖形（如下圖，每個空心小圓圈代表一個樣本值），可以很直觀地看到，這80個樣本值（20*4）各不相同，即它們存在差異。

暫時忽略其他潛在的混雜因素，這種差異的原因可能是由于它們來自不同的地區(qū)，但因為四個小組內(nèi)部的數(shù)值也都一一不同，所以，差異也可能僅僅是因為隨機誤差，通俗地理解就是人們說的運氣導致的。

不過，仔細地觀察發(fā)現(xiàn)兩個病區(qū)的數(shù)據(jù)好像明顯要低一些，這便提示地區(qū)的不同確實有可能造成了目前的差異。

為了驗證我們的猜測，就可以采用方差分析來檢驗：病區(qū)與非病區(qū)面粉硒含量的差異是否具有統(tǒng)計學意義。

這里需要再明確一點的是，我們的目標是比較這四個地區(qū)面粉中硒含量是否有差異，在實際操作中，我們比較的是四個地區(qū)硒含量的總體平均數(shù)，因此，只要總體平均數(shù)有差異，我們就說四地區(qū)硒含量有差異。

要進行方差分析，當然，我們首先要進行假設(shè)：這四組數(shù)據(jù)都沒有差異，注意是都沒有！

在這個假設(shè)下，我們可以把這四組數(shù)據(jù)看做是一個大組，即將上述80個數(shù)據(jù)視為一個整體。對于這個整體，我們可以計算一個平均數(shù)和標準差，即表1中72.22和20.00。

可是實際情況下，這80個數(shù)據(jù)是分屬于四個小組的，因此我們也可以分別計算這四個小組的平均數(shù)，即57.11、55.58、85.62、90.55。

如果假設(shè)成立（即四組數(shù)據(jù)都沒有差異），那么這四個小組的平均數(shù)應(yīng)該是圍繞著整體平均數(shù)（即72.22）上下波動的，互相差異應(yīng)該不是很大。

但現(xiàn)在我們實際獲得的數(shù)據(jù)卻顯示，小組中最低的均數(shù)為57.11、最大的為90.55，直觀上看都與72.22的距離不小，所以我們就會懷疑不能把它們看做一個整體（更嚴謹?shù)谋磉_是，它們不是來自同一個總體），從而拒絕它們相同的假設(shè)。

順著這個思路，我們獲得下面這張表格（表2）：

表2 方差分析一般結(jié)果表

上述表格中涉及的具體的計算過程大家不需要細看，只需大概了解所謂的“離均差平方和（SS）”和“均方（MS）”的計算方法。

然后對照前面我們談到的方差的概念和計算方法，你是否會發(fā)現(xiàn)，不明所以的“均方（MS）”其實可以看做是一種特殊類型的“方差”！

對照下圖方差的計算公式：左邊是離均差平方和，右邊是自由度。

我們首先看衡量“組內(nèi)變異”的均方（MS組內(nèi)），在各個小組內(nèi)部因為沒有地區(qū)差異，所以MS組內(nèi)的大小大小僅反映了隨機誤差（即運氣）的情況。

而當假設(shè)四組總體均數(shù)相同時，組間均方（MS組間）的大小也僅僅反映隨機誤差的大小。

如何理解這里的談到的“隨機誤差”呢？

對于來自同一個總體的兩個或多個樣本的差異我們可以簡單理解為隨機誤差，也就是說，雖然表面上它們不同，但這種不同并沒有意義，也不會反映額外的信息，僅僅是因為運氣導致的。

好比你和別人擲色子比大小，雖然你們獲得的點數(shù)不同，但這種不同完全是因為運氣，不能說明任何其他的問題，同時也不會出現(xiàn)一個人總贏、另一個總輸?shù)木置妗?/span>

但是如果對方悄悄在色子上做了點手腳，知道怎樣擲可以獲得大的點數(shù)，這個時候，你們點數(shù)的差異就不再僅僅因為運氣，還有色子的原因。

用統(tǒng)計學的語言來說，就是你們的差異不僅包含了隨機誤差，還包含了其他因素。所以，下次如果感覺自己總是在輸，就要看看是不是色子出現(xiàn)了問題。

回到本例，既然MS組間和MS組內(nèi)這兩者都僅反映隨機誤差的大小，那么其攜帶的信息量就應(yīng)該沒有差異（提示：方差的大小決定了數(shù)據(jù)的信息量）。

由此，在數(shù)值上MS組間與MS組內(nèi)差異不大，所以使用MS組間除以MS組內(nèi)時，所獲得F值原則上應(yīng)該在1附近。

如果現(xiàn)在我們獲得的數(shù)據(jù)計算出來的F值比1大很多（對應(yīng)的P值會很?。?，則意味著MS組間遠大于MS組內(nèi)，從而表示，MS組間攜帶了多余的信息，因此，可以證明MS組間的差異不僅僅包含隨機誤差，還包含其他因素（比如地區(qū)不同），結(jié)合本例，即意味著四地區(qū)間面粉硒含量不都相同！

將上述數(shù)據(jù)用SPSS運算后獲得結(jié)果如下表3。很顯然，F值超過46，遠大于1（注意，在正式情況下F值并非和1比較，這里僅為了方便理解），其對應(yīng)的P值遠小于0.05，由此拒絕零假設(shè)，差異有統(tǒng)計學意義，可以認為這四個地區(qū)硒含量的總體均數(shù)不全相等，也就是說至少有兩個地區(qū)總體均數(shù)不等。

表3 四地區(qū)硒含量方差分析結(jié)果表