
一、案例介紹
為研究某誘導方法對患者收縮壓的影響,將手術要求基本相同的15名患者,在手術過程中采用A麻醉誘導方法,并在T0(誘導前)、T1、T2、T3、T4五個時刻測量患者的收縮壓,試著進行方差分析。
二、問題分析
可以從案例得到,樣本數(shù)據(jù)包括誘導前T0、T1、T2、T3、T4五個時刻,所以考慮使用重復測量方差進行單因素重復測量方差,對于重復測量實驗數(shù)據(jù)的方差分析一般需要考慮組內(nèi)的影響。
三、軟件操作及結(jié)果解讀
(一) 數(shù)據(jù)導入
1.數(shù)據(jù)格式
重復測量方差數(shù)據(jù)的特殊之處在于需要有ID號,比如病例號等,以及時間點數(shù)據(jù),同一個ID(病例)會有多個不同時間點的數(shù)據(jù),比如10個樣本分別有5個相同測量時間點,一個樣本有5個,那么就有10*5=50行數(shù)據(jù)。數(shù)據(jù)格式如下:

2.導入數(shù)據(jù)
將整理好的數(shù)據(jù)上傳至SPSSAU系統(tǒng)內(nèi),如下:

上傳結(jié)果如下:

(二) 適用條件判斷
1.異常值
和大多數(shù)分析一樣,一般在分析前都需要進行一些檢驗,重復測量方差需要檢查數(shù)據(jù)是否存在異常值,以及因變量是否滿足正態(tài)分布。
異常值分析數(shù)據(jù)中不正常的值,也稱離群值。異常值的檢驗方法有很多,比如可以采用描述統(tǒng)計判斷法(一般認為介于3倍標準差外的數(shù)據(jù)),或者箱線圖等。這里使用箱線圖進行判斷。利用SPSSAU描述分析,結(jié)果如下:

從箱線圖可以粗略的判斷沒有異常值。接下來查看因變量是否滿足正態(tài)分布。
2. 正態(tài)性檢驗
判斷數(shù)據(jù)是否符合正態(tài)分布有很多種方法,比如正態(tài)性檢驗(最為嚴格),描述統(tǒng)計法(查看偏度和峰度),以及圖示法等等,這里使用直方圖進行判斷。

發(fā)現(xiàn)“收縮壓”呈現(xiàn)倒扣的鐘形,所以滿足正態(tài)性分布標準,數(shù)據(jù)滿足方差分析的前提條件,可進行重復測量方差。
(三)單因素重復測量方差分析
1. 軟件操作
重復測量方差的分析路徑為點擊【實驗醫(yī)學/研究】→【重復測量方差】然后進行分析:

2. 結(jié)果解讀
數(shù)據(jù)滿足分析條件,接下來進行重復測量方差,首先查看球形度檢驗的結(jié)果:

球形度檢驗可以看到p值大于0.05,所以通過球形度檢驗不需要校正p值?!咎貏e提示:進行球形度檢驗p值大于0.05則通過檢驗不需要校正p值,如果p值小于0.05,則沒有通過球形度檢驗需要校正p值,接下來在看球形度W值,若w值小于0.75則使用GG校正,反之使用HF校正】。
組內(nèi)分析:

根據(jù)球形度檢驗發(fā)現(xiàn)需要校正p值,從上表可以看出主效應的p值小于0.05,呈現(xiàn)0.05水平顯著,使用A麻醉誘導方法不同時間點的收縮壓有顯著性差異分析,最終結(jié)果如下:

均方(效應)
平方和SS(效應)/df(效應),例如:1110.28/4=277.57;
均方(誤差)
平方和SS(誤差)/df(誤差),例如:480.92 /4=13.359;
F
F值=均方(效應)/均方(誤差),例如:277.57/13.359=20.778;
P
P值是由F值得到的。
ges(generalized eta-squared)
平方和SS(效應)/【平方和SS(效應)+平方和SS(誤差組內(nèi))+平方和SS(誤差組間)】,例如:1110.28 /(1110.28+1136.48+480.92)=0.407;
偏Eta方
偏Eta方=平方和SS(效應)/[平方和SS(效應)+平方和SS(誤差)],例如:745176.32/(745176.32+1136.48)=0.998;
如果主效應顯著想要進一步兩兩之間比較,可以使用事后多重比較進行分析。SPSSAU提供多種事后多重比較方法,對于他們之間的差異可以查看幫助手冊,這里選擇Tukey方法(各組別的樣本數(shù)量相同時使用),假如想要了解不同時間點之間抑郁程度是否有差別。結(jié)果如下:

通過事后多重比較發(fā)現(xiàn)T0和T1、T0和T3、T0和T4等(標黃區(qū)域)的數(shù)據(jù)p值小于0.05,即存在差異。說明T0和T1等兩個時間測量的收縮壓有明顯差異。
四、結(jié)論
對案例數(shù)據(jù)進行重復測量方差分析,首先整理成正確的數(shù)據(jù)格式,然后分析前檢查數(shù)據(jù)是否滿足相關檢驗,接著對數(shù)據(jù)進行分析,分析結(jié)果得到數(shù)據(jù)不滿足球形度檢驗發(fā)現(xiàn)需要校正p值,并且發(fā)現(xiàn)使用A麻醉誘導方法不同時間點的收縮壓有顯著性差異分析,利用事后多重比較進行具體查看,發(fā)現(xiàn)T0和T1、T0和T3、T0和T4等p值小于0.05,即存在差異。
五、知識小貼士
1.SPSSAU對事后多重比較,簡單效應的表格輸出說明如下:
第一、事后多重比較規(guī)則為:針對組內(nèi)項和組間項(如果有)均輸出事后多重比較;
第二、簡單效應輸出規(guī)則為:僅針對第1個組內(nèi)項,分別與組間項(如果有)進行簡單效應表格輸出,如果組內(nèi)項超過1個,可通過互換組內(nèi)項的位置輸出其余組內(nèi)項與組間項的簡單效應表格。
2.事后多重比較與‘單獨進行事后多重比較’結(jié)果不一致?
單獨進行事后多重比較(進階方法->事后多重比較法)時,模型實質(zhì)上為單因素方差,僅考慮1個X的情況,標準誤差的計算并不一致,因此結(jié)果會不一致,但通常情況下結(jié)論會保持一致;以及此處事后多重比較使用的是邊際估計均值(偏最小二乘均值)與一般意義上的平均值有所區(qū)別,類似于SPSS軟件的EMMEANS功能。SPSSAU默認針對組內(nèi)項和組間項均進行事后多重比較。