隨機區(qū)組設計與析因設計資料的分析

1、2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,1,實習十二,隨機區(qū)組設計與析因設計 資料的分析,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,2,目的和要求,1. 掌握隨機區(qū)組設計的特點與數(shù)據(jù)分析方法。2. 掌握析因設計的特點與數(shù)據(jù)分析方法。3. 掌握單獨效應、主效應和交互效應的意義。4.了解多因素方差分析時的殘差圖判讀。,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,3,討 論 內(nèi) 容

2、,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,4,討論一,30名頭痛患者按照病情嚴重程度分為10個組，每一組內(nèi)3名患者隨機接受標準方法或兩種新方法治療，患者服用藥物后止痛時間如表16-1。,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,5,表16-1 三種方法治療頭痛患者的止痛時間（min）,問題：,1. 什么是隨機區(qū)組設計？有何特點？隨機區(qū)組設計方差分析的基本思想？ 3. 隨機區(qū)組設計方差分析的

3、基本步驟？,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,6,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,7,電腦實驗1 隨機區(qū)組設計資料方差分析（假設資料符合方差分析的前提條件）（數(shù)據(jù)文件：討論1-隨機區(qū)組.sav）,,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,8,Post Hoc Method: Post Hoc Test for:

4、 Equal Variances Assumed: LSD、SNK、Bonferroni Dunnett：Control category：firstContinue…Save Unstandardized (Predicted Values) Standardized (Residuals)

5、OK繪制殘差圖：Graphs Scatter Simple Define: Y Axis: ZRE_1, X Axis: method OK（分別將method替換為block、pre_1）,,,,,,,,隨機區(qū)組設計資料為什么不能采用完全隨機設計資料的方差分析？（假設資料滿足方差分析的前提條件）,

6、2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,9,思考,,如果隨機區(qū)組設計資料方差齊性的前提條件不滿足，應該如何分析？,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,10,思考,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,11,隨機區(qū)組設計資料的秩和檢驗(Friedman 秩和檢驗/M檢驗),2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,12,討論二,按照性別體重，將大蟾蜍配成

7、10個區(qū)組（b=10），每個區(qū)組包括4只蟾蜍，隨機將其分配到4個處理組（k=4）中：分別在蟾蜍上顎黏膜處滴加0.5ml不同的溶液并保持30min。4種溶液分別為pH 6.8 的磷酸鹽緩沖液（PBS）、總皂苷（PNS）溶液、PNS脂質(zhì)體和脂質(zhì)體。觀察記錄離體上腭黏膜纖毛運動持續(xù)的時間（min）。結果如表16-6。,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,13,表16-6 蟾蜍上顎黏膜纖毛運動時間（min）,,,,,2

8、024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,14,電腦實驗2 隨機區(qū)組設計資料秩和檢驗（數(shù)據(jù)文件：data16-4（例子16-4）.sav）,Analyze Nonparametric Test K Related Samples… Test Variables: treat1-4 Test type: Friedman Exact ：ex

9、act Continue OK,若要進行兩兩比較： 方法一：Bonfferoni思想（符號秩和檢驗 ） 方法二：編秩后進行方差分析,,,,,,,,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,15,討論二,為研究甲乙兩藥治療高膽固醇血癥的療效（膽固醇降低值mmol/L）,甲藥有用 和不用兩個水

10、平，乙藥也有用和不用兩個水平。欲研究甲和乙藥是否有降低膽固醇的作用。數(shù)據(jù)整理如下：,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,16,表16-6 兩種藥物不同搭配方案的膽固醇降低值（ mmol/L ）,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,17,某人對四個處理組的均數(shù)進行了完全隨機設計的單因素方差分析，發(fā)現(xiàn)差異有統(tǒng)計學意義(P<0.05)。（結果見后頁）,試點評該分析方法。,課題組希望考察甲藥

11、、乙藥的主效應與交互效應。,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,18,,問題：,1. 什么是析因設計？有何特點？析因設計方差分析的基本思想？ 如何理解單獨效應、主效應、交互效應？析因設計方差分析的基本步驟？,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,19,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,20,電腦實驗3 析因設計資料方差分析 （數(shù)據(jù)文件： 討論2-析因.s

12、av）,Analyze General Linear Model Univariate Dependents variable: 膽固醇降低值； Fixed factors: A、BModel ：Full Continue…,,,Plots ：Horizontal Axias: B； Separate Lines: A Add

13、Continue…Save ：Unstandardized (Predicted Values) Standardized (Residuals) Continue…Options：Display means for: A、B Compare main effects Display: descriptive statistics

14、 Continue…OK,,,,,,,,,,注意,統(tǒng)計分析時，如果存在交互效應，解釋結果時，要逐一分析各因素的單獨效應。如果不存在交互效應，則只需要考慮各因素的主效應，報告主效應及其95%CI。,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,21,課堂練習,隨機區(qū)組方差分析例16-1（教材p304）數(shù)據(jù)文件data16-1.sav析因設計例16-6（教材p310）數(shù)據(jù)文件data16-6.sav,2

15、024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,22,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,23,例16-6（p310）,某醫(yī)生研究四氧嘧啶（ALX）的劑量和造模前12h禁食對制作小鼠糖尿病模型影響。四氧嘧啶劑量（A因素）:150mg/kg(A1)、200mg/kg(A2)，飲食控制（B因素）：造模前12h禁食(B1)、不禁食(B2)。將40只小鼠隨機分配到四種實驗組合，每種組合10只小鼠，造模后72h測

16、量小鼠的血糖濃度（mmol/L）。,,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,24,,結論,FAB=2.334，P>0.05 ，接受H0，可認為四氧嘧啶劑量與是否禁食間無交互作用。FA=7.183，P<0.05，拒絕H0，可認為四氧嘧啶200mg/kg與150mg/kg劑量對應的總體均數(shù)不等。兩者平均血糖濃度之差的95%CI為(0.009，7.150)mmol/L。FB=6.629， P<0.0

17、5 ，拒絕H0，可認為禁食與不禁食對應的總體均數(shù)不等。兩者平均血糖濃度之差的95%CI為(0.830，6.990) mmol/L。,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,25,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系,26,作業(yè),1) p. 328. 2 (假設該資料滿足參數(shù)檢驗條件) 2) p. 329. 4,2024/2/27,公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系