1、第8章 正交试验设计的方差分析前面我们讨论了如何安排正交试验以及用极差分析法(即直观分析法)对试验结果进行计算分析.极差分析法简单明了，通俗易懂，计算工作量少，便于普及推广.但这种方法不能把试验中由于试验条件的改变引起的数据波动，同试验误差引起的数据波动区分开来.也就是说，不能区分因素各水平对应的试验结果间的差异，究竟是由于因素水平不同引起的，还是由于试验误差引起的，即不知道试验的精度.同时，对影响试验结果的各个因素的重要程度，既不能给出精确的定量估计，也不能提供一个标准，用来判断所考察的因素的作用是否显著.为了弥补极差分析法的不足，对试验结果的分析可采用方差分析法.8.1 正交试验方差分析的
2、基本步骤在第2章中我们已经介绍过，方差分析的基本思想是将数据的总偏差平方和(ST)分解为因素的偏差平方和(SA、SB)和误差的偏差平方和(Se)，然后将偏差平方和除以相对应的自由度(f)得到方差(VA、VB)，最后利用因素方差与误差方差之比(VA/Ve，VB/Ve)，作F检验，即可判断因素的作用是否显著.正交试验设计的方差分析也是按这样的步骤进行的，所不同的是这是考虑的是多因素试验的方差分析，而第2章中只考虑单因素和双因素试验的方差分析.一、计算1. 偏差平方和与自由度的计算方差分析的关键是偏差平方和的分解，现在以最简单的L4(23)正交表上安排的试验为例来说明(见表8-1，板书).不考虑哪些
3、因素安排在哪些列上(即表头设计时)，设试验结果为x1、x、x3和x4.总的偏差平方和: T=(x+x+x+x)-(x)2 整理后可得 () ()第1列各水平偏差平方和为S1=2=2=2= =表8-1 L4(23)正交表及计算表列号试验号123试验数据1234 1 1 2 2 1 2 1 2 1 2 2 1 x1 x2 x3 x4K1jK2jK11=x1+x2K21=x3+x4K12=x1+x3K22=x2+x4K13=x1+x4K23=x2+x3T=x1+x2+x3+x4注: Kij表示第j列第i水平的指标值之和；表示第j列第i水平的平均指标值；T表示指标值总和；表示平均指标值.同理，第2、3
4、列各水平的偏差平方和S2、S3为由此可得 ST=S1+S2+S3 (8-1)式(8-1)是正交表L4(23)的总偏差平方和的分解公式，即L4(23)的总偏差平方和等于各列偏差平方和之和.若在L4(23)正交表的第1列和第2列分别安排二水平因素A、B，在不考虑A、B因素间交互作用的情况下，则第3列(空列)是误差列.同样也可以证明ST=SA+SB+Se (8-2)上式也是总偏差平方和的分解公式，即总偏差平方和等于各列因素的偏差平方和与误差的偏差平方和之和.我们可以把上例推广到一般情况:用饱和正交表Ln(mk)安排试验(见表8-2，p160)，总的试验次数为n，每个因素的水平数为m，则每个水平作r次
5、试验，r=. 试验结果为x1，x2，x3，xn.令 则总偏差平方和为 (8-3)列偏差平方和为 (8-4)其中特别地， 当m=2(即二水平)时， 式(8-4)可表示成: (8-5)列偏差平方和Sj是第j列中各水平对应的试验数据平均值与总平均值的偏差平方和，它反映了该列水平变动所引起的试验数据的波动.若该列安排的是因素，就称Sj为该因素的偏差平方和；若该列安排的是交互作用，就称Sj为该交互作用的偏差平方和；若该列为空列，则Sj表示由于试验误差和未被考察的某些交互作用或某些条件因素所引起的波动.在正交试验设计中，通常把空列的偏差平方和作为试验误差的偏差平方和，虽然它属于模型误差，一般比试验误差大（
6、当作安全系数考虑），但用它作为试验误差进行显著性检验，可使检验结果更可靠些。总偏差平方和的自由度： fT=n-1 第j列偏差平方和的自由度： fj=mj-1 （mj：第j列水平数）此外，可以证明： ST=+ （8-6） fT=+ （8-7）式中k因、k交和k空分别为试验因素、试验考察的交互作用和空列在正交表中所占的列数。并且 k=k因+k交+k空注意：（1）当某个交互作用占有正交表的某几列时，该交互作用的偏差平方和就等于所占各列偏差平方之和，其自由度也等于所占各列的自由度之和； （2）误差的偏差平方和（Se）等于所占有空列的偏差平方和之和， 其自由度等于所有空列的自由度之和，即：， （3）上面
7、讨论的虽然是等水平饱和正交表Ln(mk)的情况，但是对于饱和的混合型正交表Ln（m1k1m2k2）也适用，不过要换上相应的m和k值。（有关“饱和正交表”与“不饱和正交表”的概念，请参见第6节（p191）！）2、方差的计算（V因，V交，Ve）方差等于各偏差平方和除以相应的自由度，即平均偏差平方和。V因= V交= Ve= 二、显著性检验数学上可以证明：在“假设H0：某因素或某交互作用对试验结果影响不显著”成立时，统计量F （8-9）服从第一自由度为f因（或f交），第二自由度为fe的F分布。对于给定的显著性水平，查F分布表得临界值F，若计算出的F值F0F 则拒绝原假设，认为该因素或该交互作用对试验的
8、结果有显著影响；若计算出的F值F0F，则接受原假设，认为该因素或该交互作用对试验结果无显著影响。经显著性检验后，可把检验结果列出方差分析表，如表8-3所示表8-3 正交试验方差分析表方差来源偏 差平方和自由度方 差F值F显著性 A BABSA=S1SB=S2SAB=S3fA=m-1fB=m-1fAB=fAfBVA=SA/fAVB=SB/fBVAB=SAB/fAB FA=VA/VeFB=VB/VeFAB=VAB/Ve查表误差eSefeVe=Se/fe 总 和STfT=n-1在进行正交试验方差分析时，应注意以下几点：（1）进行F检验时，要用到Se和fe，而Se=，fe=所以，为进行方差分析，选正交
9、表时应留出一定的空列。当无空列时，则应进行重复试验，以便求得Se2的值(见p181第5节)。（2）误差自由度fe一般不应小于2，即fe2，否则F检验的灵敏度很低，有时即使因素对试验指标有影响，用F检验也判断不出来。（3）如果fe1，为了增大fe，提高F检验的灵敏度， 在进行显著性检验之前，先比较V因和V交与Ve之间的差异程度。如果与误差方差Ve的大小相近，说明该因素或该交互作用对试验结果的影响微乎其微，其偏差平方和是由于随机误差引起的。因此，可并入误差偏差平方和Se中.通常把满足 V因(或V交)2Ve的那些因素或交互作用的偏差平方和，并入误差的偏差平方和Se中，从而得到新的误差偏差平方和S，相
10、应的自由度也并入fe中，从而得到f，然后用 (8-10)其中，校正后的误差方差为对其他因素或交互作用进行检验，这样使自由度fe扩大到f， 故可以提高F检验的灵敏度.三、最优条件的确定根据显著性检验结果，可以确定各因素对试验指标影响的主次顺序(通常根据F的大小判断)，对于显著性因素，若不考虑交互作用或交互作用不显著，则可通过比较该因素各水平对应的数据和(Kij)的大小，确定最优水平，各因素的优水平组合即为该试验的最佳水平组合，即最优条件；对于交互作用显著的某二个因素，必须先通过比较这二个因素各水平组合下试验数据之和的大小(即相当于极差分析中的二元表或搭配表)，然后再确定其最佳水平组合.最后，在最
11、优工艺条件下进行验证实验。8.2 不考虑交互作用的等水平正交试验方差分析8.2.1 二水平正交试验的方差分析（因学时有限，不讲解！只讲解三水平情况，因为三水平会，二水平自然就会！）例8-1 在双歧杆菌酸奶研制中，为选择最佳发酵条件，用L8(27)正交表安排了正交试验，试验因素与水平表见表8-4，试验方案及结果见表8-5。试对试验结果进行方差分析。表8-4 试验因素水平表因素水平葡萄糖A1#生长促进剂 B接种量C厌氧处理D基质浓度 E2#生长促进剂 F试验指标12201035充氮气不充气10120.250活菌数的对数是六因素二水平试验，且不考虑交互作用，用L8(27)最好！表8-5 试验方案及结
12、果分析表头设计ABCDEF试验数据列号试验号1234567123456781(2)1112(0)2221(1)12(0)211221(3)12(5)22211121212121(充氮气)2(不充气)1221211(10)2(12)2112211(0.25)2(0)2121127.5802.4772.6997.5682.4777.5316.6022.00057.4566.1367.28557.2756.13656.71643.5864.000K1jK2j20.32418.61020.06518.86918.65920.27519.35819.57619.81019.12419.62519.309
13、29.2819.653T=38.934K1j-K2jSj1.7140.3671.1960.179-1.6160.326-0.2180.005940.6860.05880.3160.012519.62848.157一、计算1、计算各列各水平的Kij值（K1j和K2j）各列各水平的试验数据（即指标值）之和K1j、K2j以及（K1j-K2j）的值，填入表8-5中, 如K1c=7.580+2.477+6.602+2.000=18.659 K2c=2.699+7.568+2.477+7.531=20.2752、计算各列的偏差平方和Sj及其自由度fj由式（8-5）知，Sj= (K1j-K2j)2 SA=S
14、1=1.7142=0.367 SB=S2=1.1962=0.179 SC=S3=(-1.616)2=0.326 Se=S4=(-0.218)2=0.00594 fj=m-1=2-1=1 fT=n-1=8-1=7 fe=fj=f4=1将求得的Sj值也填入表8-5中。为了判断是否计算有误，进行以下验算： ST的验算T= =38.934，CT= =38.9342=189.482QT=57.456+6.136+4.000=238.589ST=QT-CT=238.589-189.482=49.107另外 ST=SA+SB+Se=0.367+0.179+0.00594=49.107 fT的验算fT= =f
15、A+fB+fe=1+1+1=7另外 fT=n-1=8-1=7ST和fT均计算无误。3、计算方差VA=SA/fA=0.367/1=0.367 VB=SB/fB=0.179/1=0.179 Ve=Se/fe=0.00594/1=0.00594 注意： fe=12， F检验的灵敏度低！ 需要校正fe- f、Se-S、Ve-V二、显著性检验根据上述计算结果，进行显著性检验，列出方差分析表，如表8-6所示. VE=0.0125最小， VE/Ve=0.0125/0.00594=2.10因素E对试验指标的影响可忽略，故将其偏差平方和SE并入误差平方和Se中，即S=Se+SE=0.00594+0.0125=0
16、.01844， f=fe+fE=1+1=2V=S/f=0.01844/2=0.009221. 计算FjFA=VA/V=0.367/0.00922=39.8， FB=VB/V=0.179/0.00922=19.4FC=VC/V=0.326/0.00922=35.4， FD=VD/V=0.0588/0.00922=6.38FF=VF/V=48.157/0.00922=5223.12. 查FF(f，f)= F(1，2)当=0.05时，查F分布表(p324)得F0.05(1，2)=18.51当=0.01时，查F分布表(p326)得F0.01(1，2)=98.503. 显著性检验 FFF0.01（1，2
17、） F因素高度显著（用*表示）；又 F0.05(1，2)FAF0.01(1，2)，F0.05(1，2)FBF0.01(1，2)，F0.05(1，2)FCF0.01(1，2)， 因素A、B、C均为显著(用*表示)；又 FD F0.05(1，2)， D因素不显著(不用记号表示).由 F值大小可知，各因素对试验指标影响的主次顺序为:F A C B D E4.列出方差分析表表8-6 方差分析表方差来源偏差平方和自由度方差F值F显著性ABCDEF误差e误差e0.3670.1790.3260.05880.012548.1570.005940.0184111111120.3670.1790.3260.058
18、80.012548.1570.005940.0092239.819.435.46.385223.1F0.05(1,2)=18.51F0.01(1,2)=98.50*总 和49.1077三、最佳条件的确定.本例中，指标值越大越好，由K1j和K2j的大小，按各因素的主次顺序选优水平如下：F 选F1 ，A 选A1，C 选C2，B 选B1.对于不显著因素D 和E，可视具体情况而定，D 选D2(即不充气)为好，E选E1(10%)可降低成本。所以，最优工艺条件为A1 B1 C2 D2 E1 F1 (或 A1 B1 C2 E1 F1)，即最后，在最优工艺条件下进行验证实验。8.2.2 三水平正交试验的方差分
19、析（重点讲解！三水平掌握了，二水平自然就会！）例8-2 自溶酵母提取物是一种多用途食品配料。为探讨外加中性蛋白酶方法中啤酒酵母的最适合自溶条件，安排了三因素三水平试验。试验指标是自溶液中蛋白质含量（）。试验因素水平表见表8-7，试验方案及试验结果见表8-8。试对试验结果进行方差分析。表8-7 因素水平表因素水平温度（）ApH值B加酶量C1235055586.57.07.52.02.42.8本例是三因素三水平试验(见表8-7)，且不考虑因素间的交互作用，所以用L9(3 4)正交表安排试验最合适.试验方案和试验结果分析见表8-8.表8-8 试验方案及结果分析表头设计ABC试验指标Pr 列号试验号1
20、2341234567891(50)112(55)223(58)331(6.5)2(7.0)3(7.5)1231232.41(2.0)2(2.4)3(2.8)2313121233122316.254.974.547.535.545.5011.410.98.95K1jK2jK3j15.7618.5731.2525.1821.4118.9922.6521.4521.4820.7421.8722.97T=65.58248.38344.84976.56634.03458.39360.62513.02460.10461.39430.15478.30527.62一、 计算1. 计算各列各水平的Kij值(K1
21、j、K2j 和K3j)各列各水平对应的试验指标之和K1j、K1j 和K1j及其平方K、K和K，列于表8-8中。例如：K1C=6.25+5.50+10.9=22.65， K1C2=22.652=513.02K2C=4.97+7.53+8.95=21.45， K2C2=21.452=460.10K3C=4.54+5.54+11.4=21.48， K3C2=21.482=461.392. 计算各列的偏差平方和Sj 及其自由度fj根据式(8-4)可知， Sj=-CT， r=n/m=9/3=3CT=T2/n=65.582/9=477.86SA=S1=(K+ K+ K)/r-CT=(248.38+344.
22、84+976.56)/3-477.84=45.4同理可知，SB=S2=6.49， SC=S3=0.31，Se=S4=0.83fj=m-1， fA=fB=fC=fe=3-1=2验算:1. QT=6.252+4.972+8.952=530.89 ST= QT-CT=530.89-477.86=53.03另外， ST= =SA + SB + SC + Se=45.4+6.49+0.31+0.83=53.03fT=n-1=9-1=8 另外， fT= =fA +fB + fC +fe =24=8计算正确无误。3. 计算方差VA= SA /fA=45.4/2=22.7， VB= SB /fB=6.49/2
23、=3.25， VC= SC /fC=0.31/2=0.155， Ve= Se /fe=0.83/2=0.415， VC2 Ve 因素C的偏差平方和SC是由于随机误差引起的，说明因素C对试验结果的影响可忽略，故将SC并入Se中，得Se=Se+Sc=0.83+0.31=1.14fe=fe+fc=2+2=4， Ve= Se/fe=1.14/4=0.285二、显著性检验 1计算Fj 22.7/0.285=79.63.25/0.285=11.42查F 当=0.05时，查F分布表得F(f因，f)=F0.05(2，4)=6.94 当=0.01时，同理查得F0.01(2，4)=18.00.3.显著性检验FAF
24、0.01； 因素A高度显著，用*表示之。又F0.05FBF0.01； 因素B显著，用*表示之。 当然，因素C是不显著的。因素作用的主次顺序为：A B C4方差分析表表8-9 方差分析表方差来源偏差平方和自由度方差F值F显著性ABC误差e误差e45.46.490.310.831.142222422.73.250.1550.4150.28579.611.4F0.05(2,4)=6.94F0.01(2,4)=18.0*总 和53.038三、最优工艺条件本例试验的指标值越大越好。对因素A和B，通过比较Kij可知，优水平为 A3和B1，对因素C，取C1有利于节省原料。故最优水平组合为A3B1C1，即最后
25、，在最优工艺条件下进行验证实验。因为因素 C对试验指标几乎无影响，且因素C为加酶量，而酶的价格通常较高。所以最好在加酶量更低时，再次进行正交试验，以便确定最经济合理的加酶量。8.3 考虑交互作用的正交试验设计的方差分析8.3.1 考虑交互作用的二水平正交试验的方差分析（重点讲解！二水平掌握了，三水平自然也会！） 如前所述，因素间的交互作用在多因素试验设计中是经常碰到的，在一般试验中，在采取一些措施后，多数交互作用可略去。但为了使试验做得更精确，表头设计得更合理，收到更好的效果，在正交试验设计的方差分析中也要考虑因素间的交互作用。现在实例说明之。例8-3 试对例7-4（p150154）的试验结果
26、进行方差分析。表8-10 试验方案及结果分析试验号ABABCACBC吸光度xi1234567123456781(300)1112(700)2221(1800)12(2400)21122112222111(8)2(10)1212121212212112211221122121122.422.242.662.582.362.402.792.76K1K29.910.319.4210.7910.2110.010.239.9810.249.9710.1210.0910.1910.02=20.2198.01106.3088.74116.42104.24100.00104.6599.60104.8699.4
27、0102.41101.81103.84100.40Sj0.0210.02350.005510.007810.009110.0001130.00361一、计算书上为了便于计算，将试验数据x。放大10倍，即令x=10xi，变换后的数据如表8-10所示。（这种变换实际上是没有必要的，因为大家现在都用计算器或是计算机进行计算，而不是用对数计算尺进行计算！） 1.计算各列各水平的Kij值（K1j，K2j） 各列各水平对应的试验数据之和 K1j和K2j，以及（K1j-K2j）列于表8-10中2.计算各列的偏差平方和（Sj）和自由度（fj）对于二水平正交试验，由式（8-5）可知Sj=Qj-CT=- 或Sj=
28、(K1j-K2j)2 , r=n/m=8/2=4, SA=S1=0.021；同理 SB=S2=0.235； SAB=S3=0.00551 SC=S4=0.00781 SAC=S5=0.00911， SBC=S6=0.000113， Se=S7=0.00361fj=m-1=2-1=1 即fA=fB=fAB=fC=fAC=fBC=fe=2-1=1验算： ST的验算CT=20.212=51.056QT=2.422+2.242+2.762=51.337ST=QT-CT=51.337-51.056=0.281另外ST=0.021+0.235+0.00361=0.281 fT的验算fT=n-1=8-1=7
29、另外fT=fj=71=7 计算无误。3计算方差VA=SA/fA=0.021/1=0.021， VB=SB/fB=0.235/1=0.235，VAB=SAB/fAB=0.00551/1=0.00551， VC=SC/fC=0.00781，VAC=SAB/fAC=0.00911， VBC=SBC/fBC=0.00013，Ve=Se/fe=0.00361， fe=12，F检验的灵敏度低！VAB/Ve=0.00551/0.00361=1.53， VC/Ve=0.00781/0.00361=2.16VAC/Ve=0.00911/0.00361=2.52， VBC/Ve=0.000113/0.000361
30、=0.0313VAB2Ve， VBC2Ve交互作用AB和BC对试验指标的影响可以忽略， SAB和SBC主要是由于随机误差而引起的。因此，将SAB和SBC并入Se中，得S=Se+SAB+SBC =0.00361+0.00551+0.000113 =0.00923f=fe+fAB+fBC=1+1+1=3V=S/f=0.00923/3=0.00308二、显著性检验1. 计算FjFj=Vj/V， j=A、B、C、ACFA=VA/V=0.021/0.00308=6.82， FB=VB/V=0.235/0.00308=76.3FC=VC/V=0.00781/0.00308=2.54， FAB=VAC/V
31、=0.00911/0.00308=2.962.查Ff因（或f交），f=F(1，3)当=0.05时，查F分布表及F0.05(1，3)=10.13；当=0.01时，查F分布表及F0.01(1，3)=34.12.3.显著性检验FBF0.01 因素B高度显著（*）；又FA，FC，FAXC均小于F0.05，因素A和C及其交互AC对试验指标的影响均不显著。根据上述计算所D得到的方差分析表，如表8-11所示，F值的大小，可知各因素的主次顺序为：B、A、AC、C、AB、BC。表8-11 方差分析表方差来源偏差平方和自由度方差F值F显著性ABABCACBC误差e误差e0.0210.2350.005510.007
32、810.009110.0001130.003610.0923111111130.0210.2350.005510.007810.009110.0001130.003610.003086.9878.072.593.03F0.05(1,3)=10.13F0.01(1,3)=34.13*总 和0.2817三、最优条件确定因为交互作用AB和BC的影响可忽略，而AC的影响又不显著，所以确定因素的优水平时可不考虑交互作用的影响（注：当交互作用影响显著时，应根据因素之间的搭配表即二元表选取优水平）。通过比较试验数据和K1j和K2j可知，B应取B2水平，A应取A2水平。C为次要因素，可取C1或C2，视具体情况而定。这样，最优水平组合为A2B2C1或A2B2C2，即 最后，在最优工艺条件下进行验证实验。讨论：通过极差分析与方差分析的比较，可以看出方差分析的优点： （1）可以分析出试验误差的大小，从而知道试验的精度。 （2）不仅可以给出各因素及交互作用对试验指标的影响的主次顺序，而且可以分析出哪些因素的影响显著，哪些因素的影响不显著。因为主要因素不一定是显著因素，次要因素也不一定是不显著因素，所以，极差分析法不能对各种因素的主要程度给予精确的定量的估计。对于显著因素，我们在试验中取其最优水平，进行严格控制；对于不显著因素，可以视具体情况，综合考虑试验成本和操作的难易等方面确定其最适宜水平。
