双因素方差分析

1、双因素方差分析方法双因素方差分析方法l 双因素试验的方差分析双因素试验的方差分析 在实际应用中，一个试验结果（试验指标）往往在实际应用中，一个试验结果（试验指标）往往受多个因素的影响。不仅这些因素会影响试验结果，受多个因素的影响。不仅这些因素会影响试验结果，而且这些因素的不同水平的搭配也会影响试验结果。而且这些因素的不同水平的搭配也会影响试验结果。 例如：某些合金，当单独加入元素例如：某些合金，当单独加入元素A或元素或元素B时，时，性能变化不大，但当同时加入元素性能变化不大，但当同时加入元素A和和B时，合金性时，合金性能的变化就特别显著。能的变化就特别显著。 统计学上把多因素不同水平搭配对试验

2、指标的统计学上把多因素不同水平搭配对试验指标的影响称为交互作用。交互作用在多因素的方差分析影响称为交互作用。交互作用在多因素的方差分析中，把它当成一个新因素来处理。中，把它当成一个新因素来处理。 我们只学习两个因素的方差分析，更多因素的我们只学习两个因素的方差分析，更多因素的问题，用正交试验法比较方便。问题，用正交试验法比较方便。无交互作用的双因素试验的方差分析无交互作用的双因素试验的方差分析 数学模型数学模型 假设某个试验中，有两个可控因素在变化，因素假设某个试验中，有两个可控因素在变化，因素A有有a个水平，记作个水平，记作A1，A2，Aa；因素；因素B有有b个水平，个水平，记作记作B1，B

3、2，.Bb；则；则A与与B的不同水平组合的不同水平组合AiBj（i=1，2，a；j=1，2，b）共有）共有ab个，每个水平组合个，每个水平组合称为一个处理，每个处理只作一次试验，得称为一个处理，每个处理只作一次试验，得ab个观测个观测值值Xij，得双因素无重复实验表，得双因素无重复实验表双因素无重复（无交互作用）试验资料表双因素无重复（无交互作用）试验资料表因素因素 A12.bBBB1112112.baaabXXXXXX.1ajijiTX.1.2.bTTT.jjXTa.1biijjTX1XTab.1.2.bXXX因素因素 B1.aAA.iiXT b1.aTT1.aXX11abijijTX 无交

4、互作用的双因素试验的方差分析无交互作用的双因素试验的方差分析ijijijX线性统计模型线性统计模型 基本假设（基本假设（1） 相互独立；相互独立； （2） ，（方差齐性）。，（方差齐性）。ijX2,ijijXN 其中其中 111abijijab所有期望值的所有期望值的总平均总平均 11biijija11ajijjibijijijX水平水平Ai对试验结果的效应对试验结果的效应 水平水平Bj对试验结果的效应对试验结果的效应 试验误差试验误差 特性：特性： 2110; 0; 0,abijijijN11biijija11ajijjibijijijX水平水平Ai对试验结果的效应对试验结果的效应 水平水平

5、Bj对试验结果的效应对试验结果的效应 试验误差试验误差 要分析因素要分析因素A，B的差异对试验结果是否有显著的差异对试验结果是否有显著影响，即为检验如下假设是否成立：影响，即为检验如下假设是否成立：0112:0aH0212:0bH 总离差平方和的分解定理总离差平方和的分解定理仿单因素方差分析的方法，考察总离差平方和仿单因素方差分析的方法，考察总离差平方和211abTijijSSXX可分解为：可分解为：TABESSSSSSSS2.1aiAiSSbXX称为因素称为因素A的离差平方和，的离差平方和，反映因素反映因素 A 对试验指标的影响。对试验指标的影响。2.1bjBjSSaXX称为因素称为因素B的

6、离差平方和，的离差平方和，反映因素反映因素 B 对试验指标的影响。对试验指标的影响。2.11abijEijijSSXXXX称为误差平方和，反映试验误差对试验指标的影响。称为误差平方和，反映试验误差对试验指标的影响。2,ijXN 可推得：可推得：221TSSab2222,TABESSSSSSSS将将 的自由度分别记作的自由度分别记作,TABEdfdfdfdf，则，则 1 ,11AAAAEEESSdfMSFFaabSSdfMS2211ESSab 1 ,11BBBBEEESSdfMSFFbabSSdfMS若假设若假设 成立，则：成立，则：0102,HH221ASSa221BSSb对给定的检验水平对给

