随机区组设计与析因设计资料的分析

1、2022-3-20公共卫生学院 医学统计与流行病学系1实习十二实习十二 随机区组设计与析因设计随机区组设计与析因设计 资料的分析资料的分析2022-3-20公共卫生学院 医学统计与流行病学系2目的和要求目的和要求1. 掌握掌握随机区组设计随机区组设计的特点与数据分析方法。的特点与数据分析方法。2. 掌握掌握析因设计析因设计的特点与数据分析方法。的特点与数据分析方法。3. 掌握掌握单独效应、主效应单独效应、主效应和和交互效应交互效应的意义。的意义。4.了解了解多因素方差分析多因素方差分析时的时的残差图残差图判读。判读。2022-3-20公共卫生学院 医学统计与流行病学系3讨讨 论论 内内 容容2

2、022-3-20公共卫生学院 医学统计与流行病学系4讨论一讨论一 30名名头痛患者头痛患者按照病情严重程度分为按照病情严重程度分为10个个组，每一组内组，每一组内3名患者随机接受名患者随机接受标准方法标准方法或或两两种新方法种新方法治疗，患者服用药物后止痛时间如治疗，患者服用药物后止痛时间如表表16-1。 2022-3-20公共卫生学院 医学统计与流行病学系5组别号组别号止痛时间止痛时间标准方法标准方法新方法一新方法一新方法二新方法二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均均 数数8.437.687.66标准差标准差1.361.561.50表表16-1 三种方法治疗头痛患者的止痛时间（三种方法治疗头痛患者的止痛时间（min）问题：问题：1. 什么是随机区组设计？有何特点？什么是随机区组设计？有何特点？2. 随机区组设计方差分析的基本思想？随机区组设计方差分析的基本思想？ 3. 随机区组设计方差分析的基本步骤？随机区组设计方差分析的基本步骤？2022-3-20公共卫生学院 医学统计与流行病学系62022-3-20公共卫生学院 医学统计与流行病学系7电脑

4、实验电脑实验1 1 随机区组设计资料方差分析随机区组设计资料方差分析（假设资料（假设资料符合方差分析符合方差分析的前提条件）的前提条件）（数据文件：（数据文件：讨论讨论1-随机区组随机区组.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、卫生学院 医学统计与流行病学系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、机区组设计资料为什么不能采用完全随机设计随机区组设计资料为什么不能采用完全随机设计资料的方差分析？资料的方差分析？（假设资料满足方差分析的前提条件）（假设资料满足方差分析的前提条件）2022-3-20公共卫生学院 医学统计与流行病学系9思考 如果随机区组设计资料如果随机区组设计资料方差齐性的前提条方差齐性的前提条件不满足件不满足，应该如何分析？，应该如何分析？2022-3-20公共卫生学院 医学统计与流行病学系10思考2022-3-20公共卫生学院 医学统计与流行病学系11随机区组设计资料的随机区组设计资料的秩和秩和检验检验(Friedman 秩和检验秩和检验/M检验检验)2022-3-20公

7、共卫生学院 医学统计与流行病学系12讨论讨论二二 按照性别体重，将大蟾蜍配成按照性别体重，将大蟾蜍配成10个区组个区组（b=10），），每个区组包括每个区组包括4只蟾蜍，随机将其分配到只蟾蜍，随机将其分配到4个处理组个处理组（k=4）中：分别在蟾蜍上颚黏膜处滴加）中：分别在蟾蜍上颚黏膜处滴加0.5ml不同的溶不同的溶液并保持液并保持30min。4种溶液分别为种溶液分别为pH 6.8 的磷酸盐缓冲液的磷酸盐缓冲液（PBS）、总皂苷（）、总皂苷（PNS）溶液、）溶液、PNS脂质体和脂质体。脂质体和脂质体。观察记录离体上腭黏膜纤毛运动观察记录离体上腭黏膜纤毛运动持续的时间持续的时间（min）。）。结

8、果如表结果如表16-6。 2022-3-20公共卫生学院 医学统计与流行病学系13区组号区组号纤毛运动时间（秩次）纤毛运动时间（秩次）PBS溶液溶液PNS溶液溶液PNS脂质体脂质体脂质体脂质体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 蟾蜍上颚黏膜纤毛运动时间（蟾蜍上颚黏膜纤毛运动时间（min）2022-3-20公共卫生学院 医学统计与流行病学系14电脑实验电脑实验2 随机区组设计资料秩和检验随机区组设计资料秩和检验（数据文件：（数据文件：data16-4（例子（例子16-4）.sav）Analyze Nonparametric Test K Related Samples Test Variables: treat1-4 Test type: Fried

10、man Exact ：exact Continue OK若要进行两两比较：若要进行两两比较： 方法一：方法一：Bonfferoni思想（思想（符号秩和检验符号秩和检验 ） 方法二：方法二：编秩编秩后进行方差分析后进行方差分析2022-3-20公共卫生学院 医学统计与流行病学系15 讨论讨论二二 为研究为研究甲乙两药甲乙两药治疗高胆固醇血症的疗效治疗高胆固醇血症的疗效（胆固醇降低值（胆固醇降低值mmol/L）,甲药有甲药有用用 和和不用不用两两个水平，乙药也有个水平，乙药也有用用和和不用不用两个水平。欲研究两个水平。欲研究甲和乙药是否有降低胆固醇的作用。数据整理甲和乙药是否有降低胆固醇的作用。数

11、据整理如下：如下：甲药甲药(A)乙乙 药药(B) 用用(b1) 不用不用(b2)用用 (a1) 67 78 80 56 44 32 不用不用(a2) 28 30 23 17 25 18 2022-3-20公共卫生学院 医学统计与流行病学系16表表16-6 两种药物不同搭配方案的胆固醇降低值（两种药物不同搭配方案的胆固醇降低值（ mmol/L ）2022-3-20公共卫生学院 医学统计与流行病学系17 某人对四个处理组的均数进行了某人对四个处理组的均数进行了完全随完全随机设计机设计的的单因素单因素方差分析，发现差异有统方差分析，发现差异有统计学意义计学意义(P0.05 ，接受，接受H0，可认为，可认为四氧嘧四氧嘧啶剂量啶剂量与是否禁食间无交互作用。与是否禁食间无交互作用。 FA=7.183，P0.05，拒绝，拒绝H0，可认为四氧嘧啶，可认为四氧嘧啶200mg/kg与与150mg/kg剂量对应的总体均数不等剂量对应的总体均数不等。两者平均血糖浓度之差的。两者平均血糖浓度之差的95%CI为为(0.009，7.150)mmol/L。 FB=6.629， P0.05 ，拒绝拒绝H0，可认为禁食与不，可认为禁食与不禁食对应的总体均数不等。禁食对应的总体均数不等。两者平均血糖