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

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

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

3、組設(shè)計方差分析的基本思想？ 3. 隨機區(qū)組設(shè)計方差分析的基本步驟？隨機區(qū)組設(shè)計方差分析的基本步驟？2022-4-29公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系7電腦實驗電腦實驗1 1 隨機區(qū)組設(shè)計資料方差分析隨機區(qū)組設(shè)計資料方差分析（假設(shè)資料（假設(shè)資料符合方差分析符合方差分析的前提條件）的前提條件）（數(shù)據(jù)文件：（數(shù)據(jù)文件：討論討論1-隨機區(qū)組隨機區(qū)組.sav） Analyze General Linear Model Univariate: Time : Dependents variable Method : Fixed factors Block : Fixed factorsModel Cust

4、om Build terms: main effects Method : Model Block : Model Continue 2022-4-29公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系8Post Hoc Method: Post Hoc Test for: Equal Variances Assumed: LSD、SNK、Bonferroni Dunnett：Control category：firstContinueSave Unstandardized (Predicted Values) Standardized (Residuals)OK繪制殘差圖：繪制殘差圖：Graphs Scat

5、ter Simple Define: Y Axis: ZRE_1, X Axis: method OK（分別將method替換為block、pre_1）2022-4-29公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系9 隨機區(qū)組設(shè)計資料為什么不能采用完全隨機設(shè)計隨機區(qū)組設(shè)計資料為什么不能采用完全隨機設(shè)計資料的方差分析？資料的方差分析？（假設(shè)資料滿足方差分析的前提條件）（假設(shè)資料滿足方差分析的前提條件）2022-4-29公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系10 如果隨機區(qū)組設(shè)計資料如果隨機區(qū)組設(shè)計資料方差齊性的前提條方差齊性的前提條件不滿足件不滿足，應該如何分析？，應該如何分析？2022-4-29公共衛(wèi)生學院

6、 醫(yī)學統(tǒng)計與流行病學系11隨機區(qū)組設(shè)計資料的隨機區(qū)組設(shè)計資料的秩和秩和檢驗檢驗(Friedman 秩和檢驗秩和檢驗/M檢驗檢驗)2022-4-29公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系12討論討論二二 按照性別體重，將大蟾蜍配成按照性別體重，將大蟾蜍配成10個區(qū)組個區(qū)組（b=10），），每個區(qū)組包括每個區(qū)組包括4只蟾蜍，隨機將其分配到只蟾蜍，隨機將其分配到4個處理組個處理組（k=4）中：分別在蟾蜍上顎黏膜處滴加）中：分別在蟾蜍上顎黏膜處滴加0.5ml不同的溶不同的溶液并保持液并保持30min。4種溶液分別為種溶液分別為pH 6.8 的磷酸鹽緩沖液的磷酸鹽緩沖液（PBS）、總皂苷（）、總皂苷（PNS

7、）溶液、）溶液、PNS脂質(zhì)體和脂質(zhì)體。脂質(zhì)體和脂質(zhì)體。觀察記錄離體上腭黏膜纖毛運動觀察記錄離體上腭黏膜纖毛運動持續(xù)的時間持續(xù)的時間（min）。）。結(jié)果如表結(jié)果如表16-6。 2022-4-29公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系13表表16-6 蟾蜍上顎黏膜纖毛運動時間（蟾蜍上顎黏膜纖毛運動時間（min）2022-4-29公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系14電腦實驗電腦實驗2 隨機區(qū)組設(shè)計資料秩和檢驗隨機區(qū)組設(shè)計資料秩和檢驗（數(shù)據(jù)文件：（數(shù)據(jù)文件：data16-4（例子（例子16-4）.sav）Analyze Nonparametric Test K Related Samples Test

8、Variables: treat1-4 Test type: Friedman Exact ：exact Continue OK若要進行兩兩比較：若要進行兩兩比較： 方法一：方法一：Bonfferoni思想（思想（符號秩和檢驗符號秩和檢驗 ） 方法二：方法二：編秩編秩后進行方差分析后進行方差分析2022-4-29公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系15 討論二討論二 為研究為研究甲乙兩藥甲乙兩藥治療高膽固醇血癥的療效治療高膽固醇血癥的療效（膽固醇降低值（膽固醇降低值mmol/L）,甲藥有甲藥有用用 和和不用不用兩兩個水平，乙藥也有個水平，乙藥也有用用和和不用不用兩個水平。欲研究兩個水平。欲研究甲

9、和乙藥是否有降低膽固醇的作用。數(shù)據(jù)整理甲和乙藥是否有降低膽固醇的作用。數(shù)據(jù)整理如下：如下：2022-4-29公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系16表表16-6 兩種藥物不同搭配方案的膽固醇降低值（兩種藥物不同搭配方案的膽固醇降低值（ mmol/L ）2022-4-29公共衛(wèi)生學院 醫(yī)學統(tǒng)計與流行病學系17 某人對四個處理組的均數(shù)進行了某人對四個處理組的均數(shù)進行了完全隨完全隨機設(shè)計機設(shè)計的的單因素單因素方差分析，發(fā)現(xiàn)差異有統(tǒng)方差分析，發(fā)現(xiàn)差異有統(tǒng)計學意義計學意義(P0.05 ，接受，接受H0，可認為，可認為四氧嘧四氧嘧啶劑量啶劑量與是否禁食間無交互作用。與是否禁食間無交互作用。 FA=7.183，P0.05，拒絕，拒絕H0，可認為四氧嘧啶，可認為四氧嘧啶200mg/kg與與150mg/kg劑量對應的總體均數(shù)不等劑量對應的總體均數(shù)不等。兩者平均血糖濃度之差的。兩者平均血糖濃度之差的95%CI為為(0.009，7.150)mmol/L。 FB=6.629， P0.05 ，拒絕拒絕H0，可認為禁食與不，可認為禁食與不禁食對應的總體均數(shù)不等。禁食對應的總體均數(shù)不等。兩者平均血糖濃度兩者平均