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

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

2、022-3-20公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系4討論一討論一 30名名頭痛患者頭痛患者按照病情嚴(yán)重程度分為按照病情嚴(yán)重程度分為10個(gè)個(gè)組，每一組內(nèi)組，每一組內(nèi)3名患者隨機(jī)接受名患者隨機(jī)接受標(biāo)準(zhǔn)方法標(biāo)準(zhǔn)方法或或兩兩種新方法種新方法治療，患者服用藥物后止痛時(shí)間如治療，患者服用藥物后止痛時(shí)間如表表16-1。 2022-3-20公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系5組別號(hào)組別號(hào)止痛時(shí)間止痛時(shí)間標(biāo)準(zhǔn)方法標(biāo)準(zhǔn)方法新方法一新方法一新方法二新方法二18.406.906.8027.706.806.4039.609.409.204 10.10 10.30 10.6059.308.007.6069.108.808

3、.9079.006.106.4087.707.407.3098.108.008.20105.305.105.20均均 數(shù)數(shù)8.437.687.66標(biāo)準(zhǔn)差標(biāo)準(zhǔn)差1.361.561.50表表16-1 三種方法治療頭痛患者的止痛時(shí)間（三種方法治療頭痛患者的止痛時(shí)間（min）問(wèn)題：?jiǎn)栴}：1. 什么是隨機(jī)區(qū)組設(shè)計(jì)？有何特點(diǎn)？什么是隨機(jī)區(qū)組設(shè)計(jì)？有何特點(diǎn)？2. 隨機(jī)區(qū)組設(shè)計(jì)方差分析的基本思想？隨機(jī)區(qū)組設(shè)計(jì)方差分析的基本思想？ 3. 隨機(jī)區(qū)組設(shè)計(jì)方差分析的基本步驟？隨機(jī)區(qū)組設(shè)計(jì)方差分析的基本步驟？2022-3-20公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系62022-3-20公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系7電腦

4、實(shí)驗(yàn)電腦實(shí)驗(yàn)1 1 隨機(jī)區(qū)組設(shè)計(jì)資料方差分析隨機(jī)區(qū)組設(shè)計(jì)資料方差分析（假設(shè)資料（假設(shè)資料符合方差分析符合方差分析的前提條件）的前提條件）（數(shù)據(jù)文件：（數(shù)據(jù)文件：討論討論1-隨機(jī)區(qū)組隨機(jī)區(qū)組.sav） Analyze General Linear Model Univariate: Time : Dependents variable Method : Fixed factors Block : Fixed factorsModel Custom Build terms: main effects Method : Model Block : Model Continue 2022-3-20公共

5、衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系8Post Hoc Method: Post Hoc Test for: Equal Variances Assumed: LSD、SNK、Bonferroni Dunnett：Control category：firstContinueSave Unstandardized (Predicted Values) Standardized (Residuals)OK繪制殘差圖：繪制殘差圖：Graphs Scatter Simple Define: Y Axis: ZRE_1, X Axis: method OK（分別將method替換為block、pre_1） 隨

6、機(jī)區(qū)組設(shè)計(jì)資料為什么不能采用完全隨機(jī)設(shè)計(jì)隨機(jī)區(qū)組設(shè)計(jì)資料為什么不能采用完全隨機(jī)設(shè)計(jì)資料的方差分析？資料的方差分析？（假設(shè)資料滿足方差分析的前提條件）（假設(shè)資料滿足方差分析的前提條件）2022-3-20公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系9思考 如果隨機(jī)區(qū)組設(shè)計(jì)資料如果隨機(jī)區(qū)組設(shè)計(jì)資料方差齊性的前提條方差齊性的前提條件不滿足件不滿足，應(yīng)該如何分析？，應(yīng)該如何分析？2022-3-20公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系10思考2022-3-20公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系11隨機(jī)區(qū)組設(shè)計(jì)資料的隨機(jī)區(qū)組設(shè)計(jì)資料的秩和秩和檢驗(yàn)檢驗(yàn)(Friedman 秩和檢驗(yàn)秩和檢驗(yàn)/M檢驗(yàn)檢驗(yàn))2022-3-20公

