可丁寧濃度隨著懷孕週數的變化：GEE 分析

就可丁寧濃度的資料而言，可以視為重複量測（repeated-measure）的資料，而在 三次的訪視收樣中，每一次都記錄了孕婦當時的懷孕週數。在此可以利用廣義線性模 式（generalized linear model）中處理縱貫式資料（longitudinal data）的分析方法－廣 義估計式（generalized estimating equation，GEE），來分析並估計懷孕週數的增加對 可丁寧濃度的影響有多大。

1.模型與變項說明

在GEE 模式中，結果變項為可丁寧濃度，解釋變項為懷孕週數。在這裡將懷孕

週數以四種不同的尺度來分組：

(1)將懷孕週數未滿 12 週者視為基準組，以後間隔 8 週分出一組，共分出 5 組；

(2)基準組設定同(1)，但是分組間隔週數改為 4 週，共分出 9 組；

(3)將懷孕週數未滿 10 週者視為基準組，以後每間隔 10 週分出一組，共分出 5 組；

(4)基準組設定同(3)，但是分組間隔週數改為 5 週，共分出 8 組。

代表懷孕週數的變項將依照組別數目設立虛擬變項（dummy varaible），當懷 孕週數有k 組時，將利用（k-1）個虛擬變項來代表懷孕週數組別。以下為我們所 使用的廣義線性模型（generalized linear model）：

G{E(Y)}=β0+β1X1+β2X2+…+βk-1Xk-1+γZ Y：尿液或血液可丁寧濃度

X1~Xk-1：代表懷孕週數分組的（k-1）個虛擬變項 Z：代表其他需要放入模式調整的變數

G{ }‧ ：identity link

就尿液與血液可丁寧量測而言，每個孕婦最多可有三次完整的實驗分析結果，但 由於追蹤漏失的問題，使得某些孕婦會有一期或兩期的數值是遺漏值，而這些數值就 無法納入統計分析。在此GEE 分析中所用的相關結構為 AR(1)（auto-regression with first order correlation）。

2.全體孕婦可丁寧濃度與懷孕週數之 GEE 分析

不同懷孕週數之尿液、血液可丁寧濃度平均值的GEE 估計：全部樣本【表 11、表 12】

由於GEE 模式中的解釋變項為代表懷孕週數組別的虛擬變項，因此虛擬變項的 迴歸係數估計值就代表著該懷孕週數組別和基準組之間可丁寧濃度平均值的差異，而 截距項就代表著懷孕週數基準組的可丁寧濃度平均值。當基準組為懷孕週數12 週以 下時，若每間隔8 週分出一組，則可以發現懷孕週數離基準組越遠者，其可丁寧濃度 平均值的差異也就越大，而且不管是哪一組，可丁寧濃度平均值和基準組相比都有顯 著的差異（12~20 週除外）。就每一組的可丁寧濃度平均值來看，懷孕週數 36 週以上 者尿液可丁寧濃度平均值為3.283 ng/mL，懷孕週數 28 週以上但未滿 36 週者平均值

為2.802 ng/mL，懷孕週數 20 週以上但未滿 28 週者平均值為 2.389 ng/mL，懷孕週數 12 週以上但未滿 20 週者平均值為 1.652 ng/mL，懷孕週數未滿 12 週者平均值為 1.221 ng/mL，由此可以很明顯的看出懷孕週數越高，尿液可丁寧的濃度也就越大。而血液 可丁寧隨著懷孕週數而增加的趨勢和尿液可丁寧相似，不過就相同的懷孕週數相比，

尿液可丁寧濃度的平均值都比血液可丁寧來得高。

然而，當懷孕週數分組方式改為每4 週分出一組時，某些組別和基準組可丁寧平 均值的差異就不再顯著了，這有可能是因為分組太細，而有些組別人數太少所造成。

除了懷孕週數16 至 24 週這一段以外，其他組別的可丁寧濃度平均值均顯示出隨著懷 孕週數而增加的趨勢。

當基準組設定改為懷孕週數未滿10 週，而分組方式改為每 10 週分出一組以及每 5 週分出一組時，所得到的結果也是相當類似的。當分組的間隔較大時，各組可丁寧 濃度平均值隨著懷孕週數而增加的趨勢較明顯；而分組的間隔較小時，除了15 至 25 週這一段外，其他組別可丁寧濃度平均值的隨著懷孕週數增加的趨勢也和之前相似。

尿液、血液可丁寧濃度在懷孕過程中的變化趨勢：全部樣本【圖4、圖 5】

由GEE 分析而得的結果將以兩種圖形呈現：（一）將各個懷孕週數組別的可丁寧 濃度平均值標出，並且以直線相連成折線圖；（二）將各個懷孕週數組別與基準組可 丁寧濃度平均差異做成誤差長條圖。圖4 中的基準組為懷孕週數小於 12 週者，分組 方式為每8 週一組以及每 4 週一組；圖 5 中的基準組為懷孕週數小於 10 週者，分組 方式為每10 週一組以及每 5 週一組。

由這些圖可以看出可丁寧的濃度確實有隨著懷孕週數增加而上升的趨勢，然而這 樣的趨勢不完全是直線的，特別是懷孕週數的分組較細時，在懷孕週數16 至 24 週之 間有明顯轉折之處，而在懷孕週數24 週以後又急遽升高；就尿液可丁寧和血液可丁 寧相比，同一時期的尿液可丁寧濃度平均值都比血液可丁寧來得高，並且兩者隨著懷 孕週數增加而升高的趨勢相當類似。

3.三期可丁寧濃度完整收案者之 GEE 分析

前面所做可丁寧濃度與懷孕週數的GEE 分析，乃是將全體孕婦的資料均放入模

