方差分析课件.ppt

方差分析单因素试验的方差分析双因素试验的方差分析一、问题【例】某饮料生产企业研制出一种新型饮料。饮料的颜色共有四种，分别为橘黄色、粉色、绿色和无色透明。这四种饮料的营养含量、味道、价格、包装等可能影响销售量的因素全部相同。现从地理位置相似、经营规模相仿的五家超级市场上收集了前一时期该饮料的销售情况，见表1。试分析饮料的颜色对销售量是否产生影响。单因素一元方差分析的MATLAB实现 anova1函数 功能：单因素一元方差分析 调用方式： p=anova1(X) p=anova1(X,group) p=anova1(X,group,displayopt) [p,table]=anova1(…) [p,table,stats]=anova1(…)单因素试验的方差分析在MATLAB工具箱中提供了单因素方差分析函数，格式如下：P＝anova1(X)是对数据矩阵X（m×n矩阵）中各列进行单因素n水平的方差分析.n个水平有无显著差异．函数返回原假设的概率值P,概率值小于0.05认为差异是显著的，小于0.01则认为差异是极显著的．anova1函数输出:两个图表第一个为标准的方差分析表第二个为X数据的盒形图X数据的盒形图如果盒形图的中心线差别很大，则对应的F值越大，相应的概率值越小．方差分析表中有六列数据：第一列显示数据来源；第二列显示三种离差平方和（SS）；第三列显示三种离差平方和的自由度（df）；第四列显示三种离差平方和的均值平方和，即SS/df；第五列显示假设检验的F统计量，即均方值的比值；第六列显示概率值P，它由F的累积分布函数得．F越大，概率值P越小．例1有三种型号的电池A，B，C分别是三家工厂生产的．现各抽取5只评比质量，经过试验的寿命（小时）如下：A4048384245B2634302832C3940435050试问：这三种型号电池的寿命有无显著差异？程序及结果如下：X=[4048384245;2634302832;3940435050]XP=anova1(X)X=404838424526343028323940435050ans=402639483440383043422850453250P=3.0960e-004P＝anova1(X,group)该函数与P＝anova1(X)的不同是用group中的值对矩阵X各样本的盒形图进行标识．group的每一行包含的标识对应矩阵X的相应列，其长度应与矩阵X的列相等．优点是：anova1的向量输入形式不要求每个样本都具有相同数量的观察值，因此它对于不平衡数据也是适用的．例2对例1采用函数P＝anoval（X，group）进行方差分析．程序如下：x=[404838424526343028323940435050];group={A,A,A,A,A,B,B,B,B,B,C,C,C,C,C};p=anova1(x,group)输出结果和例1中输出结果不同处是：用A，B，C标出了box图中的总体.[p,table]=anova1(x,group)

以单元数组的形式返回方差分析表（包括列和行的标识）．列和行的有标识的方差分析表p=3.0960e-004table=Columns1through5SourceSSdfMSFGroups[615.6000][2][307.8000][17.0684]Error[216.4000][12][18.0333][]Total[832][14][][]Column6ProbF[3.0960e-004][][][P,tabla,stats]=anov