Matlab汇总matlab-方差分析

方差分析,概述方差分析是研究一种或多种因素的变化对试验结果的观测值有否显著影响，从而找出较优的试验条件或生产条件的一种常用的数理统计方法。人们在试验中，所考察的数量指标如产量、性能等称为观测值，影响观测值的条件称为因素，因素的不同状态称为水平，一个因素可以有多个水平，在一项试验中，可以得出一系列不同的观测值。引起观测值不同的原因是多方面的，有的是处理方式或者条件不同引起的，称为因素效应（处理效应、条件变异），有的是试验过程中偶然因素干扰或观测误差导致的，称为试验误差。方差分析的主要工作是将观测数据的总变异按照变异原因的不同分解为因素效应和试验误差。并对其做出数量分析，比较各种原因在总变异中所占的重要程度，作为统计推断的主要依据。由此确定进一步的工作方向。,例1利用四种不同的配方的材料A1,A2,A3,A4生产出来的元件，侧得其使用寿命如下表所示。问四种不同的配方下元件使用寿命有无显著差异？,单因素方差分析,分析问题：这个问题中，影响元件寿命的因素是材料配方，四种不同的配方表明因素有4种水平。这样的试验称为单因素四水平试验。从观测数据可以看出，不仅不同配方下元件的寿命不同，且同一配方下元件的寿命也不一样。数据波动的原因归纳如下：其一，同样配方下，试验条件大体相同，因此，数据的波动由于其它随机因素的干扰引起的。设想在同一配方下，元件的寿命应该有一个理论上的均值，而实测寿命数据与均值的偏离即为随机误差，此误差应该服从正态分布。其二，不同配方下，使用寿命有不同的均值，它导致不同组的元件寿命数据不同。,分析结论的推广：设试验只有一个因素A在变化，其它因素都不变。A有r个水平A1,A2,Ar，在水平Ai下进行ni次独立观测，得到试验指标列表如下,x(ij)表示在因素A的第i个水平下的第j次试验的试验结果。,将水平Ai下的试验结果x(i1),x(i2),x(I,ri)看着来自第i个正态总体XiN(mu(i),sigma2),其中，mu(i),sigma2均未知，且X(1),X(2),Xr相互独立。,Mu(i)是第i总体的均值，epsilon(ij)是相应的试验误差，比较A的r个水平的差异归结为比较这r个总体地均值，即检验假设：,统计模型1,假设（1）,Mu表示总和的均值，alpha(i)为水平Ai对指标的效应，则模型1可以等价变形到模型2,统计模型2,那么，假设（1）可以等价为如下假设（2）,假设（2）,为了导出H0的检验统计，构造：,经过计算可以证明如下平方和分解公式,称ST为总的离差平方和（总变差），描述所有观测值的离散程度SE为随机误差的影响，称为误差平方和或组内平方和SA称为因素A的效应平方和或组间平方和。由模型2和假设（2）成立条件下可以得到如下统计量及其分布,F就是H0的检验统计量，对于给定alpha，用Falpha（r-1,n-r)表示F-分布的上分位点。若FFalpha（r-1,n-r),则拒绝H0,即认为因素A的r个水平有显著差异。也可以计算P-值，PFFalpha（r-1,n-r)，显然，P-值小于alpha等价于FFalpha（r-1,n-r)，表示在显著水平alpha下的小概率事件发生了。意味着应该拒绝H0,当P-值大于alpha时，因该接受H0。,Anoval功能：单因素方差分析格式：p=anova1(X)p=anova1(X,group)p=anova1(X,group,displayopt)p,table=anova1(.)p,table,stats=anova1(.),说明：anova1(X)对样本X中的两列或多列数据进行均衡的单因素方差分析，函数返“零假设”成立的概率值。如果概率值接近于零，则零假设值的怀疑，拒绝H0。anova1(X,group)对样本X中由向量group索引的两组或多组数据进行单因素方差分析以比较个列的均值。group向量中值为整数，表示各个组别的观测值的个数。,算例，续例1x=1600,1610,1650,1680,1700,1700,1780,1500,1640,1400,1700,1750,1640,1550,1600,1620,1640,1600,1740,1800,1510,1520,1530,1570,1640,1600;N=11111112222233333333444444;p,table,stat=anova1(x,N,displayoff)p=0.1208,table=SourceSSdfMSFProbFGroups4.9212e+00431.6404e+0042.16590.1208Error1.6662e+005227.5737e+003Total2.1583e+00525,stat=gnames:4x1celln:7586source:anova1means:1.6743e+00315981.6488e+0031.5617e+003df:22s:87.0272,回答问题：根据计算结果，F=2.1659F0.05(3,22)=0.1156,且p=0.12080.05，没有理由拒绝H0,即认为四个水平没有显著差异。,单因素方差分析表,双因素方差分析,例2在一个农业试验中，考虑4种不同的种子品种A1,A2,A3,A4和三种不同的施肥方法B1,B2,B3得到产量数据表如下，试分析种子与施肥量对产量有无显著影响。,分析1：不考虑交互作用,在因素A,B的每种水平组合（A(i),B(j)下独立一次试验的到观测值x(i,j),i=1,2,r;j=1,2,s。在不考虑因素的交互作用时，样本模型为,模型3,其中,判别A,B的影响是否显著等价于检验下列假设,假设（3）,在模型3和假设（3）成立的条件下，构造检验假设（3）的统计量,同单因素方差分解道理，有,当H01成立时,且SE与SA相互对立，则,同理，当H02成立时，,通过计算，同样可以得到双因素的方差分析表,Anova2功能：双因素方差分析格式：p=anova2(X,reps)p=anova2(X,reps,displayopt)p,table=anova2(.)p,table,stats=anova2(.)说明：anova2(X,reps)进行双因素方差分析，比较数据样本X多行和多列数据得均值。不同列中数据表示单一因素变化情况，不同行中数据表示另外一个因素变化情况。如果行-列（单元）由不止一个数值，用reps来表明每个单元的多个观测值的不同标号。nnova2返回零假设成立的概率值p。如果概率接近于零，则零假设值得怀疑。通常认为，p-值小于0.05或0.1，结果较显著。算例续例2,x=325,292,316;317,310,318;310,320,318;330,370,365;p=anova2(x,1)p=0.72910.0413,根据计算结果，A因素的各个水平间有较显著差异,而B因素各个水平间差异不显著。,分析2：考虑交互作用的方差分析,例3研究种树与地理位置对松树生长的影响，对四个地区的三种同龄松树的直径进行测量得到数据如下表所示（cm)。A1,A2,A3表示不同的树种，B1,B2,B3,B4表示不同的地区。对每种组合，进行了5次测量，对此试验结果进行方差分析。,在分析1的基础上，增加sigma(i,j)表示A(i)与B(j)的交互效应。,其它同分析1计算公式,数据分解为下面的模型4,模型4,判别A,B及交互效应的影响的显著性等价于下面的假设(4),假设（4）,在分析1中，稍微修改离差的分解公式，有,为交互效应平方和,当H01成立时，有,当H02成立时，有,当H03成立时，有,有交互效应的方差分析表,算例，续例3X=232521