第九章 方差分析 讲授内容 方差分析的基本思想及应用条件 完全随机设计资料的方差分析 随机区组设计资料的方差分析 重复测量资料的方差分析

Slides:



Advertisements
Similar presentations
健康活力班 臺南市海佃國小 張雅芬護理師.
Advertisements

第6章 方差分析与试验设计 会计学2011级 主讲:王红娜.
張偉豪 三星統計服務有限公司 執行長 Amos 亞洲一哥
上海体育职业学院 祁社生 一、重视体育科研在提高竞技运动训练水平中的意义和作用
阶段总结 计量资料假设检验总结及实例分析.
The application of medical statistics methods
人群健康研究的统计方法 预防医学系 指导教师:方亚 电话:
2011年10月31日是一个令人警醒的日子,世界在10月31日迎来第70亿人口。当日凌晨,成为象征性的全球第70亿名成员之一的婴儿在菲律宾降生。 ?
《开胃理脾口服液对脾虚小鼠肠功能的影响》
第3节 体内物质的运输.
國小學童性教育介入行動研究— 以屏東縣里港國民小學五年級學生為例
第一章 建筑工程造价概述 分部分项工程量清单的编制及工程实例.
禁毒知识教育 ·.
生物統計與SAS軟體課程教學(三) 雙變項統計分析(一)
初级会计实务 第八章 产品成本核算 主讲人:杨菠.
案例: 今年7月下旬在松江中心医院住院的刘某病情危重,血小板下降到正常人的10%左右,如不立即输血小板,极可能出现颅内出血、消化道出血等,任何止血药都止不住,可血库紧缺血小板。 松江区血站第一时间向成分献血志愿者发出了短信,立刻得到很多志愿者的响应,成分献血志愿者的接力捐献挽救了病患生命 。
中考阅读 复习备考交流 西安铁一中分校 向连吾.
變異數分析 (Analysis-of-Variance簡稱ANOVA)
五至七年級數理創意營教材之研發與教學實驗 期末報告
15 簡單迴歸分析與相關分析  學習目的.
生物科 陳鳳娟老師.
幸福大讲堂 也谈老年朋友的 “老有所□” 爸妈在线专家宣讲团 ——老年朋友如何乐度后半生概述 主讲:钱锡安
预备年级体育理论 肥胖与消瘦的危害.
中央广播电视大学开放教育 成本会计(补修)期末复习
第3节 体内物质的运输.
人教版义务教育课程标准实验教科书 小学数学四年级上册第七单元《数学广角》 合理安排时间 248.
第三章 资料的统计描述 上一张 下一张 主 页 退 出.
中考语文积累 永宁县教研室 步正军 2015.9.
小学数学知识讲座 应用题.
Dr. Hamda Qotba, B.Med.Sc, M.D, ABCM 翻译:acred(DXY)
平均数检定 庄文忠 副教授 世新大学行政管理学系 SPSS之应用(庄文忠副教授) 2012/7/6.
7-1 能量的形式和轉換 1 of 12 能量是促成自然現象變化的根源,太陽能替我們將水搬到高處,人類再利用高、低水位差發電。
倒装句之其他句式.
健康體位講座 如何增進健康體適能 演講者:張書軒 老師 96.09.28.
Analysis of Variance 變異數分析
第 22 课 孙中山的民主追求 1 .近代变法救国主张的失败教训: “师夷之长技以制 夷”“中体西用”、兴办洋务、变法维新等的失败,使孙中山
传媒学院2013年度团委工作 总结分析报告
實 驗 研 究 法 多因子實驗設計 指導老師:黃萬居教授 學生:陳志鴻 m
第13章變異數分析與多變數分析  本章的學習主題 
第十六章 無母數統計.
第十六章 無母數統計 陳順宇 教授 成功大學統計系.
平均数检定 庄文忠 副教授 世新大学行政管理学系 计量分析一(庄文忠副教授) 2011/7/12.
用于计算智能研究 和人类主观测试的 统计检验 Hideyuki TAKAGI 日本 九州大学
第 13 章 實驗設計與變異數分析.
ANOVA簡介 許晉誠
第13章變異數分析與多變數分析  本章的學習主題  1. 變異數分析的應用時機 2. 變異數分析的假設前提
Test for difference among the means: t Test
統計方法的概念與應用 一、認識統計(statistics)、測驗(test)、 測量(measurement)與評價(evaluation)
本章重點: 一、集中量數的意義和種類 二、算術平均數 三、中位數(中數) 四、眾 數 五、其他集中量數 六、SPSS12.0實務操作
十、變方分析 (Analysis of Variance) (Chapter 10)
論文計畫書 國軍人員對廢軍品管理現況之研究─以空軍 後勤單位為例 指導教授:胡子陵博士 研究生:劉俊德.
庄文忠 副教授 世新大学行政管理学系 相关分析与简单回归分析 庄文忠 副教授 世新大学行政管理学系 SPSS之应用(庄文忠副教授) 2019/4/7.
MyLibrary ——数字图书馆的个性化服务
方差分析 方差分析的概念 单因素方差分析 有交互作用的双因素方差分析 无交互作用的双因素方差分析.
指導老師:蘇明俊 組員: 陳柔安 潘依蓮 張壹凱
第三篇 医学统计学方法. 第三篇 医学统计学方法 医学统计学方法 4 主讲人 陶育纯 医学统计学方法 4 主讲人 陶育纯 流行病与卫生统计学教研室
第八章 运动和力 第1节 牛顿第一定律和惯性 (第2课时  惯性).
多个样本均数的两两比较.
第八章 均值比较与检验 2019/5/10.
資本結構影響因素之探討 以上市公司為例 指導老師:包倩華老師 組 員:張安琪 連婉君.
2003年中華職棒大聯盟現場球迷參與動機與滿意度之研究 A Study of Participation Motivation and Satisfaction of the Chinese Professional League Baseball’ Scene Fan in 2003 真理大學 運動管理學系.
第四章 多组资料均数的比较 七年制医疗口腔《医学统计学》
第八章 方差分析(analysis of variance)
完全随机设计多样本资料的方差分析 赵耐青.
Part 3 初等統計與理論.
畢氏定理(百牛大祭)的故事 張美玲 製作 資料來源:探索數學的故事(凡異出版社).
实习二 方差分析 流行病与卫生统计学系
SAS 統計程序實作 PROC NPAR1WAY 無母數檢定
102年人事預算編列說明 邁向頂尖大學辦公室製作.
Presentation transcript:

