《《两因素方差分析》ppt课件》由会员分享，可在线阅读，更多相关《《两因素方差分析》ppt课件（92页珍藏版）》请在金锄头文库上搜索。
1、第三节 两因素试验的 方差分析,考查两个因素对试验指标的影响情况,3.1 交叉分组资料（crossover classification）的方差分析 设试验考察A、B两个因素，A因素分a个水平，B因素分b个水平 。 所谓交叉分组是指A因素每个水平与B因素的每个水平都要搭配 ，两者交叉搭配形成ab个水平组合即处理，试验因素A 、B在试验中处于平等地位 。如果将试验单元分成 ab 个组，每组随机接受一种处理 ，因而试验数据也按两因素两方向分组，这种试验数据资料称为两向分组资料，也叫交叉分组资料。 分无重复观测值和重复观测值两种类型。,对于A、B两个试验因素的全部ab个水平组合，每个水平组合只有一个
2、观测值（无重复）， 全试验共有ab个观测值，其数据模式如下表所示。,3.1.1 两因素无重复试验资料的方差分析,表 两因素无重复观测值的试验数据模式,注：A因素有a个水平，B因素有b个水平，共计有ab个水平组合，每一组合观测一次，有ab个观测值（表5），xij 为A的第i水平与B的第j水平组合观测值。,两因素无重复观测值试验资料的数学模型为： 式中， 为总平均数；,(5-26),i，j分别为Ai、Bj的效应； i=i-， j=j-， i、j分别为Ai、Bj观测值总体平均数， 且i=0，j=0； ij为随机误差，相互独立，且服从N(0，2),A因素的每个水平有b次重复，B因素的每个水平有a次重复
3、，每个观测值同时受到A、B 两因素及随机误差的作用。因此全部 ab 个观测值的总变异可以分解为 A 因素水平间变异、B因素水平间变异及试验误差三部分；自由度也相应分解。,离差平方和与自由度的分解如下：,矫正数 总平方和 A因素离差平方和 B因素离差平方和,各项离差平方和与自由度的计算公式为:,误差平方和 SSe=SST-SSA-SSB 总自由度 dfT=ab-1 A因素自由度 dfA=a-1 B因素自由度 dfB=b-1 误差自由度 dfe= dfT - dfA dfB =(a-1)(b-1),相应均方为,【例1】某厂现有化验员3人，担任该厂牛奶酸度（T）的检验。每天从牛奶中抽样一次进行检验，
4、连续10天的检验分析结果见表6。试分析3名化验员的化验技术有无差异，以及每天的原料牛奶酸度有无差异（新鲜牛奶的酸度不超过20 T ） 。,表1 牛奶酸度测定结果,A因素（化验员）有3个水平，即a=3；B因素（天数） 有10个水平 ，即 b =10 ， 共有ab=310=30个观测值。,1 计算各项离差平方和与自由度,2 列出方差分析表，进行F检验,表2 资料的方差分析表,结果表明，3个化验员的化验技术没有显著差异，不同日期牛奶的酸度有极显著差异。,注：F0.01(9,18）=3.60,3 多重比较 在两因素无重复观测值试验中，A因素每一水平的重复数恰为B因素的水平数b，故A因素的标准误为 ；同
5、理，B 因 素 的 标准误,对例1分析，a=3，MSe=0.0258。故,根据 dfe=18，秩次距 k=2，3 ，10，查临界 q 值 ，计算最小显著极差LSR，见表3。,表3 q值与LSR值,B因素各水平均值多重比较结果见表4,表4不同测定日牛奶酸度多重比较结果（q法）,附表：多重比较结果字母表示,结果表明，除B2与B5，B1与B9，B4与B3，B8与B3、B4，B10与B3、B4、B8差异不显著外，其余不同测定日间牛奶酸度均差异极显著或显著。酸度最高的是B7，最低的是B5和B2。从牛奶质量要求看，连续10d的牛奶酸度均在鲜奶范围内。,在进行两个因素或多个因素的试验时，除了要研究每一个因素
6、对试验指标的影响外，往往更希望知道因素之间的交互作用对试验指标的影响情况。 通过研究环境温度、湿度、光照、气体成分等环境条件对导致食品腐烂变质的酶和微生物的活动的影响有无交互作用，对有效控制酶和微生物活动，保持食品质量有着重要意义。,两个因素无重复观测值试验只适用于两个因素间无交互作用的情况； 若两因素间有交互作用， 则每个水平组合中只设 一个试验单位(观察单位)的试验设计是不正确的或不完善的。这是因为：,(1)在这种情况下，SSe,dfe实际上是A、B 两因素交互作用平方和与自由度，所算得的MSe是交互作用均方 ，主要反映由交互作用引起的变异。 (2)这时若仍按前述方法进行方差分析，由于误差
7、均方值大（包含交互作用在内），有可能掩盖试验因素的显著性， 从而增大犯型错误的概率。 (3) 每个水平组合只有一个观测值，无法估计真正的试验误差，因而不可能对因素的交互作用进行研究。,交互作用 交互作用：在多因素试验中一个因素对试验结果的影响依赖于另一因素所取的水平时，称两因素有交互作用。 在多因素对比试验中，某些因素对指标的影响往往是互相制约、互相联系的。即在试验中不仅因素起作用，而且因素间有时联合起来起作用，这种联合作用并不等于各因素单独作用所产生的影响之和，称这种联合作用为交互作用。 例：某农场对四块大豆试验田作施肥试验。每块田以不同的方式施以磷肥和氮肥，其产量如下：,可以看出 当施氮肥
8、和不施氮肥时，施以4公斤磷肥后的增产数量是不同的 当施磷肥和不施磷肥时，施以6公斤氮肥后的增产数量是不同的 若N, P分别起作用时增产为50, 30kg。但同时施时其效果并不是50+30=80kg，而是增产560-400=160kg，增加的80公斤则为交互作用的效果。,对两因素和多因素等重复试验结果进行分析， 可以研究因素的简单效应、主效应和因素间的交互作用（互作效应）。,3.1.2 交叉分组两因素等重复试验的方差分析,三种效应,1简单效应（simple effect） 是指在某一因素同一个水平上，比较另一因素不同水平对试验指标的影响。,三种效应,2主效应（main effect） 是指某一因
9、素各水平间的平均差别。它与简单效应的区别是，主效应指的是某一因素各水平间的平均差别是综合了另一因素各水平与该因素每一水平所有组合的情况。,三种效应,3. 互作效应（interaction effect） 如果某一因素的各简单效应随另一因素的水平变化而变化，而且变化的幅度超出随机波动的程度，则称两个因素间存在互作效应。,设A、B两因素，A因素有a个水平，B因素有b个水平，共有ab个水平组合，每个水平组合有n次重复试验，则全试验共有abn个观测值。试验结果的数据模式如表5所示。,两因素等重复试验的方差分析,表5 两因素等重复观测值试验数据模式,两因素等重复试验数据模式（部分）,表5中,每个组合处理
10、n 次重复之和,B因素第j水平an个数据之和,abn个数据总和,A因素第i水平bn个数据之和,其中， 为总平均数； i为Ai的效应； j为Bj的效应； () ij为Ai与Bj的互作效应。,(5-32),两因素等重复试验资料的数学模型为:,分别为Ai、Bj、Ai Bj观测值总体平均数；且,3.1.2.1 离差平方和与自由度分解,其中，SSAB，dfAB为A因素与B因素交互作用平方和与自由度。,为随机误差，相互独立，且服从N（0，2）。,若用SSAB，dfAB表示A、B水平组合间的平方和与自由度，即处理间平方和与自由度，则处理引起的变异可进一步剖分为A因素、B因素及A、B交互作用三部分，于是SSA
11、B、dfAB可分解为：,矫正数,总平方和与自由度,因素水平组合平方和与自由度,A因素平方和与自由度,各项平方和、自由度及均方的计算公式如下：,B因素平方和与自由度,交互作用平方和与自由度,误差平方和与自由度,所以，相应均方为,因素A的方差,因素B的方差,A、B互作的方差,误差方差,3.1.2.2 列方差分析表，进行F检验,FA显著，应对A因素各水平的平均数作多重比较，其平均数标准误为：,FB显著，应对B因素各水平的平均数作多重比较，其平均数标准误为：,FAB显著，应对各组合的平均数作多重比较，其平均数标准误为：,3.1.2.3 多重比较,表6 3种温度对3种pH值对酶活性的影响,【例2】现有3
12、种温度3种pH值对酶活性的影响试验结果，试作方差分析,A因素（温度）有3个水平，即a=3；B因素（pH值）有3个水平，即b=3；共有ab=33=9个水平组合；每个水平组合重复数n=3；全试验共有=333=27个观测值。,（1） 计算各项平方和与自由度,表7 方差分析表,（2）列出方差分析表，进行F检验,查临界F值： F0.05（2,18）=3.55，F0.01（2,18）=6.01； F0.01（4,18）=4.58。 因为， FAF0.05(2,18）； FBF0.05(2,18）；FABF0.01（4,18），表明不同温度、pH值与温度之间的互作对酶活性有显著或极显著影响，而pH间的差异不
13、显著。因此，应进一步进行不同温度处理均数间、各水平均数间 的多重比较。, 温度 因为A因素各水平的重复数为bn，故A因素各水平的标准误为： 对本例而言，,（3）多重比较,由dfe=18，秩次距k=2,3，从附表5中查出SSR0.05与SSR0.01的 临 界值 ，计算LSR值 ，结果列于表8。,表8 配方各水平均数比较SSR值与LSR值,表9 配方间平均数多重比较结果（SSR法）,因素A主效应分析，结果表明温度A3与A1之间差异极显著，A2与A1差异显著，A2与A3差异不显著。,因B因素各水平的重复数为an，故B因素各水平的标准误为：,在本例，B因素的影响不显著，不必进行多重比较。,以上所进行
14、的多重比较，实际上是A、B两因素主效应的检验。若A、B因素交互作用不显著，则可从主效应检验中分别选出A、B因素的最优水平，得到最优水平组合；若A、B因素交互作用显著，则应进行水平组合平均数间的多重比较，以 选出最优水平组合，同时可进行简单效应的检验。,因为水平组合数通常较大(本例ab=33=9)，采用最小显著极差法进行各水平组合平均数的比较，计算较麻烦。为了简便起见，常采用LSD法。,因为水平组合的重复数为n，故水平组合的标准误为： 本例, 各水平组合平均数间的比较,表15个水平组合平均数多重比较结果（SSR法）,分析结果表明，A3B3，A2B1，A1B1为优组合，按此组合选用温度和pH值可望
15、得到较好的酶活性。,以上的比较结果可以看出，当A、B因素的交互作用显著时，一般不必进行两个因素主效应的显著性检验（因为这时主效应的显著性在实用意义上并不重要），而直接进行各水平组合平均数的多重比较，选出最优水平组合。,简单效应的检验 简单效应实际上是特定水平组合平均数间的差数。检验尺度仍为互作效应检验中的LSD,有人设计3个罗非鱼品种A1、A2、A3和A4不同蛋白质水平饵料B1、B2、B3，每个处理配置两个鱼池进行试验。试验期内每池的产鱼量（kg）如下表。试作方差分析。,试验期内的产鱼量（kg）,解:1、提出假设, 2、平方和、自由度的分解,(2)、 方差分析表,因为互作不显著，则将互作合并到
16、误差项中,(2)、 方差分析表,（3）对B因素进行比较多重比较（SSR法）,多重比较表（字母标记法）,3.2 系统分组资料的方差分析,系统分组（hierarchical classification）设计，又称为巢式设计、树状设计、多因子嵌套设计 假设有 A、B 两个因子，这两个因子的搭配组成不再是上一节的交叉构成，而是 B 因子嵌套在 A 因子内，即 B 因子为次级因子： A 因子的某一个水平包含了 B 因子部分水平 A 因子的另一个水平包含了 B 因子的另一部分水平 即：B 因子的水平仅从属于 A 因子的一个水平 而 A 因子的水平并不包含 B 因子的所有水平,如果有第三个因子 C，则 C 因子嵌套在 B 因子内 如果有第四个因子 D，则 D 因