7、共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系12討論討論二二 按照性別體重，將大蟾蜍配成按照性別體重，將大蟾蜍配成10個(gè)區(qū)組個(gè)區(qū)組（b=10），），每個(gè)區(qū)組包括每個(gè)區(qū)組包括4只蟾蜍，隨機(jī)將其分配到只蟾蜍，隨機(jī)將其分配到4個(gè)處理組個(gè)處理組（k=4）中：分別在蟾蜍上顎黏膜處滴加）中：分別在蟾蜍上顎黏膜處滴加0.5ml不同的溶不同的溶液并保持液并保持30min。4種溶液分別為種溶液分別為pH 6.8 的磷酸鹽緩沖液的磷酸鹽緩沖液（PBS）、總皂苷（）、總皂苷（PNS）溶液、）溶液、PNS脂質(zhì)體和脂質(zhì)體。脂質(zhì)體和脂質(zhì)體。觀察記錄離體上腭黏膜纖毛運(yùn)動(dòng)觀察記錄離體上腭黏膜纖毛運(yùn)動(dòng)持續(xù)的時(shí)間持續(xù)的時(shí)間（min）。）。結(jié)

8、果如表結(jié)果如表16-6。 2022-3-20公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系13區(qū)組號(hào)區(qū)組號(hào)纖毛運(yùn)動(dòng)時(shí)間（秩次）纖毛運(yùn)動(dòng)時(shí)間（秩次）PBS溶液溶液PNS溶液溶液PNS脂質(zhì)體脂質(zhì)體脂質(zhì)體脂質(zhì)體1630（3）487（1）720（4）619（2）2621（4）387（1）601（3）567（2）3546（4）316（1）539（3）531（2）4498（4）257（1）264（2）367（3）5523（4）286（1）310（2）432（3）6531（4）367（1）431（3）422（2）7520（4）345（1）492（3）489（2）8532（4）324（2）335（3）316（1）9623

9、（4）321（1）620（3）611（2）10664（4）432（1）656（3）597（2）秩和秩和39112921平均秩次平均秩次3.91.12.92.1表表16-6 蟾蜍上顎黏膜纖毛運(yùn)動(dòng)時(shí)間（蟾蜍上顎黏膜纖毛運(yùn)動(dòng)時(shí)間（min）2022-3-20公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系14電腦實(shí)驗(yàn)電腦實(shí)驗(yàn)2 隨機(jī)區(qū)組設(shè)計(jì)資料秩和檢驗(yàn)隨機(jī)區(qū)組設(shè)計(jì)資料秩和檢驗(yàn)（數(shù)據(jù)文件：（數(shù)據(jù)文件：data16-4（例子（例子16-4）.sav）Analyze Nonparametric Test K Related Samples Test Variables: treat1-4 Test type: Fried

10、man Exact ：exact Continue OK若要進(jìn)行兩兩比較：若要進(jìn)行兩兩比較： 方法一：方法一：Bonfferoni思想（思想（符號(hào)秩和檢驗(yàn)符號(hào)秩和檢驗(yàn) ） 方法二：方法二：編秩編秩后進(jìn)行方差分析后進(jìn)行方差分析2022-3-20公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系15 討論討論二二 為研究為研究甲乙兩藥甲乙兩藥治療高膽固醇血癥的療效治療高膽固醇血癥的療效（膽固醇降低值（膽固醇降低值mmol/L）,甲藥有甲藥有用用 和和不用不用兩兩個(gè)水平，乙藥也有個(gè)水平，乙藥也有用用和和不用不用兩個(gè)水平。欲研究?jī)蓚€(gè)水平。欲研究甲和乙藥是否有降低膽固醇的作用。數(shù)據(jù)整理甲和乙藥是否有降低膽固醇的作用。數(shù)

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