第九章 方差分析 讲授内容 方差分析的基本思想及应用条件 完全随机设计资料的方差分析 随机区组设计资料的方差分析 重复测量资料的方差分析 多个样本均数间的多重比较

基本思想及应用条件 1、基本思想:根据资料设计的类型及研究目的,可将总变异分解为两个或多个部分,每个部分的变异可由某因素的作用来解释。通过比较可能由某因素所至的变异与随机误差,即可了解该因素对测定结果有无影响。

基本思想及应用条件 2、应用条件: 1)各样本是相互独立的随机样本 2)各样本来自正态总体 3)各处理组总体方差相等,即方差齐

完全随机设计资料的方差分析 一、完全随机设计 完全随机设计:(completely random design)是采用完全随机化的分组方法,将全部试验对象分配到g个处理组(水平组),各组分别接受不同的处理,试验结束后比较各组均数之间的差别有无统计学意义,推论处理因素的效应。

完全随机设计资料的方差分析 例1 某医生为研究一种四类降糖新药的疗效,以统一的纳入标准和排除标准选择了60名2型糖尿病患者,按完全随机设计方案将患者分为三组进行双盲临床试验。其中,降糖新药高剂量组21人、低剂量组19人、对照组20人。对照组服用公认的降糖药物,治疗4周后测得其餐后2小时血糖的下降值(mmol/L),结果如表9-1所示。问治疗4周后,餐后2小时血糖下降值的三组总体平均水平是否不同?

