1、第四章、多組資料均數的比第四章、多組資料均數的比較較方差分析方差分析方差分析的用途：方差分析的用途： 用于對比研究中用于對比研究中,組間（比較組組間（比較組k 2)的的樣本均數樣本均數間差別有無統計學意義間差別有無統計學意義,從而推從而推論各樣本均數的總體水平是否相同。論各樣本均數的總體水平是否相同。例講義例講義4-2：研究目的：研究目的：了解燙傷后不同時期切痂對肝臟的了解燙傷后不同時期切痂對肝臟的ATP含量的影響含量的影響 處理因素處理因素 甲組（甲組（n=10） 燙傷（對照）燙傷（對照） 30大鼠大鼠 乙組（乙組（n=10） 燙傷（燙傷（24小時切痂）小時切痂） 丙組（丙組（n =10）

2、燙傷（燙傷（96小時切痂）小時切痂） 全部在全部在 168小時處死，小時處死， 測定三組大鼠肝臟的測定三組大鼠肝臟的ATP含量，問：不同處理因素（時間）的含量，問：不同處理因素（時間）的ATP含含量是否有差別？量是否有差別？隨機隨機例：大鼠燙傷后不同時間切痂對肝臟例：大鼠燙傷后不同時間切痂對肝臟ATP含量的測量結果含量的測量結果 處理因素處理因素 對照組對照組 24小時組小時組 96小時組小時組 7.76 11.14 10.85 7.71 11.60 8.58 8.43 11.42 7.19 8.47 11.85 9.36 10.30 13.53 9.59 . . . 6.97 17.72 8

3、.69 n 10 10 10 8.04 12.76 9.25x1010N =VAR000022.001.00VAR00001201816141210864177A組B組ATP含量第一節、方差分析的基本思想第一節、方差分析的基本思想 將比較組數據變異（總變異）分解為：將比較組數據變異（總變異）分解為：總變異總變異=組間變異組間變異+組內變異組內變異 變異變異 原原 因因 方差估計量方差估計量組間變異組間變異 不同處理因素不同處理因素 MS組間組間 +隨機誤差隨機誤差組內變異組內變異 隨機誤差隨機誤差 MS誤差誤差 （隨機誤差（隨機誤差=個體變異個體變異+未知因素）未知因素）組間離均差平方和組間離