7、定的检验水平 ，F 右侧检验右侧检验 1 ,11AAAAEEESSdfMSFFaabSSdfMS 1 ,11BBBBEEESSdfMSFFbabSSdfMS 1 ,11AFFaab时，时，当当 1 ,11BFFbab时，时，当当拒绝拒绝H01，即，即A 因素的影响有统计意义。因素的影响有统计意义。拒绝拒绝H02，即，即B 因素的影响有统计意义。因素的影响有统计意义。双因素（无交互作用）试验的方差分析表双因素（无交互作用）试验的方差分析表方差来源方差来源因素因素A总和总和平方和平方和ASSBSSTSS自由度自由度AdfEdfTdf均方和均方和AAASSMSdfEEESSMSdfF 值值AAEMS

8、FMSF 值临介值值临介值(1 ,11 )Faab因素因素B误差误差ESSBdfBBBSSMSdfBBEMSFMS(1 ,11 )Fbab,ETABETABdfdfdffSSSSSSSS注意注意 各因素离差平方和的自由度为水平数减一，总平方各因素离差平方和的自由度为水平数减一，总平方和的自由度为试验总次数减一。和的自由度为试验总次数减一。双因素（无交互作用）试验的方差分析表双因素（无交互作用）试验的方差分析表,AABBSSDp SSDp简便计算式：简便计算式：其中：其中：2.1,aAiiDTb,EABTSSRDDp SSRp2.1,bBjjDTa2,pTab211abijijRX例例1 设甲、

9、乙、丙、丁四个工人操作机器设甲、乙、丙、丁四个工人操作机器、各一天，各一天， 其产品产量如下表，问工人和机器对产品产量是否有显著其产品产量如下表，问工人和机器对产品产量是否有显著 影响？影响？工人工人 A506352475442475741535848.1ajijiTX.jjXTa.1biijjTX51X 49.3 58.0 45.8机器机器 B.iiXT b612T 甲甲 乙乙 丙丙 丁丁 197 232 18316514314515955.047.748.353.0解解 基本计算如原表基本计算如原表 21131678abijijRX2.1123495aAiiDTb2.1142040.67b

10、BjjDTa231212Tpab466TSSRp114.67AASSDp318.5BBSSDp32.83ETABSSSSSSSS1 11Tdfn 1 3Adfa 12Bdfb 6EAbdfdfdf38.223AAAMSSSdf159.25BBBMSSSdf5.47EEEMSSSdf6.98AAEFMSMS29.10BBEFMSMS0.013,69.78F0.012,6BFF0.053,64.76F0.050.013,63,6AFFF0.012,610.92F 结论：工人对产品的产量有显著影响，结论：工人对产品的产量有显著影响，机器对产品的产量有极显著影响。机器对产品的产量有极显著影响。例例2：

11、某厂对生产的高速钢铣刀进行淬火工艺试验，考察回火温度：某厂对生产的高速钢铣刀进行淬火工艺试验，考察回火温度A和淬火温度和淬火温度B两个两个因素对强度的影响。今对两个因素各因素对强度的影响。今对两个因素各3个水平进行试验，得平均硬度见表：个水平进行试验，得平均硬度见表：BjAi试验结果试验结果B1（1210C）B2（1235C）B3（1250C）A1（280C） 64 66 68 A2（300C） 66 68 67A3（320C） 65 67 68假设：不同组合水平下服从正态分布、互相独立、方差相等。假设：不同组合水平下服从正态分布、互相独立、方差相等。所需要解决的问题是：所有所需要解决的问题是

12、：所有Xij的均值是否相等。的均值是否相等。方差分析表：方差分析表：方差来源方差来源 离差平方和离差平方和 自由度自由度 F值值 F0.05(2,4) F0.01(2,4) 显著性显著性因素因素A 1.56 2 FA=1.01 6.94 18.0因素因素B 11.56 2 FB=7.46 6.94 18.0 *试验误差试验误差 3.1 4总误差总误差 16.22 8 31312222.1633ijijTXSST56.133313122.iiTTSSA3122.56.113331jjTTSSB1 . 3SSBSSASSTSSE)4,2(05.0FFAA影响不显著。影响不显著。)4 , 2()4