完全随机设计资料的方差分析 分组方法:先将60名糖尿病患者从1开始到60编号;从随机数字表(附表15)中的任一行任一列开始,依次读取三位数作为一个随机数录于编号下;然后将全部随机数从小到大编序号(数据相同的按先后顺序编序号),将每个随机数对应的序号记录;规定序号1-21为甲组,序号22-40为乙组,序号41-60为丙组。

完全随机设计资料的方差分析 完全随机设计资料的方差分析 二、变异分类与表达

完全随机设计资料的方差分析 完全随机设计资料在进行统计分析时,需根据数据的分布特征选择方法,对于正态分布且方差齐同的资料,常采用完全随机设计的单因素方差分析(one-way ANOVA)或成组资料的t检验(g=2);对于非正态分布或方差不齐的资料,可进行数据变换或采用Wilcoxon秩和检验。

完全随机设计资料的方差分析 完全随机设计资料的方差分析 记总均数为 , 各处理组均数为 , 总例数为N=nl+n2+…+ng, 记总均数为 , 各处理组均数为 , 总例数为N=nl+n2+…+ng, g为处理组数。

完全随机设计资料的方差分析 完全随机设计资料的方差分析 1. 总变异 60名2型糖尿病患者的餐后2小时血糖Xij大小各不相同,与它们的总均数(overall mean)也不相同,这种变异称为总变异(total variation)。该变异既包含了随机误差(即2型糖尿病患者的个体差异和测量误差),又包含了三组用药即处理的不同,其大小用所有数据(N=60)的方差即均方MS(mean square)来表示。

完全随机设计资料的方差分析 完全随机设计资料的方差分析 1.总变异 全部测量值大小不同,这种变异称为总变异。 1.总变异 全部测量值大小不同,这种变异称为总变异。 总变异的大小可以用离均差平方和(sum of squares of deviations from mean,SS)表示,即各测量值Xij与总均数差值的平方和,记为SS总。 总变异SS总反映了所有测量值之间总的变异程度

完全随机设计资料的方差分析 完全随机设计资料的方差分析

完全随机设计资料的方差分析 完全随机设计资料的方差分析 2.组内变异 在同一处理组中,虽然每个受试对象接受的处理相同,但测量值仍各不相同,各组内2型糖尿病患者的餐后2小时血糖Xij大小各不相同,与本组的样本均数也不相同,这种变异称为组内变异(误差)。组内变异可用组内各测量值Xij与其所在组的均数的差值的平方和表示,记为SS组内, 表示随机误差(含个体差异和测量误差)的影响。又称误差变异

完全随机设计资料的方差分析 完全随机设计资料的方差分析

完全随机设计资料的方差分析 完全随机设计资料的方差分析 3.组间变异 各处理组由于接受处理的水平不同,各组的样本均数 (i=1,2,…,g)也大小不等,三组2型糖尿病患者餐后2小时血糖的样本均数各不相同,它与总均数也不相同,这种变异称为组间变异。 其大小可用各组均数与总均数的离均差平方和表示,记为SS组间,它反映了三组用药不同的影响(如处理确实有作用),同时也包括了随机误差

完全随机设计资料的方差分析 完全随机设计资料的方差分析

完全随机设计资料的方差分析 完全随机设计资料的方差分析 各组均数之间相差越悬殊,它们与总均数的差值越大,SS组间就越大,反之SS组间越小。SS组间反映了各间的变异程度。 存在组间变异的原因有: ①随机误差(包括个体变异和测量误差) ②处理的不同水平可能对试验结果的影响。

完全随机设计资料的方差分析 完全随机设计资料的方差分析 基本思想:

完全随机设计资料的方差分析 完全随机设计资料的方差分析 变异程度除与离均差平方和的大小有关外,还与其自由度有关,将各部分离均差平方和除以相应自由度,其比值称为均方差,简称均方(mean square,MS)。

