正交试验设计的方差分析例题汇总

1、8.3.2考虑交互作用的三水平正交试验的方差分析（因学时有限和正交表太大L27（313）,不讲解！只讲解二水平情况，因为二水平会，三水平自然也会！）例8-4运动发酵单细胞菌是一种酒精生产菌。为了确定其发酵培养基的最佳配方，进行了四因素三水平正交试验，试验指标为酒精浓度（g/ml）。表8-12给出了因素水平表，要求考察交互作用AxB、AXC和AXDo查附表7可得，本试验应选用L27（313）正交表，表头设计应按照“L27（313）二列间的交互作用表”进行。本例只考虑一级交互作用（p=1）,所以每个三水平交互作用应占（m-1）P=（3-1）1=2歹！J,即AXB、AXC,和AXD在L27（313）

2、正交表中各占二列。表8-12因素水平表因素水平'、葡锢糖浓度A酵母膏浓度%B培乔温度（C）C培加基pHD150255.02150.5306.03251.0357.0表头设计时应避免混杂，试验方案及试验结果见表8-13。由交互作用表可知，将因素A、B安排在第1、2列之后，第3、4歹U为AXB交互作用列；再将C安排在第5列后，AXC交互作用在第6、7歹I；最后将D安排在第9歹I，则AXD交互作用类落在第8、10列（当然也可将D安排在第8歹U,则第9、10列为AXD交互作用列）。表8-13试验方案及结果分析L27（313）试ABAXBCAXCAXDDAXD试验验号12345678910111

3、213结果V、111111111111110.20211112222222220.50311113333333330.50412221112223331.50512222223331111.10612223331112221.20713331113332221.60813332221113331.60913333332221111.201021231231231230.401121232312312310.501221233123123120.201322311232313126.301422312313121232.701522313121232314.201623121233122315.

4、901723122311233127.701823123122311236.151931321321321320.402031322132132130.302131323213213210.302232131322133211.752332132133211324.752432133211322135.302533211323212132.902633212131323217.302733213212131322.80Ku9.403.3032.7526.4019.9526.2022.6028.3016.6522.9019.7024.2022.4568.25Kj33.0527.8017.9024

5、.5526.4523.2018.8020.0023.4525.0022.4021.9024.45Kj25.8037.1517.6017.3021.8518.8526.8519.9528.1520.3526.1522.1521.35Kfj88.3610.891072.6696.96398.00685.44510.76800.89277.22524.41388.09585.64504.00K2j1092.3772.84320.41602.70699.60538.24353.44400.00549.90625.00501.76479.61597.80K3j665.641380.1309.76299.

6、29477.42355.32720.92398.00792.42414.12683.82490.62455.82Sj32.6267.9021.812.486.646.347.436.343.23一、计算(计算过程省略)1 .计算各列各水平的K值(Kij,K2j,K3j)和K2(K2j,K2j,K2j)各列各水平对应的试验数据之和Kij,K2j,K3j,及其平方和Kj,K2»K2j,列于表8-13中，例如9Ka=£Xi=0.20+0.50*2+1.50+1.10+1.20*2+1.60*2=9.40=Kii,i4K21=88.369Ka=zXi=0.40+0.50+6.15=

7、33.05=K21,=1092.30i19Ka=ZXi=0.40+0.30+2.80=25.80=K31,K3i=665.64i1表示AXB的有两列，即第3,4歹U,计算后可知K13=32.75,K23=17.90;K33=17.608=26.40;K24=24.55,K34=17.302 .计算各列的偏差平方和(Sj)及其自由度(fj)由式(8-4),可知:1m2T2-S=LKj_Qj_CTridnr=n/m=27/3=9;CT=T/n=1/27X68.252=172.52所以Sj=1/Kij2172.52=1(Kij2+Kj2+Kj2)-172.539Tj90=S=1(Kii2+Ki2+K

8、i2)-172.529i=；(88.36+i092.30+665.64)-i72.52=32.62Sb=S=67.90,S3=i6.67,8=5.i4所以3xB=S+S=2i.8iS=S=2.48,S6=3.04,S7=3.60所以3xc=S+S=6.64S8=5.i3;Si0=i.2i所以Saxd=S+S0=6.34&=S=7.43Si=.=2.33,S12=0.35,S13=0.55所以S=Si+S2+S3=3.23因为第j列的自由度为fj=m-1=3-1=2,(j=1,2,13),所以fA=fB=fC=fD=2fAXB=f3+f4=2+2=4,fAXC=f6+f7=2+2=4fA