式中，原則上每位孕婦應該都要有三次的尿液以及血液可丁寧的量測，但是由於追蹤 漏失、進入研究時為懷孕中後期以及拒絕接受採樣的緣故，並非每一位孕婦都完整測 得三次的尿液以及血液可丁寧濃度。因此在前面GEE 分析中，可丁寧濃度的資料有 蠻多的遺漏值，可能會對分析的結果造成一些偏差。

為了讓可丁寧濃度與懷孕週數的GEE 分析有較完整的資料，在此特別針對三個 孕期都測得可丁寧濃度的孕婦做進一步的分析。三期皆測得尿液可丁寧濃度的孕婦共 有130 名，但是有 1 名孕婦沒有生產期懷孕週數的資料，所以可供分析的資料共 389 筆；三期皆測得血液可丁寧濃度的孕婦共有108 名，但是有 2 名孕婦沒有生產期懷孕 週數的資料，故可供分析的資料只有322 筆。雖然可供分析的資料大幅減少，但是每 一位孕婦三期可丁寧濃度資料的完整性卻是大為提高。

可丁寧濃度在各懷孕週數組別平均值之GEE 估計：三期完整收案【表 13、表 14】

在完整收案者的GEE 分析中，只用了兩種尺度來切分懷孕週數，這兩種尺度都 是以懷孕週數未滿12 週者為基準組，而分組方式分別為每 8 週一組以及每 4 週一組，

這些懷孕週數組別的可丁寧濃度平均值都和懷孕週數未滿12 週的基準組做比較。

就尿液可丁寧來看，當懷孕週數分組為每8 週一組時，以懷孕週數未滿 12 週者 的可丁寧濃度平均值為基準，各組的尿液可丁寧濃度平均值都比基準組的平均值來得 高，這樣的差異有隨著懷孕週數增加而變大的趨勢，除了懷孕週數20 週以上但未滿 28 週的組別尿液可丁寧濃度平均值和基準組差異不顯著外，其他三組的尿液可丁寧 平均值和基準組都有相當顯著的差異。若和前面尿液可丁寧濃度GEE 分析的結果比 較（表11），由尿液可丁寧三期完整收案的資料進行 GEE 分析所得的參數估計值是 比較大的，除了懷孕週數未滿12 週之組別尿液可丁寧濃度平均值比之前的結果來得 低外，其他懷孕週數組別的尿液可丁寧濃度平均值都比之前的結果高，特別是懷孕週 數36 週以上的組別，尿液可丁寧濃度平均值高達 4.033 ng/mL。當懷孕週數分組為每 4 週一組時，尿液可丁寧濃度隨著懷孕週數增加而提高的趨勢就有了一些變化，尿液 可丁寧濃度平均值在懷孕週數16 週至 24 週之間以及 28 週至 32 週之間有降低的趨 勢，其中懷孕週數16 週至 28 週這一段，三個懷孕週數組別的尿液可丁寧濃度平均值

和基準組相比都沒有顯著差異。

就血液可丁寧來看，當懷孕週數分組方式為每8 週一組時，若和懷孕週數未滿 12 週者相比，其他各組的血液可丁寧濃度平均值都比基準組高，並且有隨著懷孕週 數增加而可丁寧濃度平均值升高的趨勢。其中只有兩組（懷孕週數12 週以上未滿 20 週、20 週以上未滿 28 週）血液可丁寧濃度平均值與基準組的差異沒有統計上顯著意 義，懷孕週數在28 週以上的兩組和基準組的差異都相當顯著。若和先前血液可丁寧 GEE 分析的結果相比（表 11），除了懷孕週數 20 週以上但未滿 28 週的組別外，其他 各組的血液可丁寧濃度平均值都比先前的數值高。當懷孕週數分組方式為每4 週一組 時，可以發現血液可丁寧濃度平均值在20 週至 24 週之間以及 40 週以上這兩個時期 有小幅的降低，其他時期都呈現可丁寧濃度平均值隨著懷孕週數增加而升高的趨勢。

其中有兩組（懷孕週數12 週以上未滿 16 週、24 週以上未滿 28 週）血液可丁寧濃度 平均值和基準組沒有統計上顯著差異，其他各個組別的可丁寧濃度平均值和基準組的 差異都相當顯著。

尿液、血液可丁寧濃度在懷孕過程中的變化趨勢：三期完整收案【圖6、圖 7】

尿液、血液可丁寧濃度在懷孕期間的變化趨勢，可以更清楚從圖中觀察出來。當 懷孕週數分組的尺度較大時（圖6[A]），可以看出尿液可丁寧和血液可丁寧濃度平均 值都很明顯的隨著懷孕週數增加而升高，特別是懷孕後期增加的幅度最大。若和之前 用全體孕婦的資料所畫出的趨勢圖（圖4[A]）相比，用完整收案資料所畫出的趨勢 圖變化的幅度較大，特別是尿液可丁寧濃度的趨勢。

當懷孕週數分組的尺度變小時，就可以看出可丁寧濃度並不是隨著懷孕週數而線 性遞增的，尿液可丁寧濃度平均值在懷孕週數28 週之前有一個大幅降低的趨勢，也 可以發現懷孕週數24 週以上但未滿 28 週這個組別的可丁寧濃度的變異非常大（圖 6[B]）。而血液可丁寧濃度平均值雖然在兩個時期有小幅的降低，但是整體而言仍有 隨著懷孕週數而增加的趨勢（圖6[D]）。

若可丁寧濃度隨著懷孕週數的趨勢能依照不同的菸害暴露狀態呈現，則可以提供

若可丁寧濃度隨著懷孕週數的趨勢能依照不同的菸害暴露狀態呈現，則可以提供