完全随机设计资料的方差分析 完全随机设计资料的方差分析 如果各组样本的总体均数相等(H0:…),即各处理组的样本来自相同总体,无处理因素的作用(处理效应),则组间变异同组内变异一样,只反映随机误差作用的大小。组间均方与组内均方的比值称为F统计量

完全随机设计资料的方差分析 完全随机设计资料的方差分析

完全随机设计资料的方差分析 完全随机设计资料的方差分析 F值(Fisher)接近于l,就没有理由拒绝H0;反之,F值越大,拒绝H0的理由越充分。数理统计的理论证明,当H0成立时,F统计量服从F分布。F分布有两个自由度, 分子自由度为1,分母自由度为2,记为F~ 。

完全随机设计资料的方差分析 完全随机设计资料的方差分析

完全随机设计资料的方差分析 完全随机设计资料的方差分析 通过上述变异的分解,可以看出,方差分析的基本思想就是根据试验设计的类型,将全部测量值总的离均差平方和及其自由度分解为两个或多个部分,除随机误差作用外,每个部分的变异可由某个因素的作用(或某几个因素的交互作用)加以解释,如组间变异SS组间可由处理因素的作用加以解释。通过比较不同变异来源的均方,借助F分布做出统计推断,从而推论各种研究因素对试验结果有无影响。

应用条件: 1)各样本是相互独立的随机样本 2)各样本来自正态总体 3)各处理组总体方差相等,即方差齐 完全随机设计资料的方差分析

完全随机设计资料的方差分析 完全随机设计资料的方差分析 分析步骤 (1)建立检验假设,确定检验水准 H0:三个总体均数全相等,即μ1=μ2=μ3 H1:三个总体均数不全相等,亦即至少有两个总体均数不等。 即μ1≠μ2≠μ3或μ1=μ2≠μ3或μ1=μ3≠μ2或μ2=μ3≠μ1 α=0.05

完全随机设计资料的方差分析 完全随机设计资料的方差分析 (2)计算检验统计量 可根据下表的公式和前面表9-1下半部分数据来计算。也可用统计软件包如SAS或SPSS等进行计算,直接获得表9-4的方差分析表。

完全随机设计资料的方差分析 完全随机设计资料的方差分析 本例的资料是按完全随机设计方法获得的试验结果, 可将总变异分解成组间变异和组内变异, 并列方差分析表.

完全随机设计资料的方差分析 完全随机设计资料的方差分析 按表中的公式计算各离均差平方和SS、自由度、均方MS和F值。

完全随机设计资料的方差分析 完全随机设计资料的方差分析 (3)确定P值,作出推断结论 以求F值时分子的自由度ν1=ν组间、分母的自由度ν2=ν组内查附表3的F界值表得P值。若F≥Fα(ν1, ν2),则P≤α,按α水准,拒绝H0,接受H1,有统计学意义。 本例:ν1=3−1=2,ν2=60−3=57。因附表3中ν2无57,故取最接近者ν2=60,得P<0.01。按α=0.05水准,拒绝H0,接受H1,有统计学意义。可以认为2型糖尿病患者经药物(新药和标准药物)治疗4周,其餐后2小时血糖的总体平均水平不全相同,即三个总体均数中至少有两个不同。

完全随机设计资料的方差分析 完全随机设计资料的方差分析 注意:方差分析的结果若拒绝H0,接受H1,不能说明各组总体均数两两间都有差别。如果要分析哪些两组间有差别,要进行多个均数间的多重比较(见本章第六节)。当g=2时,方差分析的结果与两样本均数比较的t 检验等价,有 。

随机区组设计资料的方差分析 随机区组设计 随机区组设计(randomized block design)又称为配伍组设计,是配对设计的扩展。具体做法是:先按影响试验结果的非处理因素(如性别、体重、年龄、职业、病情、病程等)将受试对象配成区组(block),再分别将各区组内的受试对象随机分配到各处理或对照组。

