《方差分析二析因正交重复测量.ppt》由会员分享，可在线阅读，更多相关《方差分析二析因正交重复测量.ppt（63页珍藏版）》请在三一办公上搜索。
1、1,方差分析（二）analysis of variance ANOVA,温州医学院环境与公共卫生学院 叶晓蕾,2,前面内容回顾,1、完全随机设计的ANOVA2、随机区组设计的ANOVA3、交叉设计的ANOVA,所关心的问题：一个处理因素不同处理水平间的均数有无差异？,以上第2个设计中，设立单位组（区组）的目的是控制混杂因素。使混杂因素在各处理水平间达到均衡，提高检验效率。,第3个设计呢？,3,六.析因设计（factorial design）的ANOVA,是将两个或多个因素的各个水平进行排列组合，交叉分组进行试验，检验各因素各水平之间的差异有无统计学意义，而且可以检验因素间的交互作用。完全交叉分
2、组试验设计,优点：全面均衡、高效。缺点：工作量大。实验组合数=每个因素水平数的乘积。如四因素三水平的实验组合数=34=81。多因素多水平的实验可采用正交设计(非全面实验)。,4,例 A、B两药治疗12名贫血病人，性别、年龄一致，随机分成四组，治疗后一个月测得血中红细胞增加数结果如下表，A、B两药的治疗效果如何？两药是否存在交互效应？,5,析因设计资料的特点,2个或以上（处理）因素（factor）(分类变量)，每个因素有2个或以上水平（level）；每一组合涉及全部因素，每一因素只有一个水平参与；每一组合中至少有 2个或以上的观察值；观测值为定量数据（需满足随机、独立、正态、等方差的ANOVA条
3、件）。,6,7,一、概念,因子是指影响试验结果或观察指标的因素；因素的不同状态为水平。,1.因子（因素factor）和水平(levels)：,8,2.单独效应（simple effect)：,指其它因素固定在一个水平时，余下的一个因素不同水平之间均数的差别。,9,2.单独效应（simple effect)：,指其它因素固定在一个水平时，余下的一个因素不同水平之间均数的差别。例：用B药水平时，A药的单独效应1.10 不用B药水平时，A药的单独效应0.40,10,3.主效应（main effect)：,指在某一因素各水平间的平均差别。与单独效应不同的是，它综合了其它因素各水平与该因素每一水平所有组
4、合的情况。每个因子对观察指标的作用。有几个因子就应该有几个主效应。,11,3.主效应（main effect)：,指在某一因素各水平间的平均差别。与单独效应不同的是，它综合了其它因素各水平与该因素每一水平所有组合的情况。每个因子对观察指标的作用。有几个因子就应该有几个主效应。例：A因素的主效应：0.75。即A因素单独效应的平均值。,12,4.交互效应（interaction)：因素之间联合对观察指标产生的影响。当某因素的各单独效应随另一因素变化而变化时，则称这两个因素间存在交互效应。,13,4.交互效应（interaction)：因素之间联合对观察指标产生的影响。当某因素的各单独效应随另一因素
5、变化而变化时，则称这两个因素间存在交互效应。在B药存在时，A药增加Hb作用为：2.1-1.0=1.1单纯用A药增加Hb作用为：1.2-0.8=0.4B药有加强A药升Hb作用为：1.1-0.4=0.7,14,一级交互作用：两个因子间（2-way)的交互作用。,二级交互作用：三个因子间(3-way)的交互作用。,析因分析中，应先重点考察各因素间的交互作用，若不存在交互作用，则分析某一因素的作用只需要考察该因素的主效应；若存在交互作用，则各因素的主效应检验结果已无实际意义，应当按各因素各种水平的组合来研究，即逐一分析各因素的单独效应。,15,6.随机效应型：该因素中各水平是由总体中随机抽样所得的。如
6、要研究什么温度下催化剂的效果最好，随机抽取了30、40、50三个水平。,5.固定效应型：该因素中各水平是固定的。如性别因素分为男、女两水平。,16,二、用途：,作多因素分析。1、分析每个因素的单独效应和主效应；2、考察因素之间的交互作用。,22析因设计目的：分析A、B两个因素对结果有无影响（主效应），A与B有无交互作用。,17,两因素析因设计：,三、析因设计资料方差分析中变异的分解,18,四、两因素析因设计方差分析变异的计算,19,1.建立检验假设，确定检验水准：（1）对A药：H0：用与不用血中红细胞增加数的总体均数相等。H1：用与不用血中红细胞增加数的总体均数不相等。（2）对B药：H0：用与
7、不用血中红细胞增加数的总体均数相等。H1：用与不用血中红细胞增加数的总体均数不相等。（3）对交互作用 H0：A药与B药无交互作用。H1：A药与B药有交互作用。2.求检验统计量：3.确定P，作出统计推断结论,五、析因设计资料方差分析的基本步骤,20,六、两因素析因设计资料方差分析F值的计算,21,SPSS操作演示：（1）“A、B两药治疗12名贫血病人”（2）P.194例11-1（3）P.198例11-2,22,23,24,七、正交试验设计的方差分析(orthogonal design),正交试验设计是利用一套规格化的正交表，将个试验因素、各水平之间的组合进行均匀搭配，合理安排，是一种高效的、多因
8、素试验设计方法。适用情况：当实验涉及的因素在三个或三个以上，而且因素间可能存在交互作用时。,25,（一）正交试验设计基本方法 1、正交表是正交设计的核心，是试验设计中合理安排试验并进行统计分析的主要工具。常用的正交表：L4(23)，L8(27)，L9(34)，L8(4124)，等。,26,混合正交表：L16(41，212)正交表：最多可安排13个因素，其中1个因素为4水平，12 个因素为2水平，要作16次试验。,27,正交表的特点：,1）每一列中，不同数字（水平）出现的次数相同；2）任意两列中，将同一横行的两个数看成有序的数对，每种数对出现的次数相等。所以，正交表充分体现了：分布均衡，整齐可比
9、，具有充分的代表性。,28,2、交互作用表用于安排有交互作用的试验。对应于每一张正交表都有一张交互作用表。,29,30,31,若A、B两因素安排于1、2列，如果A、B存在交互作用，就表现在第3列上，此时，第3列上就不能安排第3因素了。若在第3列安排了因素C，则第3列既是C的主效应，又是A、B的交互作用，导致效应的混杂。若A、B不存在交互作用，第3列可安排因素C，第3列也可空着，作为误差的计算来源。,32,3.正交试验的基本步骤：1、明确试验目的，确定试验指标。最好是定量指标。2、确定试验的因素和水平：3、选正交表安排实验计划：即作表头设计。根据研究的具体问题，实验次数，水平多少选用合适的正交表
10、。4、根据安排的计划进行实验，测定各个实验结果。5、对实验结果进行计算分析。,33,例11-4 P.202 研究雌螺产卵的最优条件，在20cm的泥盒里饲养同龄雌螺10只，试验条件有4个因素，每个因素2个水平。试在考虑温度与含氧量对雌螺产卵有交互作用的情况下安排正交试验。,本题：试验目的提高雌螺产卵数量；试验指标产卵数量;因素A，B，C，D；AB 水平2个,雌螺产卵条件因素与水平,34,因素影响试验指标的因素。已知对指标影响不大的因素，可不予考虑；对指标影响大小已经了解的因素，可固定在相应水平上。主要应考察的因素：重要的但影响尚不清楚的因素，以及因素间可能存在的不容忽视的交互作用（一般只考察一阶
11、交互作用）。水平各因素可能处的状态。最好水平数相同，一般取2水平。重要的因素可多取些水平。,35,表头设计 选择适当的正交表，将因素放入正交表表头的工作。,首先根据试验因素的水平数选相同水平个数的正交表，然后使要考察的因素个数加上要考察的交互作用个数少于（至多等于）正交表的列数，并且注意要避免使因素和交互作用或不同交互作用在同一列上重叠（“混杂”），而不考察的交互作用可以混杂。本例：水平数=2；因素4+交互作用1=5 故选取L8（27）。注意：若采用方差分析，则需要留有空列，或作重复测定，以便计算误差。,36,操作口诀：“因素顺序上列，水平对号入座，实验横着做”。,37,正交试验结果的分析：,
12、有直观法和方差分析法。方差分析法：,38,H0:各因素的作用及有关交互作用都不存在H1:各因素的作用及有关交互作用存在=0.05,SS总=SSA+SSB+SSC+SSD+SSAB+SS误差FA=MSA/MS误差FB=MSB/MS误差FAB=MSAB/MS误差,39,SPSS,40,故：取C2A2B1。,41,42,1、作表头设计：,因素 A B C AC H D AD BD E CD F G,(交互)AB AC BC ABC AD BD ABD CD ACD BCD ABCD,列号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15,选L16(215)正交表,因素及效应：A
13、、B、C、D、E、F、G、AC、AD、BD、CD,43,因素及效应：A、B、C、D、E、F、G、AC、AD、BD、CD,44,45,H0:各因素的作用及有关交互作用都不存在H1:各因素的作用及有关交互作用存在=0.05,SS总=SSA+SSB+SSC+SSD+SSE+SSF+SSG+SSH+SSAC+SSAD+SSBD+SSCD+SS误差FA=MSA/MS误差FB=MSB/MS误差FCD=MSCD/MS误差,SPSS,46,47,48,49,胆汁滴数增量越大越好，可见A2C2、B2D2、C2D2组合较好。综上所述，茵陈胆道汤A：金钱草 B：大黄 C：木香 D：黄芩 E：茵 陈 G：栀子 有利胆
14、作用。,50,八、重复测量资料的方差分析(repeated measurement data),重复测量资料是同一受试对象的同一观察指标在不同的时点上进行多次测量所得的资料，常用来分析该观察指标在不同时间上的变化特点。,例12-3 将手术要求基本相同的15名患者随机分3组，在手术过程中分别采用A，B，C三种麻醉诱导方法，在T0（诱导前）、T1、T2、T3、T4五个时相测量患者的收缩压，数据记录见表12-17。试进行方差分析。,51,不同麻醉诱导时相患者的收缩压（mmHg）,与随机区组设计资料的区别：（1）重复测量资料中同一受试对象（看成区组）的数据高度相关。（2）重复测量资料中的处理因素在受试
15、对象间为随机分配，但受试对象内的各时点往往是固定的。,52,两因素重复测量资料的单变量方差分析：,变异分解：,SS总=SS受试对象间 SS受试对象内=(SS处理SS个体间误差)(SS时间SS处理与时间交互SS个体内误差)总=受试对象间 受试对象内=(处理 个体间)(时间处理与时间交互个体内),53,重复测量资料方差分析的前提条件:,（1）一般方差分析的条件；（2）协方差阵（covariance matrix）的球形性(球对称)。,球对称性常用Mauchly检验(Mauchlys test)来判断。若资料不满足球对称性，即按=0.10水准拒绝H0，接受H1，则应对受试对象内所有变异的自由度（包括
16、时间效应、处理与时间交互效应以及个体内误差）均进行校正。否则会增大第一类错误的概率。,54,重复测量资料方差分析的基本步骤,建立检验假设，确定检验水准：对于处理因素：H0:不同诱导方法的收缩压的总体均数相同 H1:不同诱导方法的收缩压总体均数不全相同对于时间因素：H0:不同麻醉诱导时相收缩压的总体均数相同 H1:不同麻醉诱导时相收缩压的总体均数不全相同对于交互作用：H0:诱导方法与诱导时相无交互效应 H1:诱导方法与诱导时相有交互效应均取=0.05,55,2.计算检验统计量：利用统计软件SPSS。,受试对象间：,SPSS,56,受试对象内：,57,球对称检验：,58,各次重复测量间变化趋势分析：,59,60,结论：不同诱导方法的收缩压总体均数存在差别；不同麻醉诱导时相收缩压的总体均数不全相同；诱导方法与诱导时相有交互效应。,61,62,这是什么设计类型的资料?,63,是重复测量资料吗？,
