南医大统计学之数据方差分析.pptVIP

方差分析(ANOVA) －－多个均数的比较 Analysis of Variance 南京医科大学 卫生统计学 问题的提出 某研究者欲探讨某一治疗晚期胃癌的新药A的治疗效果，共纳入120名受试者，随机分为三组，每组40例： A组：安慰剂组 B组：阳性对照组（经典药物，较贵） C组：试验组（新药，较便宜） 当试验药物疗效（生活质量评分 ）比安慰剂好，同时比阳性对照药差的不太多，则可批准上市。 因素和水平 因素(factors)：将试验对象随机分为若干个组，加以不同的干预，称为处理因素。 在相同的因素下的不同干预，称为不同的水平(level) 问题的提出 t检验实例 南京医科大学的学工部门想知道2008年南京市进入南医大就读的学生中，市区的学生和郊区的学生成绩是否有差异。 因素：学生所来自的区域 水平：市区，郊区 问题的提出 t检验的局限性 单因素两水平 问题的提出 t检验实例 阿卡波糖的降血糖效果。分别给对照组和试验组服用拜唐平胶囊和阿卡波糖胶囊，观察8周后血糖下降值的差异。 因素：不同的药物 水平：阿卡波糖，拜唐平 问题的提出 江苏、安徽、浙江三省的平均入学成绩 问题的提出 假如 安徽和江苏比，犯一类错误的概率是0.05， 安徽和浙江比，犯一类错误的概率是0.05， 浙江和江苏比，犯一类错误的概率是0.05。 那么这三次相关的比较，犯一类错误的概率是 问题的提出 一种新的降血脂药，120人分为安慰剂组，用药组1(2.4g)，用药组2(4.8g)，用药组3(7.2g)。实验结束后观察血脂水平。 问题的提出 假如每次t检验犯一类错误的概率是0.05，那么要完全地进行比较，犯一类错误的概率是 方差分析 方差分析，又称变异数分析。 Analysis of Variance，简写为ANOVA。 由英国统计学家R.A.Fisher提出。 例 三组血红蛋白增加量(g) 变异的分解 A组内个体间血红蛋白增加量不同 组内变异原因：个体编译 A、B、C三组间血红蛋白平均增加量不同 组间变异原因：可能是三种药物治疗效果有差别（处理因素间有本质差别）、随机误差 所有个体血红蛋白增量几乎都不相同 这些差异的总和就是总变异 总变异 = 组内变异 + 组间变异 总变异 SS总 Sum of squares about the mean of all N values. 组间变异 SS组间 Sum of squares between groups 组内变异 SS组内 Sum of squares within groups 总变异的分解 SS总＝SS组间＋SS组内 67.6685＝8.2930+59.3755 SS的本质是离均差平方和 方差分析的原理 方差分析表 方差分析表 2 方差分析的基本思想 组内变异：抽样误差 组间变异：组间本质差别＋抽样误差 如果组间无本质差别，则组间变异＝组内变异 或： RA Fisher Rothamste Agricultural Station 3 F 分布 方差比的分布！ F 分布 F 分布 F(2,57)的F分布及界值 方差分析表 方差分析假设检验步骤 Step 1: 建立假设 H0:三种药物治疗贫血效果相同 H1:三种药物治疗贫血效果 不同或不全相同 检验水准 0.05 Step 2: 计算统计量F Step 3: 计算P Step 4: 下统计学结论及专业结论 例子的结论是： 三种药物中至少有两种药物间的治疗效果不同A-B? A-C? B-C?  方差分析后的 两两比较 多重比较(multiple comparison) 多组间的两两比较为什么不能用 t 检验？ 进行一次假设检验，犯I类错误的概率：? 进行多次(k)假设检验，犯I类错误的概率： 1－(1－?)k 组数为3, k=3,　1－(1－0.05)k=0.1426 组数为4, k=6,　1－(1－0.05)k=0.2649 组数为5, k=10, 1－(1－0.05)k=0.4013 SNK法 SNK(Student-Newman-Keuls)法 用于多组间的两两比较。 SNK法步骤 H0: 相比较的两总体均数相等； H1: 相比较的两总体均数不等。 ?＝0.05。 计算检验统计量: q 组次 　1 　2 　3 均数 1.840 1.415 0.930 组别 A B C a=2 a=2 a=3 SNK法步骤 均数 1.840 1.415 0.930 组别 A B C 其它两两比较 LSD 法 其余组与 某一组 比较 Bonferroni