9、XD=f8+f10=2+2=4,fe=f11+f12+f13=2X3=6验算：Sr的验算nQ="xi2=0.22+0.52+-+2.82=320.98i1St=Qt-CT=320.98-172.52=148.46k另外Sr=ZSj=32.62+67.90+0,55=148.45j=ifT的验算fr=n-1=27-1=26另外fr=13fi=13X(m-1)=13x(3-1)=26.计算过程无误.3 .计算方差V=SfjVA=S/fa=32.62/2=16.31同理可得V=33.95,Vc=1.24,Vd=3.72Vaxb=5.45,Vaxc=1.66,Vaxd=1.59Ve=0.53

10、8:业均大于2,且fe=6>1,.无需校正Ve!Ve:、显著性检验(计算过程省略)1 .计算FjFa=VA/Ve=16.31/0.538=30.32,.Fj=V/Ve,同理可得，Fb=63.10,Faxb=10.13,Fc=2.30,Faxc=3.092 .查FaFf因,fe)=f2,6),Ff交,fe)=F44,6)当a=0.05时，查得F0.05(2,6)=5.14,F0.05(4,6)=4.53;当a=0.01时，查得F0.01(2,6)=10.92,F0.01(4,6)=9.15.3 .显著性检验Fa=30.32和Fb=63.10均大于Fo.01（2,6）=5.14,因素所口融匀

11、高度显著（用*表示）；又Faxb=10.13>F0.0i（4,6）=9.15,交互作用AX他高度显著（用*表示）；又.Fd=6.91,.介于Fo.o5（2,6）和Fo.oi（2,6）之间，因素D显著（用*表示）；又Fc=2.30<Fo.o5（2,6）=5.14,以及Faxc=3.09和Faxd=2.96均小于Fo.05（4,6）=4.53,因素C及交互作用AX5f!AxD匀不显著.根据F值大小，可知各因素及交互作用对试验指标影响的主次顺序为：B,A,AXB,D,AXC,AXD,C.4 .列方差分析表表8-14方差分析表、.、.广.力左来源偏差平方和自由度、.、.广.力左F值巳显著性

12、A32.62216.5130.32F0.05(2,6)=5.14*B67.90233.9563.10F0.0i(2,6)=10.92*AXB21.8145.4510.13*C2.4821.242.30F0.05(4,6)=4.53AXC6.6441.663.09F0.0i(4,6)=9.15D7.4323.726.91*AXD6.3441.592.96误差e3.2360.538总和148.4526三、最优工艺条件确定因素A、B及交互作用AXB都高度显著，但因在主次顺序中，AXB排在A、B之后，因此应优先考虑AB的优水平。A和B的优水平确定了，其搭配也就随之确定，不必再通过A、B的二元表确定A与

13、B搭配。通过比较试验指标和K值大小，可知A和B的优水平，分别为A和因素D作用显著，但D与A的交互作用AXD不显著，故可不考虑交互作用，通过比较K值可知D的优水平为d因素C作用不显著，可以降低成本和操作方便等方面来考虑选取最适水平。对本例通过比较K值确定C的优水平为G。因此，最优水平组合为ARGR,即最优工艺条件为葡萄糖浓度15%酵母膏浓度1.0%、培养温度30c和培养基pH值7.0。最后，最好能在最佳条件ARGC3下，再实施一次试验，测定试验指标值(即酒精浓度)，在L27(313)正交表中，没有ARGC3这一组试验。在正常情况下，ABGQ组合条件下的试验指标值，应大于表8-13中的最大Xi值，

14、即第17号试验的xi7=7.70。8.4混合型正交表的方差分析混合型正交表的方差分析与等水平正交表的方差分析无本质的区别，只是用公式时，要注意各列水平数的差别。例8-5试对例7-2试验数据进行方差分析。课本中为简化计算，对表7-5(p144)的试验数据Xi作了线性变换，实际上没有必要。在不对Xi作变换的情况，请同学们自己再做一次方差分析，作为课外作业去完成。且求S时，对二水平因素用通式和简化式分别计算！总的偏差平方和n_nS产v(Xi-X)2=vXi2-CT1 1i1fT=n-1因素的偏差平方和S分两种类型进行计算：1、对于四水平因素1m2Sj=%Kij-CTrij(对二水平因素，也可用这一通

15、式计算，建议全部用通式计算，以免产生混乱！)m=4,r=n/m=8/4=2fj=m-12、对于二水平因素，简化计算公式为：Sj=1(Kij-K2j)2,n=8nfj=m-1,m=2(方差分析和显著性检验，见书上p181)讨论：(1)方差分析法与极差分析法得出的各因素主次顺序相同，都是A、C、B;(2)由方差分析可知，因素A显著，因素C不显著，而因素B对试验结果无影响(即将Sb并入S中，及S9=Se+SB,fe=fe+fb)；(3)主要因素A的优水平A;不显著因素C,可根据具体情况确定其水平，为缩短加工时间，可选G水平，但从指标值看，还的选G为好；对试验结果无影响的因素B,选B或B均可，从试验的

