SPSS詳細(xì)操作：樣本均數(shù)間的多重比較

下面我們用例子來探討：方差分析得出各組總體均數(shù)間的差異有統(tǒng)計學(xué)意義之后，如何進(jìn)行樣本均數(shù)間的多重比較。

一、問題與數(shù)據(jù)

為調(diào)查A、B、C三種治療措施對患者谷丙轉(zhuǎn)氨酶（ALT）的影響，某科室將45名患者隨機(jī)分為三組，每組15人，分別采取A、B、C三種治療措施，治療后ALT水平（U/L）如下。試問應(yīng)用三種治療措施后，患者的ALT水平是否有差異？如果有差異，又是哪兩個組間的ALT水平存在差異？

表1. 三組患者治療后的ALT水平（U/L）

二、對數(shù)據(jù)結(jié)構(gòu)的分析

整個數(shù)據(jù)資料涉及三組患者，每組15人，測量指標(biāo)為血常規(guī)報告的ALT水平，因此屬于多組設(shè)計的定量資料。

我們已經(jīng)知道可以采用單因素方差分析的方法，比較三組總體均數(shù)之間的差異是否具有統(tǒng)計學(xué)意義，但如果差異有統(tǒng)計學(xué)意義，并不意味著任意兩組之間的均數(shù)差異都有統(tǒng)計學(xué)意義。想進(jìn)一步了解哪兩個組間的ALT水平存在差異，還需要做樣本均數(shù)間的多重比較。 

三、SPSS分析方法

1. 數(shù)據(jù)錄入SPSS

2. 進(jìn)行單因素方差分析：

（1）Analyze → General Linear Model → Univariate；

（2）主對話框設(shè)置：分析變量（ALT）送入Dependent Variable 框中→分組變量（Group）送入Fixed Factor(s) 框中；

（3）Options設(shè)置：勾選Descriptive statistics（統(tǒng)計描述）和Homogeneity tests（方差齊性檢驗(yàn)）→Continue。

3. 均數(shù)間的多重比較：

（1）點(diǎn)擊Post Hoc按鈕，將Group送入Post Hoc Tests for框中→在Equal Variances Assumed（方差齊）框中，勾選幾種常用的多重比較方法：

LSD、S-N-K、Bonferroni、Tukey、?idák、Scheffé's、Dunnett

（2）在Dunnett選項(xiàng)下面，Control Category選擇First（以A組作為對照組），Test框中勾選2-sided→Continue→OK。

四、結(jié)果解讀

Tests of Between-Subjects Effects表格給出了方差分析的結(jié)果。在方差齊的條件下，Group一行結(jié)果顯示，F(xiàn)值=68.810， P（Sig.）<0.001。

Multiple Comparisons表格給出了部分方法的多重比較結(jié)果，分別列出了每個組和其他組比較的均數(shù)的差值（Mean Difference (I-J)）、標(biāo)準(zhǔn)誤（Std. Error）、P值（Sig.）和均數(shù)差值的95%置信區(qū)間（95% Confidence Interval）。檢驗(yàn)水準(zhǔn)α設(shè)為0.05，組間差異有統(tǒng)計學(xué)意義的結(jié)果已用*標(biāo)出。

不同多重比較方法的選擇，需要結(jié)合研究設(shè)計和每個方法各自的特點(diǎn)及適用條件。我們以Bonferroni法和Dunnett法的結(jié)果為例，進(jìn)行解讀：

Bonferroni法結(jié)果顯示，A組與B組的ALT水平相比，Mean Difference=-15.160 U/L，P（Sig.）<0.001；A組與C組相比，Mean Difference=1.133 U/L，P（Sig.）=1.000；B組與C組相比，Mean Difference=16.293 U/L，P（Sig.）<0.001。

Bonferroni法和?idák法的檢驗(yàn)原理為根據(jù)需要比較的次數(shù)來調(diào)整檢驗(yàn)水準(zhǔn)。以Bonferroni法為例，調(diào)整的方法有兩種：1、將α’=α/m作為檢驗(yàn)水準(zhǔn)的調(diào)整值（m為共需比較的次數(shù)，假設(shè)m=C32=3，α’=0.05/3=0.017），兩兩比較得出的P值與α’進(jìn)行比較，即P<0.017時才能拒絕零假設(shè)；2、將計算得到的P值擴(kuò)大為原來的m倍為P’，然后將P’與0.05比較，即P’<0.05時拒絕零假設(shè)。SPSS采取的是第二種方法，直接為我們給出了Bonferroni法和?idák法調(diào)整后的P’值，故將調(diào)整后的P’值與檢驗(yàn)水準(zhǔn)0.05比較即可，不用再自行計算。

Dunnett法結(jié)果顯示，B組與A組的ALT水平相比，Mean Difference=15.160 U/L，P（Sig.）<0.001；C組與A組相比，Mean Difference=-1.133 U/L，P（Sig.）=0.688。

Homogeneous Subsets表格同樣給出了部分方法的多重比較結(jié)果。我們以Student-Newman-Keuls（S-N-K）法的結(jié)果為例，進(jìn)行解讀：

C組ALT水平的均數(shù)為12.147 U/L，A組均數(shù)為13.280 U/L，C組和A組被分到一個亞組（Subset 1），P（Sig.）=0.469；B組均數(shù)為28.440 U/L，被單獨(dú)分到另一個亞組（Subset 2），P（Sig.）>0.999。

五、撰寫結(jié)論

A組患者ALT水平為（13.28 ± 4.39）U/L，B組患者ALT水平為（28.44 ± 3.65）U/L，C組患者ALT水平為（12.15 ± 4.64）U/L，A、B、C三種治療措施對患者ALT水平的影響差異具有統(tǒng)計學(xué)意義（F=68.810，P<0.001）。其中，B組與A組，B組與C組間患者的ALT水平差異均具有統(tǒng)計學(xué)意義（P<0.001），A組與C組間患者的ALT水平差異無統(tǒng)計學(xué)意義（P>0.05）。

六、延伸閱讀

SPSS中均數(shù)間多重比較方法可以分為兩大類，包括多重兩兩比較的檢驗(yàn)（Multiple Comparisons）和對極差進(jìn)行亞組同質(zhì)性的檢驗(yàn)（Homogeneous Subsets），二者結(jié)果的表現(xiàn)形式有所不同，如本例結(jié)果所示。

兩兩比較檢驗(yàn)給出了每個組和其他組比較的結(jié)果。亞組同質(zhì)性檢驗(yàn)則先按照每個組均數(shù)的大小進(jìn)行排序，將組間均數(shù)差異無統(tǒng)計學(xué)意義的均數(shù)分到一個亞組，這樣會形成若干個亞組，同一亞組內(nèi)的均數(shù)差異無統(tǒng)計學(xué)意義，不同亞組間的均數(shù)差異有統(tǒng)計學(xué)意義。

本例中，C組和A組被分到一個亞組，說明兩組間均數(shù)差異無統(tǒng)計學(xué)意義；B組被單獨(dú)分到一個亞組，說明B組與A組、B組與C組間的均數(shù)差異都有統(tǒng)計學(xué)意義。