4.6 析因設計一元定量資料方差分析

4.6.1 問題與數據

【例4-6】某研究者欲研究IL-11藥對5.5Gy照射小鼠骨髓造血細胞周期（G0/G1期）的影響，選取45只小鼠并將其完全隨機地均分成3組，每組15只，分別在5.5Gy劑量照射前給IL-11藥、照射后給IL-11藥和照射對照（即不給IL-11藥）。每一組中的15只小鼠隨機地等分成3組，分別在照后6h、12h、24h三個時間點上處死，測量其骨髓造血細胞周期（G0/G1期）。實驗結果如表4-6所示，請進行相應的統計分析。

4.6.2 對數據結構的分析

實驗中涉及兩個實驗因素：“給藥時機”和“處死時間”，受試對象為小鼠。從表面上看，此設計像具有一個重復測量的兩因素設計，這樣看是把“處死時間”當成了重復測量因素。其實，在不同時間點上獲得的數據沒有關聯，不是來自于同一只小鼠，因為每只小鼠身上只能測到一個數據。因此，此資料不是具有一個重復測量的兩因素設計定量資料。由于實驗條件為兩因素各水平的全面組合，研究者在每種實驗條件下均進行了5次獨立重復實驗，在沒有專業依據認為這兩個實驗因素對觀測結果的影響重要程度不同的情況下，此資料應為兩因素析因設計定量資料。

4.6.3 分析目的與統計分析方法的選擇

對本資料，應采用兩因素析因設計定量資料方差分析來處理。這樣既可分析兩實驗因素的主效應，也可分析二者之間的交互效應。

4.6.4 SAS程序

SAS程序名為SASTJFX4_6.SAS。

程序第1步為建立數據集，treat代表“給藥時機”，rep代表“受試對象個體”，time代表“處死時間”，outcome代表“骨髓造血細胞周期（G0/G1期）（單位）”。第2步為進行兩因素析因設計定量資料方差分析，其中“treat|time”等價于“treat time treat*time”，TDIFF和PDIFF則分別用來輸出兩兩比較的t值和相應的P值。

4.6.5 主要分析結果及解釋

這是輸出結果的第1部分。由方差分析結果可以看出，兩個實驗因素“給藥時機”、“處死時間”及其交互作用均有統計學意義，即不同給藥時機、不同處死時間上骨髓造血細胞周期的值是不完全相同的。

這是輸出結果的第2部分。本部分首先給出了各實驗組在兩兩比較中的編號，見“LSMEAN Number”列。隨后給出了各實驗組兩兩比較的結果，行與列的交叉處即是這兩組比較的結果，上行為t值，下行為對應的P值，形式上看兩兩比較用的是t檢驗，實際上是方差分析計算出來的。另外，由于篇幅排版限制，對輸出結果進行了四舍五入，所以本部分結果比SAS實際輸出結果保留的小數點后位數要少一些。由兩兩比較結果可知：除實驗組1與實驗組7、實驗組2與實驗組8、實驗組3與實驗組6、實驗組3與實驗組9、實驗組5與實驗組8、實驗組6與實驗組9之間觀測指標的差異無統計學意義外，其他各組間觀測指標的差異均有統計學意義。