随机区组设计资料的方差分析 随机区组设计资料的方差分析 与完全随机设计相比,随机区组设计的特点是随机分配的次数要重复多次,每次随机分配都对同一个区组内的受试对象进行,且各个处理组受试对象数量相同,区组内均衡。在进行统计分析时,将区组变异离均差平方和从完全随机设计的组内离均差平和中分离出来,从而减小组内平方和(误差平方和),提高了统计检验效率。若将区组作为另一处理因素的不同水平,随机区组设计等同于无重复观察的两因素设计。

随机区组设计资料的方差分析 例9-2 为探索丹参对肢体缺血再灌注损伤的影响,将30只纯种新西兰实验用大白兔,按窝别相同、体重相近划分为10个区组。每个区组3只大白兔随机采用A、B、C三种处理方案,即在松止血带前分别给予丹参2ml/kg、丹参1ml/kg、生理盐水2ml/kg,在松止血带前及松后1小时分别测定血中白蛋白含量(g/L),算出白蛋白减少量如下表9-6所示,问A、B两方案分别与C方案的处理效果是否不同?

随机区组设计资料的方差分析 随机区组设计资料的方差分析 方法:先将小白鼠的体重从轻到重编号,体重相近的3只小白鼠配成一个区组,共10个区组。在随机数字表(附表15)中任选一行一列开始的2位数作为1个随机数;在每个区组内将随机数按大小排序;各区组中内序号为1 的接受A方案、序号为2的接受B方案、序号为3的接受C方案.

随机区组设计资料的方差分析 随机区组设计资料的方差分析 符合随机区组设计的资料在进行数据统计分析时也需根据数据的分布特征选择方法,对于正态分布且方差齐同的资料,应采用两因素方差分析(two-way ANOVA)或配对t检验(g=2)。当不满足方差分析和t检验条件时,可对数据进行变换或采用随机区组设计资料的Friedman M检验。

随机区组设计资料的方差分析 随机区组设计资料的方差分析 变异分解 为说明随机区组设计方差分析的变异分解和计算分析过程,将按随机区组设计的试验结果用符号表示整理成表。将第j (j=1,2,…,n)区组的受试对象随机分配接受处理因素第i(i=1,2,…,g)水平的处理,试验结果用Xij表示,整理成表4-7形式。

随机区组设计资料的方差分析 随机区组设计资料的方差分析

随机区组设计资料的方差分析 随机区组设计资料的方差分析 (1) 总变异SS总:反映所有观察值之间的变异,计算见公式(4-1)。

随机区组设计资料的方差分析 随机区组设计资料的方差分析 (4) 误差变异:完全由随机误差产生的变异,记为SS误差。

随机区组设计资料的方差分析 随机区组设计资料的方差分析

随机区组设计资料的方差分析 随机区组设计资料的方差分析 分析步骤 (1)建立检验假设,确定检验水准 对于处理组, H0:三个总体均数全相等,即A、B、C三种方案的效果相同 H1:三个总体均数不全相等,即A、B、C三种方案的效果不全相同 对于区组, H0:十个总体均数全相等 H1:十个总体均数不全相等 均取α=0.05

随机区组设计资料的方差分析 随机区组设计资料的方差分析 (2)计算检验统计量 变异来源 SS df MS F P 处理组 13.7018 2 6.8509 32.639 <0.01 区组 1.5577 9 0.1731 0.825 >0.05 误差 3.7790 18 0.2099 总 19.0385 29

随机区组设计资料的方差分析 随机区组设计资料的方差分析 按 水准,1=2、 2=18查附表3的F界值表,得F0.05(2,18)=3.55, F0.01(9,18)= 3.60,F=32.64>F0.01(2,18),P<0.01。按a水准,拒绝H0,认为三种方案的处理效果不全相等,还不能认为十个区组的总体均数不全相同。

随机区组设计资料的方差分析 随机区组设计资料的方差分析 注意:方差分析的结果拒绝H0,接受H1,不能说明各组总体均数间两两都有差别。如果要分析哪些两组间有差别,可进行多个均数间的多重比较(见本章第六节)。当g=2时,随机区组设计方差分析与配对设计资料的t 检验等价,有 。