13、, 2(01. 005. 0FFFBB影响显著，由于影响显著，由于高速钢洗刀的硬度越大越好，因此因素高速钢洗刀的硬度越大越好，因此因素B可取可取B3水平，即淬火温度水平，即淬火温度1250C为好，因素为好，因素A水平的确定，应考虑经济方便，取水平的确定，应考虑经济方便，取A1水平为好。水平为好。不同品牌的彩电在各地区的销售量数据不同品牌的彩电在各地区的销售量数据 品牌品牌(因素因素A) 销售地区销售地区( 因素因素B )B1B2B3B4B5A1A2A3A4365345358288 350368323280 343363353298 340330343260 323333308298 四、双因素

14、方差分析例题1、对因素A提出的假设为H0: 1=2=3=4 (品牌对销售量没有影响)H1: i (i =1,2, , 4) 不全相等 (品牌对销售量有影响)2、对因素B提出的假设为H0: 1=2=3=4=5 (地区对销售量没有影响)H1: j (j =1,2,5) 不全相等 (地区对销售量有影响)四、双因素方差分析例题差异源差异源SSdfMSFP-valueF crit 行行(因素因素A)13004.5534334.8518.10777 9.46E-053.4903 列列(因素因素B)2011.74502.9252.100846 0.1436653.2592 误差误差2872.712239.3

15、917 总和总和17888.9519例例1的上机操作的上机操作对应例对应例1 的数据输入方式的数据输入方式原始数据，行因素水平，列因素水平原始数据，行因素水平，列因素水平*0.010.0220.050.0010.01在在 下接受，在下接受，在 下否决下否决0.010.05在在 下否决下否决0.01（A）（B）工人对产品产量有显著影响，而机器对产品产量的影响工人对产品产量有显著影响，而机器对产品产量的影响极极显著。显著。 有交互作用的双因素试验的方差分析有交互作用的双因素试验的方差分析ijkijijkijX线性统计模型线性统计模型 基本假设（基本假设（1） 相互独立；相互独立； （2） ，（方差

16、齐性）。，（方差齐性）。ijX2,ijijXN 有检验交互作用的效应，则两因素有检验交互作用的效应，则两因素A，B的不同水的不同水平的搭配必须作重复试验。平的搭配必须作重复试验。 处理方法：把交互作用当成一个新因素来处理，处理方法：把交互作用当成一个新因素来处理，即把每种搭配即把每种搭配AiBj看作一个总体看作一个总体Xij。观测值观测值总平均总平均 因素因素A的效应的效应 因素因素B的效应的效应 交互作用交互作用的效应的效应 试验误差试验误差 有交互作用的双因素试验的方差分析有交互作用的双因素试验的方差分析线性统计模型线性统计模型 其中其中 111abijijab所有期望值的总平均所有期望值

17、的总平均 11biijija11ajijjibijkijkijX水平水平Ai对试验结果的效应对试验结果的效应 水平水平Bj对试验结果的效应对试验结果的效应 试验误差试验误差 ijkijijkijXijijij交互效应交互效应 特性：特性： 110; 0; abijij 要判断因素要判断因素A，B及交互作用及交互作用A B对试验结果是否对试验结果是否有显著影响，即为检验如下假设是否成立：有显著影响，即为检验如下假设是否成立：0112:0aH0212:0bH2110; 0; 0,abijkijijijN03:0 1,2, ;1,2,ijHia jb 总离差平方和的分解定理总离差平方和的分解定理仿单

18、因素方差分析的方法，考察总离差平方和仿单因素方差分析的方法，考察总离差平方和2111abnTijkijkSSXX可分解为：可分解为：TABA BESSSSSSSSSS SSA称为因素称为因素A的离差平方和，反映因素的离差平方和，反映因素 A 对试验对试验指标的影响。指标的影响。 SSB称为因素称为因素B的离差平方和，反映因素的离差平方和，反映因素 B 对试验指标的影响。对试验指标的影响。SSA B称为交互作用的离差平方和，称为交互作用的离差平方和，反映交互作用反映交互作用A B对试验指标的影响。对试验指标的影响。SSE称为误差平称为误差平方和，反映试验误差对试验指标的影响。方和，反映试验误差对

19、试验指标的影响。2,ijkXN 若若“各因素、各水平及其交互作用的影响无统计意义各因素、各水平及其交互作用的影响无统计意义”的假设的假设 成立，则成立，则 则则1,1AAAAEEESSdfMSFF aab nSSdfMS可推得：可推得：2222221 ,1 ,1ABTSSSSSSababn222211 ,1A BESSSSabab n1,1BBBBEEESSdfMSFF bab nSSdfMS11 ,1A BA BA BA BEEESSdfMSFFabab nSSdfMS 由由 作右侧假设检验来考察各因素及因素作右侧假设检验来考察各因素及因素间的交互作用对试验指标的影响力间的交互作用对试验指标的影响力.,ABA BFFF双因素有重复（有交互作用）试验资料表双因素有重复（有交互作用）试验资料表因素因素 A12.bBBB1111211 111121.bnnbnXXXXXX.因素因素 B1AaA1121112.aaaba na nabnXXXXXX双因素（有重