方差分析-SAS.ppt

《方差分析-SAS.ppt》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、方差分析方差分析主要过程说明单因素方差分析双因素无交互作用方差分析方差分析主要过程说明ANOVA过程格式及使用说明GLM过程格式及使用说明ANOVA过程格式及使用说明对于平衡数据资料（各水平下等重复，数据没有丢失），一般用ANOVA过程.过程格式：PROCANOVA[选项]；CLASS处理因素；MODEL因变量=效应表[/选择项]；MEANS效应表[/选择项]；PROCANOVA语句的选项主要有：DATA=数据集名指明要分析的SAS数据集，缺省时SAS将使用最近建立的数据集.OUTSTAT=输出数据集指定分析计算结果

2、输出的数据集名.CLASS语句指明分类变量，是ANOVA过程的必需语句，并且必须出现在MODEL语句之前.分类变量可以为数值型或字符型，分类变量的个数表示方差分析的因素个数.MODEL语句定义分析所用的效应模型，即方差分析的因变量和效应变量.在方差分析过程中，关键在于定义线性数学模型，常用的模型定义语句有：MODELy=a单因素一元方差分析MODELy=ab双因素无交互作用一元方差分析MODELy=aba*b双因素有交互作用一元方差分析MEANS语句用来计算该语句所列的每个效应所对应的因变量均值，其选项用于设定多重比

3、较的方法.过程说明：GLM过程格式及使用说明GLM即广义线性模型（GeneralLinerModel）过程，对于非平衡数据，应采用GLM过程.它使用最小二乘法对数据拟合广义线性模型.该过程功能强大，可用于多种不同的统计分析中.GLM过程用于方差分析时，主要语句和使用格式与上述ANOVA过程类似.单因素方差分析单因素等重复方差分析(ANOVA过程)单因素不等重复的方差分析（GLM过程）单因素等重复方差分析(ANOVA过程)应用实例源程序运行结果应用实例一个工厂用三种不同的工艺生产某种电池.从三种工艺生产的电池中分别抽取

4、5个样品，测得样品寿命的数据如下（单位小时）：A重复工艺1工艺2工艺314026392463440338304344228485443250平均423044我们要研究的指标是电池的寿命，工艺是影响寿命的一个因素，三种工艺分别是该因素的三个水平.在试验中我们假设其它因素都处于相同的状态.这里我们希望利用上面得到的数据来考察“工艺”的不同是否对“寿命”这个指标有影响？源程序Dataexam;/*建立数据集*/Dotrt=1to3;/*3个处理（trt）分别为1、2、3*/DoI=1to5;/*每个处理下5次重复*/Inp

5、utx@@;Output;End;End;Cards;404638424426343028323940434850;Procanova;/*调用方差分析过程*/Classtrt;/*定义处理为分类变量*/Modelx=trt;/*定义效应模型*/Title'方差分析';Run;运行结果AnalysisofVarianceProcedureDependentVariable:XSourceDFSumofSquaresMeanSquareFValuePr>FModel2573.33333333286.6666666719

6、.770.0002Error12174.0000000014.50000000CT14747.33333333R-SquareC.V.RootMSEXMean0.7671729.8479823.8078865538.66666667以上结果相当于方差分析表，F值为19.77，显著性水平为0.0002，小于0.01，说明各处理间的均值差异极显著.单因素不等重复的方差分析（GLM过程）应用实例源程序运行结果某食品公司对一种食品设计了四种新包装.为了考察哪种包装受欢迎，选取了十个有近似相同销售量的商店作试验，其中两种包装各

7、指定两个商店销售，另两种包装各指定三个商店销售.在试验期间，各商店的货架排放位置，空间都尽量一致，营业员也采用相同的促销方法.一段时间的销售量记录如表所示应用实例包装类型商店1121419242181217303132130395754源程序Dm"LOG;CLEAR;OUTPUT;CLEAR;";Datanew;Inputstr$x@@;/*x表示含销售量，str表示包装处理*/Cards;A112A118A214A212A213A319A317A321A424A430;Procprint;Title"单因素不等重复

8、方差分析";Procglm;Classstr;Modelx=str;Run;运行结果GeneralLinearModelsProcedureDependentVariable:XSourceDFSumofSquaresMeanSquareFValuePr>FModel3258.0000000086.0000000011.220.0071Erro