16、指标可知，A3BG为最佳工艺条件，(即5号试验)。因此，此时指标值最大。(251cm/100g)(极差分析结果：ABG或ABO方差分析结果：ABG或ARG)8.5重复试验和重复取样的方差分析(因时间有限，不讲解！)在实际工作中，用正交表安排试验时，为了提高试验及其统计分析的精确性和可靠性，往往采取重复试验和重复取样，在安排试验时，将同号试验重复做若干次，从而得到在同一条件下若干次试验的数据，叫做重复试验，若在一个试验中，同时抽取若干个样品进行测试，则叫做重复取样。8.5.1 重复试验的方差分析在用正交表安排试验时，若表上各列已被因素及交互作用占满，没有空列，也无经验误差。这时，为了估计试验误差

17、，一般选用更大的正交表以外，还可以重复试验，由于正交本身的需要，有时虽然正交表的所有列并未被因素及交互作用占满，但也要做重复试验。重复试验的方差分析与无重复试验的方差分析比较，有以下几点不同：(1)假设每号试验重复数为S,在计算K%KjKj时，是以各号试验下“S个试验指标数据之和”进行计算；(2)重复试验时，总偏差平方和与及其自由度fT按下式计算：ns2Sr=.”xit2-小4nsfT=ns-1式中：n-试验条件数，即正交表的总试验号;s-各号试验重复数Xit-第i号试验第t次重复试验数据(i=1,2,n;t=1,2,s);T-所有试验数据之和(包括重复试验)；nsT=二,二X,imp(3)重

18、复试验时，各列偏差平方和(S)计算公式中的“水平重复数”改为“水平重复数乘以试验重复数”，修正项CT也有变化，S的自由度fj仍为水平数减1。m2Si=-zKij2-CT,CT=,r=-rsinsmfj=m-1(4)重复试验时，总误差平方和包括空列误差S1和重复试验误差Se2,即Se=$1+S2其总的自由度fe等于$1的自由度fe1与Se2的自由度fe2之和，即：fe=fe1+fe2Se及fe2的计算公式如下nsnsSe2=:二xit2xit)2i=1t=1si1t=1fe2=n(S-1)重复试验时，用Ve=Se/fe检验各因素及其交互作用的显著性。当正交表的各列都已排满因素及交互作用而无空列时

19、(即Se=0和fe1=0)用Ve=Se/fe2来检验因素及交互作用的显著性。例8-6(p183)四因素四水平正交试验，每号试验重复三次，由附表7可知，对四因素四水平试验，选Li6(45)正交表最合适，本例不考虑因素间的交互作用，因素水平如表8-17所示，而表8-18为试验方案与试验结果计算表。一、计算(简略)1 .计算各列水平K值(Kj,Kj,Kj,Kj)和Kj如Ku=6+12.5+17.5+19.2=55.2Kj=55.22=3047.04心=19.2+19.5+18.9+19.2=76.8K45=76.82=5898.242计算各列偏差平方和(S)及其自由度(fj)m2Sj=zKij2-C

20、T,CT=L,r=-rsynsm=1?23032_)2=Kj-=Kij-1912.6943y16312y如Sa=S=：1(K21+K21+K312+K412)-1912.691=-X(3047.04+6528.64+7656.25+6320.25)-1912.6=49.9912同理可得Sb=S=33.42Sc=S3=29.01Sd=S=13.54Se1=S=9.65nsnsSe=mx/，xit)2i3t=1si1111631163=.;xit2'C、Xit)i二tJSi111=(22+22+6,52+6.92)-1/3乂(62+12.52+20,42)=2050.32-2048.31=

21、2.01所以Se=Se+Se=9.65+2.01=11.66fj=m-1=4-1=3fA=fB=fC=fD=3fe1=f5=4-1=3fe2=n(S-1)=16x(3-1)=32fe=fe1+fe2=3+32=35验算:Srns2St，'Xi2-3032163iatans163、Xi2imtt=2050.32-1912.69=137.625另外St=-SjSe2=SaSbScSdS.Se2j1=49.99+33.42+29.01+13.54+9.65+2.01=137.62fTfT=nS-1=163-1=475另外fT=-fj-fe2=fAfBfcfDfe1fe2j1=3532=473

