1、12一、析因设计资料的方差分析一、析因设计资料的方差分析(Factorial experimental design) 是一种将是一种将两个或多个因素两个或多个因素的各水平交叉分组的各水平交叉分组进行实验的设计。它不仅可检验各因素内部不进行实验的设计。它不仅可检验各因素内部不同水平间有无差异,还可检验同水平间有无差异,还可检验2个或多个因素个或多个因素间是否存在间是否存在交互作用交互作用。多因素实验。多因素实验 3析因设计析因设计 各因素各水平的全面组合各因素各水平的全面组合处理组合数处理组合数 g = 各因素水平数之积。各因素水平数之积。常见的设计模型常见的设计模型: 22 222 32 2
2、23 注意:注意:析因设计时,分析析因设计时,分析的因素数和水平数不宜过的因素数和水平数不宜过多,一般因素数不超过多,一般因素数不超过4,水平数不超过水平数不超过3。4实验对象的分组:实验对象的分组:完全随机设计:各组随机分配完全随机设计:各组随机分配 n 个试验个试验 对象,总对象数为对象,总对象数为 gn。随机区组设计:随机区组设计: n 个区组,每个区组个区组,每个区组 g 个个 试验对象随机分配。试验对象随机分配。5 例例11-1 将将20只家兔随机等分只家兔随机等分4组,每组组,每组5只,只,进行神经损伤后的缝合试验。处理由进行神经损伤后的缝合试验。处理由A、B两因素组两因素组合而成
3、,合而成,因素因素A为缝合方法,有两水平为缝合方法,有两水平:外膜缝合外膜缝合a1、束膜缝合束膜缝合a2;因素因素B为缝合后时间,亦有两水平:缝为缝合后时间,亦有两水平:缝合后合后1月月b1、缝合后、缝合后2月月b2。试验结果为家兔神经缝。试验结果为家兔神经缝合后的轴突通过率合后的轴突通过率(%)(注:(注:测量指标,视为计量资测量指标,视为计量资料料),见表),见表11-1。欲用析因分析比较不同缝合方法。欲用析因分析比较不同缝合方法及缝合后时间对轴突通过率的影响及缝合后时间对轴突通过率的影响。6 表11-1 家兔神经缝合后的轴突通过率(%)2740 /2027380, 34800273807
4、420CSS总A(缝合方法) 外膜缝合( a1) 束膜缝合 (a2) B(缝合后时间) 1 月( b1) 2 月( b2) 1 月(b1) 2 月(b2) 合 计 10 30 10 50 10 30 20 50 40 70 30 70 50 60 50 60 10 30 30 30 iX 24 44 28 52 iT 120 220 140 260 740(X) 2iX 4400 11200 4800 4400 34800(X2) 7 B 因素 ( 2 水平 ) A 因素 (2 水平) 缝合后1 月 缝合后2 月 (b1) (b2) 外膜缝合(a1) 24 (a1b1) 44 (a1b2) 束
5、膜缝合(a2) 28 (a2b1) 52 (a2b2) 图11-1 2因素2水平析因试验示意图 将表11-1的4组数据的均数整理成图11-1,现分析A因素不同水平、B因素不同水平的单独效应、主效应和交互作用。8 B因素 A因素 b1 b2 平 均 b2b1 a1 24 44 34 20 a2 28 52 40 24 26 48 22 平均 a2a1 4 8 6 表11-2 2因素2水平析因试验的均数差别 B的主效应的主效应B的单独效应的单独效应A的单独效应的单独效应A的主效应的主效应9B因素B1B2平均 B2-B1A因素A1A1B1A1B2S1A2A2B1A2B2S2平均SA2-A1T1T2T
6、B的主效应的主效应B的单独效应的单独效应A的单独效应的单独效应A的主效应的主效应10 1. 1. 单独效应单独效应(simple effect) 指其他因素的水指其他因素的水平固定时,同一因素不同水平间的平固定时,同一因素不同水平间的差别差别 2. 2. 主效应主效应 (main effect)指某一因素各水平间的指某一因素各水平间的平平均差别均差别120 B= 24AB=2时, 因素的单独效应 因素固定在时, 因素的单独应水平效水平14 A=B8 A=2时, 因素的单独效应 因素固定在时, 因素的单效应水独平水平A的主效应的主效应: (4+8)/ 2=611221 22 11 1()()2(
7、84) 22ABa baba bab2 22 11 21 1()()2(2420) 22BAa ba babab即即ABBA。 3. 3. 交互作用交互作用(interaction) 当某因素的各个单独效应随另一因素变化而变化时,则称这两个因素间存在交互作用。12缝合2月 (b2)缝合1月 (b1) 4个均数可作线图,若两条直线几乎相互平行, 则表示两因素交互作用很小;若两条直线相互不平行, 则说明两因素可能存在交互作用。134方差分析方差分析 表11-2中,A因素(缝合方法)的主效应为6%,B因素(缝合时间)的主效应为22%,AB的交互作用表示为2%。以上都是样本均数的比较结果,要推论总体均
8、数是否有同样的特征,需要对试验结果作假设检验即方差分析后下结论。14变异来源 自由度 SS MS F P 总变异 19 7420 (处理组间) (3) (2620) A 1 180 180 0.60 0.05 B 1 2420 2420 8.07 0.05 AB 1 20 20 0.07 0.05 误 差 16 4800 300 表11-5 表11-1析因试验结果方差分析表 析因设计的析因设计的方差分析方差分析: 总变异的离均差平方和总变异的离均差平方和SS及自由度及自由度分解为各因分解为各因素的作用、两因素的交互作用及误差等部分。素的作用、两因素的交互作用及误差等部分。15 结合样本均数的比
9、较结果,A因素的主效应为6%,AB的交互作用为2%,均不具有统计学意义,仅B因素(缝合后时间)的主效应22%有统计学意义。 结论:尚不能认为两种缝合方法对神经轴突通过率有影响;可以认为缝合后2月与1月相比,神经轴突通过率提高了。16析因设计析因设计优点优点:是一高效的实验设计方法,:是一高效的实验设计方法,不仅能分析各因素内部不同水平间有无差别,不仅能分析各因素内部不同水平间有无差别,还可分析各因素间的还可分析各因素间的交互作用交互作用。缺点缺点:与正交实验设计相比,属于全面试验,:与正交实验设计相比,属于全面试验,因此,研究的因素个数和因素的水平数不宜因此,研究的因素个数和因素的水平数不宜过
10、多。过多。 17原则原则:若存在交互作用,需逐一分析各因素:若存在交互作用,需逐一分析各因素的单独效应,此时分析主效应无意义;的单独效应,此时分析主效应无意义;若不存在交互作用,则两因素的作用相互独若不存在交互作用,则两因素的作用相互独立,只需考虑各因素的主效应即可。立,只需考虑各因素的主效应即可。18正交实验设计正交实验设计:利用一套规范化的:利用一套规范化的正交表正交表,使,使每次试验的各因素及其水平得到合理安排的高每次试验的各因素及其水平得到合理安排的高效多因素实验设计。效多因素实验设计。 多因素实验多因素实验 常用于常用于:寻找疗效好的药物配方、医疗仪器多:寻找疗效好的药物配方、医疗仪
11、器多个参数的优化组合、生物体的培养条件等最优个参数的优化组合、生物体的培养条件等最优搭配方案的研究。搭配方案的研究。二、正交设计资料的方差分析二、正交设计资料的方差分析(Orthogonal experimental design) 19与析因实验相比,与析因实验相比,正交设计:正交设计:是析因试验的部分实施是析因试验的部分实施(可减少多因素实验可减少多因素实验的次数的次数)。可成倍的减少实验次数可成倍的减少实验次数(以牺牲分析各因素(以牺牲分析各因素的交互作用为代价)。的交互作用为代价)。 g=22222=32 g=16 或或 g=8只分析只分析有意义的主效应有意义的主效应和和部分部分重要因
12、素的重要因素的一一级交互作用。级交互作用。 20 列列 号号处理号处理号12345671111111121112222312211224122221152121212621221217221122182212112表表11-13 L11-13 L8 8(2 27 7)正交设计表)正交设计表 21(一一) 正交表的几个基本概念:正交表的几个基本概念:1. 正交表的符号正交表的符号:Ln(Km) L:表示正交表:表示正交表(orthogonal layout) n:表示正交表有:表示正交表有n行,代表实验次数行,代表实验次数 m:表示正交表有:表示正交表有m列,表示列,表示最多最多允许安允许安排的
13、因素排的因素(及其交互作用及其交互作用)的个数的个数 K:表示各因素的水平数,即每列中的数:表示各因素的水平数,即每列中的数字为字为1、2k22L8(27):表示最多可安排表示最多可安排7 7个个2 2水平的因素要水平的因素要做做8 8次实验的正交表次实验的正交表L L1616(2(21515) )L L3232(2(23131) )L L9 9(3(34 4) )L L1818(3(37 7) )23 列列 号号实验号实验号12345671111111121112222312211224122221152121212621221217221122182212112表表11-13 L11-13
14、 L8 8(2 27 7)正交设计表)正交设计表 242.正交表的性质正交表的性质: 每列中不同数字出现的次数相同;每列中不同数字出现的次数相同; 任两列同一横行的任两列同一横行的2个有序数对出现的次数个有序数对出现的次数相等。相等。3. 正交表的交互作用表正交表的交互作用表每个正交表均有其对应的交互作用表。每个正交表均有其对应的交互作用表。25 列列 号号列号列号1 12 23 34 45 56 67 71 13 32 25 54 47 76 62 21 16 67 74 45 53 37 76 65 54 44 41 12 23 35 53 32 26 61 1表表11-14a L11-1
15、4a L8 8(2(27 7) )交互作用表交互作用表 264.正交表分类正交表分类: 相同水平正交表相同水平正交表: 即各列的水平数相同即各列的水平数相同 混合水平正交表混合水平正交表: 即各列的水平数不即各列的水平数不(全全)等等27(二二)正交表的选用正交表的选用1.根据研究目的和专业知识,确定实验的因根据研究目的和专业知识,确定实验的因素个数,并明确主要因素。素个数,并明确主要因素。2.根据各因素的水平数,确定选用哪类正交根据各因素的水平数,确定选用哪类正交表(相同水平或混合水平)。表(相同水平或混合水平)。3.再根据因素个数、可能存在的交互作用,再根据因素个数、可能存在的交互作用,确
16、定选多少列即多大的正交表确定选多少列即多大的正交表 。28(三三)表头设计表头设计 利用所选用的利用所选用的正交表正交表及其相应的及其相应的交互作用交互作用表表,将各因素及要分析的可能存在的交互,将各因素及要分析的可能存在的交互作用安排进所选正交表各列的过程。作用安排进所选正交表各列的过程。29 例例11-4 P247研究雌螺产卵的最优条件,在研究雌螺产卵的最优条件,在20cm2的泥盒里饲养同龄雌螺的泥盒里饲养同龄雌螺10只,试验条件有只,试验条件有4个因素个因素(表(表11-15),每个因素),每个因素2个水平。试在考虑温度与个水平。试在考虑温度与含氧量对雌螺产卵有交互作用的情况下安排正交试
17、含氧量对雌螺产卵有交互作用的情况下安排正交试验。验。 表表11-15 雌螺产卵条件因素与水平雌螺产卵条件因素与水平因素因素水平水平A因素因素温度温度()B因素因素含氧量含氧量(%)C因素因素含水量含水量(%)D因素因素pH值值150.5106.02255.0308.030正交设计的步骤: 1.正交表的选择: L8(27) 2.表头设计: 3.确定各组的试验条件:31注:此表用于安排各因素于各列注:此表用于安排各因素于各列 列列 号号列号列号1 1 A A2 2 B B3 34 4 C C5 56 67 7 D D1 13 32 25 54 47 76 62 21 16 67 74 45 53
18、37 76 65 54 44 41 12 23 35 53 32 26 61 1表表 L L8 8(2(27 7) )交互作用表交互作用表 A AB B表头设计:表头设计:32 表表11-14 L11-14 L8 8(2 27 7)正交设计表的表头设计)正交设计表的表头设计 因素因素 实施实施 列列 号号 个数个数 比例比例1 12 23 34 45 56 67 74 41/21/2A AB BABABC CD D 表头设计:表头设计:33 列列 号号试验试验序号序号1 A2 B3 4 C567 D111111112111222231221122412222115212121262122121
19、7221122182212112表表11-13 L11-13 L8 8(2 27 7)正交设计表)正交设计表 注:此表用于安排实验注:此表用于安排实验确定试验条件:确定试验条件:34正交设计的统计分析正交设计的统计分析: 直接分析:直接分析:P250 方差分析方差分析:总变异的离均差平方和:总变异的离均差平方和SS及自及自由度由度分解为各因素不同水平间、两因素的交分解为各因素不同水平间、两因素的交互作用及误差等部分。互作用及误差等部分。35处理组:非各因素各水平的全面组合,而是各因处理组:非各因素各水平的全面组合,而是各因素按隶属关系系统分组,各因素水平没有交叉。素按隶属关系系统分组,各因素水
20、平没有交叉。 a1 a2 a3 b1 b 2 b 3 b4 b 5 b 6 b7 b 8 b 9 三、嵌套试验三、嵌套试验(nested design)36析因设计析因设计:g 个处理全部都作用于个处理全部都作用于同一级别同一级别的实验单位的实验单位。裂区设计裂区设计:A 因素的因素的 I 个水平作用于一级实个水平作用于一级实验单位,验单位, B 因素的因素的 J 个水平作用于个水平作用于二级实验二级实验单位单位。两因素析因设计的特殊形式。两因素析因设计的特殊形式。五、五、 裂区设计裂区设计(split-plot design)37在相同试验条件下,通过改进实验设计方法可以提高实验效率。 注意多因素试验与多向分类方差分析的区别,如随机区组试验和两因素析因试验,前者是单因素试验,后者是两因素试验,但数据分析都是采用双向分类方差分析。