01方差分析.docx
方差分析专题单因素试验的方差分析（一）单因素试验在科学试验和生产实践中，影响一事物的因素往往是许多的。例如，在化工生产中，有原料成分、原料剂量、催化剂、反应温度、压力、溶液浓度、反应时间、机器设备及操作人员的水同等因素。每一因素的变更都有可能影响产品的数量和质量。有些因素影响较大，有些较小。为了使生产过程得以稳定，保证优质、高产，就有必要找出对产品质量有显著影响的那些因素。为此，我们需进行试验。方差分析就是依据试验的结果进行分析，鉴别各个有关因素对试验结果影响的有效方法。在试验中，我们将要考察的指标称为试验指标.影响试验指标的条件称为因素因素可分为两类，一类是人们可以限制的（可控因素）；一类是人们不能限制的。例如，反应温度、原料剂量、溶液浓度等是可以限制的，而测量误差、气象条件等一般是难以限制的。以下我们所说的因素都是指可控因素。因素所处的状态，称为该因素的水平（见下述各例）。假如在一项试验中只有一个因素在变更称为单因素试验,假如多于一个因素在变更称为多因素试验例1设有三台机器，用来生产规格相同的铝合金薄板。取样，测量薄板的厚度精确至千分之一厘米。得结果如表9.1所示。表9.1铝合金板的厚度机器I机器11机器In0.2360.2570.2580.2380.2530.2640.2480.2550.2590.2450.2540.2670.2430.2610.262这里，试验的指标是薄板的厚度。机器为因素，不同的三台机器就是这个因素的三个不同的水平。我们假定除机器这一因素外，材料的规格、操作人员的水同等其它条件都相同。这是单因素试验。试验的目的是为了考察各台机器所生产的薄板的厚度有无显著的差异。即考察机器这一因素对厚度有无显著的影响。例2下面列出了随机选取的、用于计算器的四种类型的电路的响应时间（以毫秒计）。表9.2电路的响应时间类型I类型11类型In类型IV192016182221152220331819182726154017这里，试验的指标是电路的响应时间。电路类型为因素，这一因素有4个水平。这是一个单因素试验。试验的目的是为了考察各种类型电路的响应时间有无显著差异。即考察电路类型这一因素对响应时间有无显著的影响。例3一火箭运用了四种燃料，三种推动器作射程试验。每种燃料与每种推动器的组合各放射火箭两次，得结果如下（射程以海里计）。表9.3火箭的射程推动器（B）B?458.256.265.352.641.260.8A249.154.151.6燃料（八）42.850.548.4As60.158.370.973.239.240.775.858.248.771.55141.4这里，试验的指标是射程，推动器和燃料是因素，它们分别有3个、4个水平。这是一个双因素的试验。试验的目的在于考察在各种因素的各个水平下射程有无显著的差异，即考察推动器和燃料这两个因素对射程是否有显著的差异。本节限于探讨单因素试验，我们就例1来探讨。在例1中，我们在因素的每一水平下进行了独立试验，其结果是一个随机变量。表中数据可看成来自三个不同总体（每个水平对应一个总体）的样本值。将各个总体的均值依次记为从，?,3。按题意须要检验假设%：M=42=M3HI:4“G外不全相等现在进而假设各总体均为正态变量，且各总体的方差相等，那么这是一个检验同方差的多个正态总体均值是否相等的问题。下面所要探讨的方差分析法，就是解决这类问题的一种统计方法。现在起先探讨单因素试验的方差分析。设因素有S个水平A，A2,4,在水平Aj（J=I,2,s）下，进行%（nj2）次独立试验，得到如下表的结果。我们假定：各个水平4（=1,2,s）下的样原来F,迎”，再“自具有相同方差2,均值分别为j（7=1,2,5）的正态总体N（j,/）,勺与/未知。且设不同水平A.下的样本之间相互独立。由于马N（勺,O?）,即有为/N（0,c2）,故勺-勺可看成是随机误差。记xij-j-ij»则勺可写成xij=+%,i=12,%;j=l,2,马N（0,/）,各与独立，,（1.1）其中勺与人均为未知参数。（1.1）式称为单因素试验方差分析的数学模型。这是本节的探讨对象。方差分析的任务是对于模型（1.1）,I0检验S个总体N（MQ2）,n（2q2）,N（m,q2）的均值是否相等，即检验假设H0x=2=sH1:必,2,，4不全相等。（1.2）2°作出未知参数必,2,，$，人的估计。为了将问题（1.2）写成便于探讨的形式，我们将外,2,，4的加权平均值勺勺ny=l记为，即"=叫勺（13）n7=1其中二£叼。4称为总平均。再引入；=1M=j-R,j=,2、,s（1.4）此时有用+%4+凡。.=0，%表示水平4下的总体平均值与总平均的差异，习惯上将称为水平4的效应。利用这些记号，模型（1.1）可改写成Xij=+S,÷ij,*ErljBj=O,/=1,2,wy;J=1,2,5,;=|JjN(0,2),各%独立而假设(1.2)等价于假设(1.2)'“0:心=&=瓦=OHi:"0,瓦不全为零。这是因为当且仅当必=2=4时J=，即与=0,(/=1,2,s)。(二)平方和的分解下面我们从平方和的分解着手，导出假设检验(1.2)'的检验统计量。引入总平方和ST=(x2C/=1*=1(1.6)其中途色nZ=I=是数据的总平均。Sr能反映全部试验数据之间的差异，因此S7又称为总变差。又记水平Aj下的样本平均值为工八即我们将»写成s=(-jv)2=k-X)+(-X)2y=l=l7=1/=I5zry_Xn,_.Vnj_=(-)2+U-)2+2X(Xy-Xj×Xj-X)J=I/=I；=1<=1J=I<=留意到上式第三项(即交叉项)SJj、-Xj)(xj-x)=2(x-x)E(Xij-xj)=07=1/=I=l1.i=I于是我们就将s7分解成为S=Se+S,(1.8)其中SE=(-x.j)21(1.9)=1Z=I3_'_'、-22SA=ZZ(X1.X)2=Z%(Xj7)2=Z%KjiJ(1.10)jli三Ijj上述SE的各项(-1j)2表示在水平A,下，样本视察值与样本均值的差异，这是由随机误差所引起的。SE叫做误差平方和。SA的各项(TJ-幻2表示为水平下的样本平均值与数据总平均的差异，这是由水平Aj引起的。SA叫做因素A的效应平方和。(1.8)式就是我们所须要的平方和分解式。(三)Se，S八的统计特性为了引出(1.2)'的检验统计量，我们依次来探讨SA的一些统计特性。(1) SE的统计特性将SE写成SE=X(_彳1)2+(Xi2_工，2)2+(xis-Xs')2(1.11)J=I=li=ltli留意到£(七一工/)2是总体N(j,b2)的样本方差的叫一倍，于是有/=1丐_£(%f)2-/51)J因各马独立，故(1.lI)式中各平方和独立。由力?分布的可加性知,%2恪(厂1)即-f("s),(1.12)由(1.12)式还可知，SE的自由度为一S。且有E(Sf)=(n-s)2(1.13)(2) SA的统计特性我们看到SA=ZZ(X./-x)2=nj(.j-X)2是S个变量"J(xj-1)；=1i=l7=1(/=1,2,/)的平方和，它们之间仅有一个线性约束条件名MIM(X1.X)=SJ=I7=1nj(x.j-X)=O故知SA的自由度为5-1。再由(1.3),(1.6)及勺的独立性，知X-NwE)n即得E(SA)=En.".j-nx(1.14)E(.j)-nE(x)t>U;)+辰町)fl-11O()+反)2=,z；=1%7=1£1+2+【/i>-n<+-n(厅，=s2+n2+2Elnjbj-2-nj=lJ=由(1.1)'式，知2%=0,故有7=1(1.15)(1.16)E(SA)=(S-I)M+Z叩号；=i进一步还可以证明SA与SE独立，且当为真时SJ-z2(-i)证略。思索:当”。为真时，整个样原来自什么总体？(四)假设检验问题的拒绝域现在我们可以来确定假设检验问题(1.2)'的拒绝域了。由(1.15)式知，当"o为真时=2(1.17)5-1Sv即T是的无偏估计。而当为真时，Z勺b；0,此时S-IJ=I力凤鸟_)=22(1.18)5-15-1又由(1.13)式知=2(1.19)n-s即不管”0是否为真，£都是的无偏估计。n-s综上所述，分式SAS-ISEn-s的分子与分母独立，SE的分布与40无关，分母的数学期望总是当，0为真时，分子的数学期望为2,而当修为真时，由(1.18)式分子的取值有偏大的趋势。故知检验问题(1.2)'的拒绝域具有形式n-s其中上由预先给定的显著性水平。确定。由(1.12),(1.16)式及SE与SA的独立性知，当HO由此得检验问题(1.2)'的拒绝域为SAF=sFa(s-,n-s)(1.20)n-s上述分析的结果可排成表95的形式，称为方差分析表表9.5单因素试验£方差分析表方差来源平方和自由度均方尸比因素ASA5-1SaOA=S-1SaFSE误差SEn-sQSEn-s总和Stn-_._a,表中S2j,S分别称为S.,SE的均方。5-1n-s思索：当”。为真时，均方的数学期望分别是什么？因此均方又可以称什么？另外，由于在S7中个变量/-7之间仅满意一个约束条件(1.6),故57的自由度为H-Io例4如上所述，在例1中须要检验假设"：1=2=3H：»2，3不全相等试取=0.05,完成这一假设检验。解：表9.6例4的方差分析表方差来源平方和自由度均方尸比因素A0.0010533320.0005266732.92误差0.00019200120.00001600总和0.0012453314因”os(2,12)=3.89<32.92,故在水平0.05下拒绝修，认为各台机器生产的薄板厚度有显著的差异。例5设在例2中的四种类型电路的响应时间的总体均为正态，且各总体的方差相同。又设各样本相互独立。试取=0.05
