臨床場景:一位言語治療師想知道,12 週的舌頭阻力訓練(TSE)能否改善社區老年人的吞嚥功能。她在第 0、4、8、12 週分別測量每位個案的最大舌壓。問題來了:這四個時間點的測量值,都來自同一個人,它們之間絕對不是「各自獨立的觀察」。
重複測量為什麼違反獨立性假設?
從第2章到第17章,我們使用的大部分統計方法,背後都有一個核心假設:每一筆觀察值都是彼此獨立的。也就是說,第 1 號個案的資料不應該影響第 2 號個案的分析結果。
但當我們對同一個人重複測量時,同一個人在不同時間點的分數,必然會互相關聯。天生舌頭力量較強的人,在第 0、4、8、12 週的測量值都傾向偏高;反之,基線功能較差者,各次測量也都傾向偏低。這種組內相關(within-person correlation)組內相關(within-person correlation)指來自同一個體不同時間點的觀察值之間的相關性。由於這些數值共享同一個人的「個人特性」,它們之間必然存在正相關。這違反了普通迴歸的獨立性假設,若忽略它,標準誤將被低估,導致假陽性增加。若被忽略,後果是標準誤被低估,p 值偏小,假陽性率升高。
核心問題
「我能不能直接把所有時間點的資料疊在一起,然後用普通迴歸分析?」
不行。這樣做會讓統計程式「以為」你有更多獨立觀察值,進而低估標準誤、誇大統計顯著性。結論可能看起來很漂亮,但其實是統計上的假象。
直式資料(Long Format)vs 橫式資料(Wide Format)
在進行 GEE 分析之前,必須先確認資料格式。大多數統計方法(包括 GEE)要求資料採用直式格式(Long Format),即每一列只記錄一個時間點的觀察值。
橫式資料(Wide Format)
每位個案只有一列,不同時間點用不同欄位表示
| 個案編號 | 舌壓_W0 | 舌壓_W4 | 舌壓_W8 | 舌壓_W12 |
|---|---|---|---|---|
| 001 | 18.2 | 21.4 | 25.1 | 27.8 |
| 002 | 15.6 | 16.9 | 19.3 | 22.0 |
| 003 | 22.1 | 23.5 | 26.7 | 29.4 |
直式資料(Long Format)
每個時間點占一列,GEE 需要這種格式
| 個案編號 | 測量週數 | 舌壓 |
|---|---|---|
| 001 | 0 | 18.2 |
| 001 | 4 | 21.4 |
| 001 | 8 | 25.1 |
| 001 | 12 | 27.8 |
| 002 | 0 | 15.6 |
| 002 | 4 | 16.9 |
在橫式資料中,如果某位個案某次測量缺失,整列就很難處理。而直式資料只需刪除那一列,其他時間點的資料仍然完整可用。這正是 GEE 的一大優勢。
GEE 如何解決問題?
廣義估計方程式(GEE,Generalized Estimating Equations)GEE 由 Liang 與 Zeger 於 1986 年提出。它是一種半母數(semi-parametric)方法,不需要對結局變項的分布做嚴格假設,只需指定平均值與共變數結構之間的關係,以及一個「工作相關矩陣」來描述重複測量之間的相關性。是由 Liang 與 Zeger(1986)提出的半母數方法。它的核心思想很直觀:
群體平均效應
GEE 估計的是群體層次(population-averaged)的效應。例如:「整體而言,接受 TSE+CTAR 訓練的族群,其舌壓在 12 週後比基線高出多少?」這對臨床政策制定非常有用。
半母數方法
GEE 不要求結局變項必須符合常態分布。它可以處理連續變項、二元變項、計數變項、順序變項,彈性極高。
缺失資料處理
只要缺失資料的機制是完全隨機缺失(MCAR),GEE 仍然可以使用剩餘的觀察值進行有效估計,不需要強制每位個案都有完整的所有時間點資料。
本章學習目標
- 說明重複測量資料違反獨立性假設的原因,以及其統計後果。
- 比較直式(Long)與橫式(Wide)資料格式,並能進行格式轉換。
- 理解 AR(1)、Exchangeable、Unstructured 等工作相關矩陣的意涵與適用情境。
- 以 QIC 指標選擇最適工作相關矩陣。
- 正確解讀 GEE 中「時間主效應」、「組別主效應」與「時間×組別交互作用」的臨床意義。
- 在 SPSS 中完成 GEE 分析,並將結果對應至 Chou 等人(2025)的研究發現。