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

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

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

3、基本步驟？,2024/2/27,公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系,6,2024/2/27,公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系,7,電腦實(shí)驗(yàn)1 隨機(jī)區(qū)組設(shè)計(jì)資料方差分析（假設(shè)資料符合方差分析的前提條件）（數(shù)據(jù)文件：討論1-隨機(jī)區(qū)組.sav）,,2024/2/27,公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系,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）,,,,,,,,隨機(jī)區(qū)組設(shè)計(jì)資料為什么不能采用完全隨機(jī)設(shè)計(jì)資料的方差分析？（假設(shè)資料滿足方差分析的前提條件）,

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

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

8、024/2/27,公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系,14,電腦實(shí)驗(yàn)2 隨機(jī)區(qū)組設(shè)計(jì)資料秩和檢驗(yàn)（數(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,若要進(jìn)行兩兩比較： 方法一：Bonfferoni思想（符號(hào)秩和檢驗(yàn) ） 方法二：編秩后進(jìn)行方差分析,,,,,,,,2024/2/27,公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系,15,討論二,為研究甲乙兩藥治療高膽固醇血癥的療效（膽固醇降低值mmol/L）,甲藥有用 和不用兩個(gè)水

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

11、、乙藥的主效應(yīng)與交互效應(yīng)。,2024/2/27,公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系,18,,問題：,1. 什么是析因設(shè)計(jì)？有何特點(diǎn)？析因設(shè)計(jì)方差分析的基本思想？ 如何理解單獨(dú)效應(yīng)、主效應(yīng)、交互效應(yīng)？析因設(shè)計(jì)方差分析的基本步驟？,2024/2/27,公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系,19,2024/2/27,公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系,20,電腦實(shí)驗(yàn)3 析因設(shè)計(jì)資料方差分析 （數(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)計(jì)分析時(shí)，如果存在交互效應(yīng)，解釋結(jié)果時(shí)，要逐一分析各因素的單獨(dú)效應(yīng)。如果不存在交互效應(yīng)，則只需要考慮各因素的主效應(yīng)，報(bào)告主效應(yīng)及其95%CI。,2024/2/27,公共衛(wèi)生學(xué)院 醫(yī)學(xué)統(tǒng)計(jì)與流行病學(xué)系,21,課堂練習(xí),隨機(jī)區(qū)組方差分析例16-1（教材p304）數(shù)據(jù)文件data16-1.sav析因設(shè)計(jì)例16-6（教材p310）數(shù)據(jù)文件data16-6.sav,2

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

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

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