随机区组设计资料的方差分析 随机区组设计资料的方差分析 随机区组设计确定区组因素应是对试验结果有影响的非处理因素。区组内各试验对象应均衡,区组之间试验对象具有较大的差异为好,这样利用区组控制非处理因素的影响,并在方差分析时将区组间的变异从组内变异中分解出来。因此,当区组间差别有统计学意义时,这种设计的误差比完全随机设计小,试验效率得以提高。

重复测量资料的方差分析 例9-4 为研究减肥新药盐酸西布曲明片和盐酸西布曲明胶囊的减肥效果是否不同,以及肥胖患者服药后不同时间的体重随时间的变化情况。采用双盲双模拟随机对照试验,将体重指数BMI27的肥胖患者40名随机等分成两组,一组给予盐酸西布曲明片+模拟盐酸西布曲明胶囊,另一组给予盐酸西布曲明胶囊+模拟盐酸西布曲明片。所有患者每天坚持服药,共服药6个月(24周),受试期间禁用任何影响体重的药物,而且受试对象行为、饮食及运动与服药前的平衡期均保持一致。分别于平衡期(0周)、服药后的8周、16周、24周测定肥胖患者的体重(kg)得表9-13的资料。

受试 对象j 剂型 k 服药后测定时间i(周) 8 16 24 1 84.4 82.2 83.0 21 2 64.4 61.4 61.8 62.0 105.0 100.8 97.4 96.6 22 91.0 88.4 87.4 89.6 3 63.8 61.6 60.4 23 76.0 76.2 72.8 71.6 4 86.2 85.5 81.8 71.0 72.0 69.8 68.4 5 75.6 73.4 74.0 73.0 25 69.4 66.6 62.8 60.8 6 61.2 60.2 26 89.9 92.6 95.5 7 67.8 66.0 63.4 63.6 27 66.8 62.6 77.2 73.6 72.6 28 9 73.2 72.2 74.6 29 70.0 67.6 10 65.4 30 86.6 84.0 81.4 78.0 11 80.0 77.0 72.4 31 90.4 77.4 12 74.4 75.2 32 74.8 76.6 13 82.6 80.4 81.2 79.6 33 67.4 61.0 58.2 14 68.6 65.0 63.2 34 80.2 75.4 15 79.0 73.8 72.5 35 76.5 78.5 36 83.2 17 68.2 70.2 37 68.7 65.8 63.0 66.4 18 38 78.4 19 39 66.5 20 76.4 40 64.6 64.2

重复测量资料的方差分析 重复测量资料和随机区组设计资料的区别: (1)重复测量资料中同一受试对象(看成区组)的数据高度相关,无论哪位受试对象服用盐酸西布曲明片剂或是胶囊,其服药后8周、16周和24周的体重均和前面时间点(含服药前的0周)的体重相关。表9-14为分不同剂型后使用统计软件包计算得到的各时点简单相关系数r,从中可以看出,不同时点间相关系数介于0.850 ~0.989之间,其P值全为0.000,均有统计学意义,说明不同时点数据其相关性较强。

(k=1时)服药后测定时间i 服药后测定时间i(周) (k=2时)服药后测定时间i 8 16 24 0.989 0.971 0.939 0.944 0.850 0.986 0.966 0.961 0.880 0.985 0.958

重复测量资料的方差分析 (2)重复测量资料中的处理因素在受试对象(看成区组)间为随机分配,但受试对象(看成区组)内的各时间点往往是固定的,不能随机分配;随机区组设计资料中每个区组内的受试对象彼此独立,处理只在区组内随机分配,同一区组内的受试对象接受的处理各不相同。 本节主要介绍两因素重复测量资料的单变量方差分析方法。

