《chap5方差分析》PPT课件.ppt

2019/9/19,1,Chap5方差分析（analysis of variance，ANOVA）P6884 6学时 学目的与要求 ：14, 58, 掌握：基本思想 、单因素完全随机设计与配伍设计F检验 、及两两比较。 熟悉：用改变水平法进行两两比较。 了解：组内分组、析因、拉丁方设计F检验 。 教学内容提要 ： 重点讲解：基本思想 、单因素完全随机与配伍设计F检验 、及两两比较。 讲解：两两比较的思想和方法。 介绍：组内分组、析因、拉丁方设计F检验 。 重点：单因素完全随机与配伍设计F检验 。 难点：组内分组、析因、拉丁方设计F检验 。,2019/9/19,3,1 ANOVA基本思想 P68 方差分析（analysis of variance）即变异数分析：变异的来源及大小。 用途：多样本均数比较；因素间交互作用（狐假虎威，201109陈光标演唱会）；回归方程线性检验等。,特点：不直接比较平均数，而是利用变异的关系来进行判别。 基本思想: H0:1=2=n ， SS总多个组分，以随机误差为基础，按F分布规律推断。,2019/9/19,4,P71 【例5-1】 不同中药对小白鼠E-SFC的影响 3种变异（离均差平方和）： 总变异（total variation SS总） ：39个数据不等。 组内变异（within group variation SS组内）：个体、观测误差。 组间变异（between groups variation SS组间）：处理、随机误差。,2019/9/19,5,完全随机设计：独立性（independency）、正态性（normality）和方差齐性（homogeneity of variance）。 均衡设计（各n相等）：正态，s2max/s2min3可按方差齐性处理。 无重复数据的资料：配伍组设计、交叉设计、正交设计、拉丁方设计。因每个单元格子中只有1个观察数据，不需考虑正态性和方差齐性。,2019/9/19,6,2 完全随机设计资料 P69 一、完全随机设计 两或多n，独立性、正态性和方差齐性单因素方差分析（one-way ANOVA）。 1. 建立H0:1=2=n 2. 计算F ：FMS组间 / MS组内 （5-6） 3.确定P： F ,P 。 4.判断结果：P 不拒绝H0 结束。 P 拒绝H0 多重比较 。 one-way ANOVA是两样本均数t 检验的推广，两组比较时是等价的。,2019/9/19,7,二、多重比较（multiple comparison） k2，t test 两两比较第1类错误的概率 无差别判为有差别。 例：4个样本均数间的两两比较有42=4！/2！（4-2）！6次对比，若每次比较的检验水准=0.05，则6次比较均不犯第一类错误的概率为（1-0.05）6，第一类错误总概率1(10.05)60.26490.05。 t test还无法研究实验因素之间的交互作用。,2019/9/19,8,1.满足正态性和方差齐性 （1）多实验组与1对照组比较：Dunnett法。 （2）每两个均数比较： LSD（最小显著差异法） SNK（Student-Newman-Keuls即q检验）法 Tukey（可靠显著差异）法 Bonferroni调整法（校正最小显著差异）：最保守 ， 调整=比较次数m，组K7时不能用。 mK（K1）/2=21, 犯错误P0.05*211.05不符合逻辑。,2019/9/19,9,2.不满足正态性和方差齐性 数据变换方差分析。 非参数检验法：秩和检验。 近似检验：Tamhanes T2，Dunnetts T3，Games-Howell，Dunnetts C等。,2019/9/19,10,P71 【例5-1】 研究单味中药对小白鼠细胞免疫机能的影响，39只小白鼠随机分为4组，雌雄尽量各半，用药15天后，进行E-玫瑰花结形成率（E-SFC）测定，结果如表5-1。 表5-1 不同中药对小白鼠E-SFC的影响,2019/9/19,11,【SPSS操作】以药物和e为变量名，建立39行2列的数据文件L5-2.sav： （1）正态性与方差齐性检验 Analyze Descriptive Statistics explore e 入Dependent variable,药物入Factor plots Nomality plots with tests UntransformedOK,2019/9/19,12,2019/9/19,13,2019/9/19,14,（2）单因素方差分析 Analyze Compare Means One-way ANOVA E入Dependent variable，药物入Factor Post Hoc(多重比较)， LSD， Dunnett，Control Categories下拉栏选First Continue Options，Descriptive、Homogeneity of variance test ContinueOK。,2019/9/19,15,2019/9/19,16,2019/9/19,17,2019/9/19,18,3 配伍组设计资料的方差分析及多重比较 P72 一、配伍组设计资料：每格1个数据，无重复。 两因素（Two-way ANOVA）：处理因素、配伍因素。配对 t 的推广。 【例5-3】,2019/9/19,19,建立H0：处理组：不同疗法治愈天数相等； 配伍组：不同年龄治愈天数相等。 2. 计算F ： F处理MS处理/MS误差， F配伍MS配伍/MS误差 3.确定P值： F ,P 。 4.判断结果：,2019/9/19,20,【SPSS操作】以疗法（标签13）、年龄（标签15）和天数为变量名，建立3列15行的数据文件L5-3.sav Analyze General Linear Models Univariate单变量 ，天数入Dependent variable，疗法、年龄入Fixed Factor（s） Model，Custom,在Build Terms下选Main effects，年龄、疗法从Factors & Covariates入Model ContinuePost Hoc，年龄、疗法入Post Hoc Tests，LSDOK,2019/9/19,21,2019/9/19,22,2019/9/19,23,2019/9/19,24,二、完全随机设计与配伍设计方差分析的比较 1. SS 配伍设计：预估配伍因素对结果影响较大，配伍组差异可能有显著意义。 2. df 完全随机设计：不具备配伍条件，估计配伍组差异小者，n大。,2019/9/19,25,4 拉丁方设计 P76 一、拉丁方设计：三因素ANOVA 总变异处理间、区组间、序列间、误差。 【例5-4】A、B、C、D4种降脂中药，4窝大鼠各取4只体重相近的雄鼠，用 、 、 、 4种饲料，结果见表5-4。作方差分析。,表5-4 4种中药对主动脉硬化斑块面积影响（mm2） 窝别 饲料（序列） （区组） 1 C85.00 D102.00 A90.00 B100.0 2 D108.00 A88.00 B110.00 C80.00 3 A95.00 B98.00 C82.00 D112.0 4 B105.00 C78.00 D104.00 A92.00,2019/9/19,26,【SPSS操作】以斑块面积和中药（标1A，2B，3C，4D）、窝别、饲料为变量名，建立4列16行数据文件L5-4.sav。 Analyze General Linear Models Univariate，斑块面积入Dependent variable，中药、窝别、饲料入Fixed Factors Model，Custom将中药、窝别、饲料入Model，在Build Term（分析效应选项）选Main effects ContinuePost Hoc，中药、窝别、饲料入Display Means for， LSD， Bonferroni， Continue OK。,2019/9/19,27,2019/9/19,28,2019/9/19,29,2019/9/19,30,2019/9/19,31,2019/9/19,32,2019/9/19,33,2019/9/19,34,2019/9/19,35,2019/9/19,36,2019/9/19,37,二、拉丁方设计与配伍组设计比较 1. 拉丁方SS误差 F效能。 2. 拉丁方 df误差 F F效能 。 凡具备二因素实验条件，且预估第三因素（序列）间差异显著者用拉丁方设计，否则用配伍组设计。,2019/9/19,38,5 析因设计的方差分析 P77 总变异各因素、因素间的交互作用和误差。 【例5-5】 12例病人随机分为4组：第1组一般疗法，第2组一般加A药，第3组一般加B，第4组一般加A和B。1个月后观察胆固醇降低量（mg%）。A、B药有无降胆固醇作用？两药有无交互作用？,表5-5 两药降胆固醇量析因试验结果,2019/9/19,39,A、B各有用和不用2水平22析因设计 药物间：H0：A、B两药无疗效 H1：A、B两药有疗效 交互作用间：H0：A、B无交互作用 H1：A、B有交互作用,2019/9/19,40,【SPSS操作】以a、b、胆固醇为变量名，分别在a，bValues中标0=不用，1=用，建立3列12行的数据文件L5-5.sav。 Analyze General Linear Models Univariate，胆固醇入Dependent variable，a、b入Fixed Factor（s） option a、b 、 a*b display means for descriptive statistics continue OK,2019/9/19,41,2019/9/19,42,2019/9/19,43,6 交叉设计资料的方差分析 P78 两阶段交叉设计：三因素方差分析，变异处理间、阶段间、个体间与误差。 【例5-6】某中药（A）与胃泌素（B）对巴氏小胃盐酸分泌量（mol/h）的影响，用12只小狗进行交叉设计，清洗期为两天，试验结果见表5-6：,表5-6 12只小狗交叉设计试验巴氏小胃盐酸分泌量,2019/9/19,44,【SPSS操作】以狗号、阶段、药物（1中药，2胃泌素）、试验结果为变量名，建立4列24行的数据文件L5-6.sav。 Analyze General Linear Models Univariate，试验结果入Dependent variable，狗号、阶段、药物入Fixed Factors Model，选Custom，狗号、阶段、药物入Model；Build Term(s)下拉框选Main effects Continue Options，将阶段、药物入Display Means for，Compare main effects，Descriptive statistics ContinueOK。,2019/9/19,45,2019/9/19,46,2019/9/19,47,2019/9/19,48,2019/9/19,49,2019/9/19,50,2019/9/19,51,2019/9/19,52,2019/9/19,53,2019/9/19,54,2019/9/19,55,7 组内分组资料的方差分析 P79 基本思想：假设各总体均数相等。 SS总 SS大组 、SS小组、SS误 。 组内分组设计资料，不分析交互作用。 SS误 = SS总SS大组SS小组 df误=df总-df大组-df小组 MS = SS /df （5-11） F大组= MS大组MS小组 F小组= MS小组MS误差 （5-12）,2019/9/19,56,【例5-7】 ，不同产地的当归不同部位（头为1号，身为2号，尾为3号）的M物质含量（mg /10g），结果见表5-7。大组为产地，小组为部位，组内分组。 两种分析： 产地； 部位。,2019/9/19,57,【SPSS操作】以产地、部位、观测指标为变量名建立3列17行的数据文L5-7.sav。 Analyze General Linear Models Univariate ，M入Dependent Variable，产地、部位入Fixed Factors Model，Custom，按大组、小组的顺序将产地、部位逐个选入Model，Build Term(s)下选Main effects，选Type Continue Options，产地、部位入Display Means for， Compare main effects ContinueOK。,2019/9/19,58,2019/9/19,59,2019/9/19,60,2019/9/19,61,2019/9/19,62,2019/9/19,63,2019/9/19,64,2019/9/19,65,2019/9/19,66,2019/9/19,67,8 定量反应结果时样本含量的估计 P80 一、意义 n反映重复性原则：n 抽样误差。 n 实验过程控制难度系统误差。 n的估计是在保证研究结论具有一定可靠性（控制CI宽度和检验效能）的前提下的最小n。,2019/9/19,68,1. 影响样本含量的因素,2019/