22、.计算方差V尸?fjV-Sa49.99Va-=16.66,fA3_Sc29.01QA7Vc9.67,fc3Ve-Se=11.66=0.33fe35Vb-邑=3342=11ifB3Vd-包=1354=4.51,fD3二、显著性检验1.计算F值VjVeFaVaVe16.66=50.480.33FbVbVe11.41Q33.760.33FcVe0.33坦=29.30FdVdVe451=13.670.332 .查Fa值F：.=(f大，fe)=F：.(3,35)当a=0.05时,F0.05(3,35)=2.88(由附表5内查得到a=0.01时,F0.01(3,35)=4.41(由附表5内查得到3 .显著

23、性检验因为，Fa、Fb、Fc、Fd均大于Fo.01(3,35)D四个因素均高度显著，方差分析表如表8-19所示(P186)三、确定最优条件四个因素的作用均高度显著，且由F值大小可知因素作用的主次顺序为A、B、CD通过比较K值，可知各因素的优水平为AB、G、口，故最优水平组合为ABGQ,表8-18的试验方案中无该水平组合的试验，所以应在最优水平组合下，再安排实施一次试验，并且其试验指标值应大于表8-18中的最大指标值。8.5.2重复取样的方差分析由于重复试验使试验次数成倍增加而增加试验费用，故在实际工作中，更常用的是采用重复取样方法来提高试验的可靠性，重复取样与重复试验在误差偏差平方和的计算上完

24、全一样，但重复取样的误差，反映的是原材料和产品的不均匀性与试样的测量误差，即局部（试验）误差；而重复试验的误差，反映的是整个试验过程中的各种干扰引起的误差，即整体误差。通常，局部误差比整体误差要小，原则上不能用来检验各因素水平间是否存在差异，否则，会得到几乎全部因素及交互作用都是显著的不正确结论。但是，若符合下面两种情况，则可以把重复取样得到的局部误差Se2当作试验误差Se,进行统计检验。（1）正交表中各列已排满，无空列提供一次误差（Sei）,这时，为了少做试验而用重复取样误差（Se2）作为试验误差（Se）,检验各因素交互作用的显著性，若检验结果有一半左右的因素及交互作用不显著，就可以认为这种

25、检验是合理的；（2）若重复取样得到的局部（试验）误差（Se2）与整体（试验）误差（$1）相差不大，也就是说，要求两类误差的F值：匚SeiF二对于给定的信度，有F<Fa（fel,fe2）,说明Sei与Se2的误差不显著，这时，就可以将Se2和Se合并作为试验误差，即Se=Se+Se;Fe=fei+fe2但是，若F>E（fe1,fez）,则两类误差有显著差异，不能合并使用。例8-7三因素三水平正交试验，不考虑交互作用，因此，选用L9（34）正交表最合适。因素水平表见表8-20,试验方案见表8-21（seep188）。重复取1三次，即s=3.解：一、计算1 .计算各列水平的K值（Kj,K

26、j和Kj）和Kj2如:K11=0.655+0.657+0.787=2.099,K21=2.0992=4.406K13=0.760+1.305+0.657=2.722,K23=2.7222=7.40913K和K的计算结果，列于表8-21中.2 .计算各列偏差平方和（S）及自由度（fj）检验：Sj“kj工、飞一迎。k12.866jsiJns33iJ939Sa=Si=1/9X(7.684+15.413+4.406)-2.866=0.190同理可得：Sb=S2=0.00889,Sc=S3=0.0188Sel=S4=0.00622nsnsSe2=Xit2-1/S-(%期)2n=9,s=3itJit493

27、9s='、'、xit2-1/3%xit)2i4tJitJ=(0.2782+0.2512+0.2592)-1/3*(0.7602+1.1622+_2+0.787)=3.100-3.089=0.0110fj=m-1=3-1=2fa=fB=fc=2fel=f4=3-1=2fe2=n(s-1)=9X(3-1)=18验算：Stnsn&=一Xit2-T/ns=3.100-2.866=0.234i3tT4ST八SjSe2=3+Sb+Sc+Sei+Se2=0.234jiQfTfT=ns-1=9X3-1=264fT="ff=fA+fB+fc+fel+fe2=26j£j

28、e2所以计算无误3 .比较两类误差F=(Sei/fei)/(Se2/fe2)=(0.00622/2)/(0.0110/18)=5.09因为F0.01(2,18)=6.01>5.09,所以两类误差可以合并使用。Se=Se+Se=0.00622+0.0110=0.01722Fe=fe1+fe2=2+18=204 .计算方差V=S/fiV=S/fA=0.190/2=0.095同理,Vb=0.00445,Vc=0.0094,Ve=0.000861二、显著性检验Fj=V/Ve1 .求FFa=VA/Ve=0.095/0.000861=110.34FB=M/Ve=0.00445/0.000861=5.168Fc=VC/Ve=0.0094/0.000861=10.922 .查F%F%=(f因,fe尸F%(2,20)%=0.01时,F。.。1(2,20)=5.850c=0.05时,F0.05(2,20)=3.49因为Fa>F0,01所以因素B高度显著；又因为Fo.05<Fb<Fo.01，所以B显著；又因为Fc>F0.01所以因素C高度显著方差分析表见表8-22(p.190)三、最优工艺条件的确定由表8-22可见，因素A,C高度显著，因素B显著，根据Fj值的大小可知，因素的主次顺序为ACB通过比较Kj值可