重复测量资料的方差分析 离均差平方和与自由度的分解 两因素重复测量资料的总变异包括两部分: 横向分组的受试对象间(between subjects)的变异 纵向分组的受试对象内(within subjects)的变异。 其中横向分组受试对象间的变异又分为处理因素K(在此为剂型)的变异和个体间误差的变异两部分;而纵向分组受试对象内的变异则可分为时间因素I的变异、处理K和时间I的交互作用(KI)以及个体内误差的变异三部分

重复测量资料的方差分析 重复测量资料方差分析的基本步骤 重复测量资料的方差分析步骤仍为三步,本例如下: (1)建立检验假设,确定检验水准 对于处理因素K H0:不同剂型(片剂和胶囊)的减肥效果相同 H1:不同剂型(片剂和胶囊)的减肥效果不同

重复测量资料的方差分析 对于时间因素I H0:服用减肥药前后不同时间体重的总体均数全相等 H1:服用减肥药前后不同时间体重的总体均数不全相等 对于交互作用KI H0:药物剂型K和时间I无交互效应 H1:药物剂型K和时间I有交互效应 均取α=0.05

重复测量资料的方差分析 (2)计算检验统计量 使用统计软件包SAS或SPSS等进行计算。对本例可得到表9-15的方差分析表 以求F值时分子自由度ν1、分母自由度ν2查附表3的F界值表得相应P值,或直接由计算机所给P值作出推断结论。本例,按α=0.05水准,减肥药剂型K(片剂和胶囊),剂型K与时间I的交互效应KI均不拒绝H0,无统计学意义,还不能认为盐酸西布曲明不同剂型的减肥效果不同,也还不能认为剂型K与时间I间有交互效应。而时间因素I拒绝H0,接受H1,有统计学意义,可认为服用减肥药盐酸西布曲明前后不同时间(8周、16周和24周)的平均体重不全同。

变异来源 SS df MS F P (受试对象间) (13163.9810) (39) 处理K 5.9290 1 0.017 0.897 个体间误差 13158.0520 38 346.2645 (受试对象内) (904.6500) (120) 时间I 384.5300 3 128.1767 28.213 0.000 交互作用KI 2.1940 0.7313 0.161 0.922 个体内误差 517.9260 114 4.5432 总 14068.6310 159 88.4820

重复测量资料的方差分析 重复测量资料方差分析的前提条件 进行重复测量资料的方差分析,除需满足一般方差分析的条件外(详后),还需特别满足协方差阵(covariance matrix)的球形性(sphericity / circularity)或复合对称性(compound symmetry)。Box(1954)指出,若球形对称性质不能满足,则方差分析的F值是有偏的,因为它增大了第一类错误的概率。球对称性通常采用Mauchly检验(Mauchly’s test)来判断

 多个样本均数间的多重比较 当方差分析的结果为拒绝H0,接受H1时,只说明g个总体均数不全相等。若想进一步了解哪些两个总体均数不等,需进行多个样本均数间的两两比较或称多重比较(multiple comparison)。 若用上一章的两样本均数比较的t检验进行多重比较,将会加大犯Ⅰ类错误(把本无差别的两个总体均数判为有差别)的概率。

 多个样本均数间的多重比较 第六节           多个样本均数间的多重比较

多个样本均数间的多重比较 第六节 多个样本均数间的多重比较  多个样本均数间的多重比较 第六节           多个样本均数间的多重比较 下面介绍三种多重比较方法:Bonfferoni t检验 、Dunnett-t检验和SNK-q检验。 Bonfferoni t检验 (1)建立检验假设,确定检验水准 H0:μA=μB,即任两对比组的总体均数相等 H1:μA≠μB,即任两对比组的总体均数不等

 多个样本均数间的多重比较 第六节           多个样本均数间的多重比较 (2)计算检验统计量 (3)确定P值,作出推断结论

