晨晰統計部落格新站（統計、SPSS、BIG DATA討論園地）

上一篇曾經介紹過介入性研究常使用的統計方法，有獨立樣本t檢定、成對樣本t檢定、單因子共變數分析、二因子混合設計變異數分析等，除了這些方法之外，本篇再介紹一個最近比較受歡迎的方法—廣義估計方程（Generalized estimating equation, GEE）。

嚴格來說，GEE只能算是一種估計方法而非統計方法，且建立在廣義線性模式（Generalized Linear Models, GedLM）來進行分析。那麼本篇將以一個2×2的介入性研究（兩組前後測）為例子，分享如何在SPSS以GEE進行分析。

一、資料介紹

以下圖來說，第一個欄位是受試者的編號，第二個欄位為組別（1=實驗組；0=對照組），第三個欄位為time（1=前測；0=後測），第四個與第五個欄位則為outcome。

SPSS要跑GEE的資料，一定要為長資料（long form）的形式，何謂長資料的格式呢?也就是每一列僅代表每一位受試者的某一次資料，以本篇所舉的例子來說，因為每一位受訪者都有2次測驗，因此可以看到在編號欄位都會每2筆是一樣的編號，所以編號欄位是相當重要的，因為這樣在跑分析時，SPSS才有辦法辨別哪幾筆資料是來自於同一位受試者，另外再去搭配一個欄位（time）去呈現此筆資料是屬於該受試者的第幾次資料。

另一個要注意的點是，當某一組的編號完成後，要續編下一組時，流水號要接著往下編（如下圖編號17的地方），也就是在SPSS裡，每一位受試者的編號都要不同，這樣SPSS才有辦法辨別，當然您也可以歷用身份証字號、學號、病歷號來替代。

二、操作介紹

1.點選分析→廣義線性模式→廣義估計方程。

（一）重覆

2.將ID欄位丟入受試者變數

3.共變異數矩陣選到穩健估計（預設值）

4.工作相關矩陣選到AR(1) or 可交換式（Exchangeable）

（二）反應變項

5.將依變項丟入依變項欄位

6.分配選到常態（預設值），連結函數選到單位（預設值），此為線性函數的意思

（三）預測變項

7.將組別與時間丟到因子欄位，因子欄位會自動將類別變項轉成虛擬變數

8.到因子的選項中，將下方選到遞減選項，此代表組別與時間這兩個變數，將以數字較小作為參照組，也就是說用對照組與前測作為參照組

（四）模式

9.將組別與時間以主效果丟入右邊模式欄位

10.將組別與時間再以交互作用效果丟入右邊模式欄位

（五）EM平均數

11.將左邊的因子全部選到右邊顯示平均數，並將比對方式全部選成配對比較

12.可依照研究者需求調整配對比較方法，即可按下確定進行分析