4、均差平方和kiiixxn12)(SS組間組間 =組內離均差平方和組內離均差平方和SS組內組內 =總離均差平方和總離均差平方和 kinjijixx112)(SS總總 =變異估計量（變異估計量（離均差平方和）離均差平方和）及分解及分解：組間組間組間SSMS 組內組內組內SSMS 組間組間組間/SSMS組內組內組內/SSMS其中：其中：組間組間=組數組數-1=K-1組內組內=總例數總例數-組數組數=N-K組內組間總SSSSSS注：注：MS=S2 （方差）（方差） 理論上，如果處理因素無作用，則理論上，如果處理因素無作用，則（組間變異（組間變異=組內變異）組內變異） ， F 或波動在或波動在1左右。左

5、右。 如果如果F ，或，或FF ，P ，說明處理因素有作用。說明處理因素有作用。方差分析：統計檢驗：方差分析：統計檢驗：假設：假設：H0： 1= 2= 3=. k H1：各：各不等不等 ,研究總研究總體對象體對象干預前干預前A組組B組組C組組干預后干預后不切痂不切痂24小時切痂小時切痂96小時切痂小時切痂ATP含量含量04. 8 x76.12 x25. 9 xMS組間=MS組內MS組間=MS組內 ？圖圖4-1組間誤差與組內誤差示意圖組間誤差與組內誤差示意圖第二節、完全隨機設計的方差分析第二節、完全隨機設計的方差分析 一、一、完全隨機分組設計（例：表完全隨機分組設計（例：表4-1資料）資料） 處

6、理 指標 N個個 隨機化分組隨機化分組 甲處理（甲處理（n1） （試驗對象）（試驗對象） 乙處理乙處理 （n2） 丙處理（丙處理（n3） 優點：設計簡單，各組例數可等和不相等優點：設計簡單，各組例數可等和不相等xxx二、數據結果的方差分析步驟二、數據結果的方差分析步驟 完全隨機分組設計數據的變異分解完全隨機分組設計數據的變異分解 總變異分解為：總變異分解為： 組間變異（不同處理組間變異（不同處理+誤差）誤差） 組內變異（誤差）組內變異（誤差）完全隨機分組設計數據的完全隨機分組設計數據的SS計算公式計算公式 變異來源變異來源 SS =離均差平方和離均差平方和 總變異總變異 處理間處理間 組內組內

7、 (誤差）誤差） CXxxSSijnjki 2）（總CnxkikiiiSSinjxxn 21)(11)(組間組間總組內SSSSSS NXC/)(2校正數校正數C:55頁公式頁公式4-1（總變異）、（總變異）、4-2（組間）、（組間）、4-3（誤）（誤）表表4-2 完全隨機分組設計數據的完全隨機分組設計數據的方差分析計算方差分析計算 變異變異 SS MS F值值 總總 N-1 處理處理 k-1 組內組內 N-k (誤差）誤差） 總SS組間SS組內SS組間組間SSMS 組內組內SSMS 組內組間MSMS表表4-1大鼠燙傷后不同時間切痂對大鼠燙傷后不同時間切痂對肝臟肝臟ATP含量的測量結果含量的測量

8、結果 對照組對照組 24小時組小時組 96小時組小時組 7.76 11.14 10.85 7.71 11.60 8.58 8.43 11.42 7.19 . . . 6.97 17.72 8.69 合計合計 n 10 10 10 30 8.04 12.76 9.25 10.02 80.43 127.55 92.49 300.47 676.32 1696.96 868.93 3242.21xx 2x 完全隨機分組的方差分析計算步驟完全隨機分組的方差分析計算步驟 1.數據準備（計算表數據準備（計算表4-1） 2.建立假設建立假設 H0: 1= 2= 3，即不同時期即不同時期ATP含量含量的總體均數

9、相等的總體均數相等 H1: 各時期該指標均數不等或不全等各時期該指標均數不等或不全等 =0.05 3.計算各離均差平方和（計算各離均差平方和（SS）2）計算組間離均差平方和）計算組間離均差平方和3）計算組內離均差平方和）計算組內離均差平方和1）計算總離均差平方和）計算總離均差平方和 SS估計量的計算方法和公式：估計量的計算方法和公式：802.232407.300921.32422 CXSS總83.1191049.9255.12743.80)(22221 CCnxSSiiki組間97.112831.119802.232 組間總組內SSSSSSNXC/)(2 表4-2 方差分析表 變異來源變異來源

10、 SS MS F P 組間組間 119.83 2 59.916 14.32 0.05 組內組內 112.97 27 4.184 總總 232.8026 29 5.求求P值和結論值和結論 本例本例 v組間組間=2, v誤差誤差=27 F0.05,2,27=3.35 F=14.323.55 , P0.05。 結論結論:在在 =0.05水平水平,p0.05，拒絕，拒絕H0，三組，三組ATP含量差異有統計意義，含量差異有統計意義，不同時間切痂對肝臟不同時間切痂對肝臟ATP含量不同含量不同.4.列方差分析表列方差分析表第三節、隨機單位組設計的方差分析第三節、隨機單位組設計的方差分析 隨機單位組設計隨機單

11、位組設計 又稱隨機區組、配伍組設計又稱隨機區組、配伍組設計 （the randomized block design）一、隨機單位組設計的方法一、隨機單位組設計的方法 設計：設計： 將性質大致相同的實驗對象組成一將性質大致相同的實驗對象組成一個單位組（區組），每個單位組內的個單位組（區組），每個單位組內的實驗對象數實驗對象數=處理個數。共有處理個數。共有b個單位個單位組。組。 單位組是配對設計的擴大。單位組是配對設計的擴大。例：研究三種劑量雌激素對大鼠子宮重例：研究三種劑量雌激素對大鼠子宮重量的影響量的影響 處理因素：三種劑量（處理因素：三種劑量（0.2、0.4、0.8（ug） 研究指標：不同

12、處理后的子宮重量研究指標：不同處理后的子宮重量 實驗對象及例數：大鼠實驗對象及例數：大鼠4個種系共個種系共12只只 控制因素：處理前不同種系大鼠子宮重量控制因素：處理前不同種系大鼠子宮重量 實驗設計：隨機區組設計實驗設計：隨機區組設計 方法：方法：1.將同種系的將同種系的3只大鼠為一個單位組，只大鼠為一個單位組，共共4個單位組。個單位組。 2.在單位組內隨機分配處理因素于大鼠。在單位組內隨機分配處理因素于大鼠。三個處理組的隨機區組設計方法三個處理組的隨機區組設計方法 . 單位組單位組 實驗單位（小鼠編號）實驗單位（小鼠編號） (區組區組) 1 2 3 1 甲甲 丙丙 乙乙 2 丙丙 乙乙 甲甲

13、 3 乙乙 甲甲 丙丙 4 丙丙 乙乙 甲甲實驗單位的隨機化分組在區實驗單位的隨機化分組在區組內進行。組內進行。甲甲=0.2乙乙=0.4丙丙=0.8隨機單位設計隨機單位設計 表表4-4 大鼠注射不同劑量雌激素后大鼠注射不同劑量雌激素后 子宮重量數據子宮重量數據 種類種類 雌激素劑量（雌激素劑量（ug100g -1） 0.2 0.4 0.8 A 106 116 145 B 42 68 115 C 70 111 133 D 42 63 87 隨機區組隨機區組(單位組）設計的方差分析單位組）設計的方差分析 有兩種類型資料可采用：有兩種類型資料可采用： 1.隨機單位設計實驗（單位組內的實驗隨機單位設計

14、實驗（單位組內的實驗對象為異體）的幾組數據比較。對象為異體）的幾組數據比較。 2單位組內的實驗對象為同一個體，非單位組內的實驗對象為同一個體，非隨機分組的數據。隨機分組的數據。 （例：（例：放置不同時間放置不同時間血糖濃度的結果血糖濃度的結果）例：例：8名受試者放置不同時間血名受試者放置不同時間血糖濃度的結果糖濃度的結果 受試者受試者 放置時間（分）放置時間（分） 編號編號 0 45 90 135 1 5.27 5.27 4.94 4.61 2 5.88 5.83 5.38 5.00 3 5.44 5.38 5.27 5.00 8 5.21 5.11 5.00 4.44二、隨機區組設計試驗數據

15、的二、隨機區組設計試驗數據的方差分析的原理方差分析的原理 （一）總變異分解（一）總變異分解 組間變異組間變異(處理因素處理因素+隨機誤差隨機誤差) 單位組變異單位組變異(個體間變異個體間變異+隨機誤差）隨機誤差） 誤差變異（隨機因素）誤差變異（隨機因素） 組內組內變異變異 組間變異（組間變異（MSB） 誤差變異（誤差變異（MSE）F=設計優點：將個體變異從組內變異中分離，設計優點：將個體變異從組內變異中分離，比完全隨機分組設計更容易查出處理的差別比完全隨機分組設計更容易查出處理的差別。表表4-4 大鼠注射不同劑量雌激素后大鼠注射不同劑量雌激素后子宮重量子宮重量 種類種類 雌激素劑量（雌激素劑量

16、（ug100g -1） 0.2 0.4 0.8 合計合計 A 106 116 145 367 B 42 68 115 225 C 70 111 133 314 D 42 63 87 192 65.0 89.5 120.0 91.5 260 358 480 1098 19664 34370 59508 113542 xx 2x 隨機區組設計方差分析計算步驟隨機區組設計方差分析計算步驟 例：表例：表4-4資料（數據準備）資料（數據準備） 1.建立檢驗假設建立檢驗假設 H0: 1= 2= 3 即不同劑量雌激素對子即不同劑量雌激素對子宮重量的總體均數無影響宮重量的總體均數無影響 H1: 不同劑量的雌激

17、素對子宮重量均數不同劑量的雌激素對子宮重量均數不等或不全等不等或不全等 =0.012.隨機區組隨機區組SS估計量的計算公式和步驟估計量的計算公式和步驟: 1307512/109811354222 CXSS總1）總離均差平方和）總離均差平方和2）組間（處理）離均差平方和）組間（處理）離均差平方和60744480358260)(22221 CCnXSSiki處理自由度自由度 v總=N-1=12-1V處理處理=K-1=3-14）誤差離均差平方和）誤差離均差平方和67.645767.6457607413075 區組處理總SSSSSSSSE3）區組（單位組）間離均差平方和）區組（單位組）間離均差平方和6

18、7.64573192314225367)(2222211 CCkXSSkinj區組V單位單位= n 1表表4-6 方差分析表方差分析表 變異來源變異來源 SS MS F P 總總 13075 11 處理組間處理組間 6074 2 3037 33.54 0.01 單位組間單位組間 6457 3 2152 23.77 0.01 誤差誤差 543.33 6 90.55 誤差誤差=總總-處理組處理組-單位組單位組 結論結論:在在 =0.01水平上，水平上，P0.01),拒絕拒絕H0假假設，即不同劑量的雌激素對子宮重量不同。設，即不同劑量的雌激素對子宮重量不同。78. 9,92.106, 3 ,01.

19、06, 2,01. 0 FFFF單位處理第四節、均數間的多重比較第四節、均數間的多重比較 一、簡介一、簡介 當當FF，P0.05時時，應做均數間的多重比較。做均數間的多重比較。 注意：多個均數的組間比較采用注意：多個均數的組間比較采用t檢驗會增加檢驗會增加，型誤差的概率，多重比較的統計方法可避免。型誤差的概率，多重比較的統計方法可避免。 研究者的目的研究者的目的 常用多重比較方法常用多重比較方法 1.處理組兩兩之間處理組兩兩之間 SNK-q檢驗檢驗 的比較的比較 2.各處理組僅與對照組各處理組僅與對照組 Dunnett-t檢驗檢驗 做比做比較較1. SNK-q檢驗檢驗(Student-Newm

20、an-keuls） 用于均數間的兩兩組間比較（用于均數間的兩兩組間比較（85頁）頁）)11(2BAXXnnMSSBA誤差05. 0,05. 0 pqq查查q界界值表值表BAXXBASXXq公式公式4-5檢驗檢驗公式公式如如例：大鼠燙傷后不同時間切痂對肝臟例：大鼠燙傷后不同時間切痂對肝臟ATP含量的測量結果含量的測量結果 方法與步驟（見講義方法與步驟（見講義59頁）：頁）： 1）建立檢驗假設）建立檢驗假設 H0： A= B任意比較兩組總體均數相等，任意比較兩組總體均數相等， H1： , 任意兩組任意兩組 A B ， =0.05 2）將均數依大小排列（確定）將均數依大小排列（確定a） 組別組別 對

21、照對照 96小時小時 24小時小時 均數均數 8.04 9.25 12.76 排序排序 1 2 3 a為對比組所為對比組所包含的組數包含的組數例：大鼠燙傷后不同時間切痂對肝臟例：大鼠燙傷后不同時間切痂對肝臟ATP含量測量均數的兩兩比較含量測量均數的兩兩比較 比較組比較組 a q q0.05 q0.01 P 對與對與96 -1.21 2 1.87 2.92 3.96 0.05 對與對與24 -4.72 3 7.29 3.53 4.54 0.01 96與與24 -3.51 2 5.42 2.92 3.96 0.01 3）列均數兩兩比較表，計算）列均數兩兩比較表，計算q值值 4）結論：）結論：在在 =0.05水準上，對照組與水準上，對照組與96小小時的時的ATP含量無差別，含量無差別，24小時與對照和小時與對照和96小時切痂的小時切痂的ATP含量差別有統計意義。含量差別有統計意義。24小時切痂