多个样本均数间的多重比较 第六节 多个样本均数间的多重比较 对比组 A与B (1) (2) (3) t P (4) (5) 1与3  多个样本均数间的多重比较 第六节         多个样本均数间的多重比较 对比组 A与B (1) (2) (3) t P (4) (5) 1与3 3.7652 1.2483 3.016 0.002~0.005 1与2 3.3952 1.2650 2.684 0.005~0.01 2与3 0.3700 1.2800 0.289 >0.50 两均数之差 标准误

 多个样本均数间的多重比较 第六节         多个样本均数间的多重比较 Dunnett- t 检验 适用于g-1个实验组与一个对照组均数差别的多重比较 ,检验统计量为t, 亦称t检验。 计算公式为: ν=ν误差

 多个样本均数间的多重比较 第六节         多个样本均数间的多重比较

H0:μT=μC,即任一实验组与对照组的总体均数相等 H1:μT≠μC,即任一实验组的总体均数低于对照组的总体均数 α=0.05  多个样本均数间的多重比较 第六节         多个样本均数间的多重比较 对例9-2资料,本例也可用LSD-t检验,但一般用Dunnett-t 检验。 (1)建立检验假设,确定检验水准 H0:μT=μC,即任一实验组与对照组的总体均数相等 H1:μT≠μC,即任一实验组的总体均数低于对照组的总体均数 α=0.05

多个样本均数间的多重比较 第六节 多个样本均数间的多重比较 (2)计算检验统计量 对比组 T与C (1) (2) tD P A与C  多个样本均数间的多重比较 第六节         多个样本均数间的多重比较 (2)计算检验统计量 对比组 T与C (1) (2) tD P A与C -1.5900 -7.760 <0.01   B与C -1.1940 -5.827 两均数之差

多个样本均数间的多重比较 第六节 多个样本均数间的多重比较 (3)确定P值,作出推断结论  多个样本均数间的多重比较 第六节         多个样本均数间的多重比较 (3)确定P值,作出推断结论 将表中tD取绝对值,并以计算MS误差时的自由度ν误差=18和实验组数a=k−1=2(不含对照组)查附表5的Dunnett t界值表得P值,按α=0.05水准,A方案与C方案、B方案与C方案均拒绝H0,接受H1,有统计学意义。可以认为A方案与C方案、B方案与C方案大白兔血中白蛋白的减少量不同。

 多个样本均数间的多重比较 第六节         多个样本均数间的多重比较 SNK-q检验 SNK(Student-Newman-Keuls)检验,亦称q检验,适用于多个样本均数两两之间的全面比较。检验统计量q的计算公式为: ν=ν误差

 多个样本均数间的多重比较 第六节         多个样本均数间的多重比较

 多个样本均数间的多重比较 第六节         多个样本均数间的多重比较 例9-5 请对例9-1资料中治疗4周后,餐后2小时血糖下降值的三组总体均数进行两两比较。 (1)建立检验假设,确定检验水准 H0:μA=μB,即任两对比较组的总体均数相等 H1:μA≠μB,即任两对比较组的总体均数不相等

多个样本均数间的多重比较 第六节 多个样本均数间的多重比较 (2)计算检验统计量 首先将三个样本均数由小到大排列,并编组次  多个样本均数间的多重比较 第六节         多个样本均数间的多重比较 (2)计算检验统计量 首先将三个样本均数由小到大排列,并编组次 (3)确定P值,作出推断结论 组别 高剂量组(i=1) 低剂量组(i=2) 对照组(i=3) 9.1952 5.8000 5.4300 组次 1 2 3

多个样本均数间的多重比较 第六节 多个样本均数间的多重比较 对比组 A与B (1) (2) , (3) 对比组内 包含组数a q临界值 P  多个样本均数间的多重比较 第六节         多个样本均数间的多重比较 两均数之 差标准误 对比组 A与B (1) (2) , (3) 对比组内 包含组数a q临界值 P 0.05 0.01 (5) (6) (7) (8) 1与3 3.7652 0.8827 4.266 3 3.40 4.28 0.01~0.05 1与2 3.3952 0.8945 3.796 2 2.83 3.76 <0.01 2与3 0.3700 0.9051 0.409 >0.05 两均数之差 q