[理学]第4章 方差分析1.pptVIP

第四章 方差分析目录 方差分析 §4.1 单因素方差分析 §4.2 两因素方差分析 §4.3 正交实验方差分析 §4.4 方差分析条件的检验 §4.5 方差分析在医学案例中的应用 方差分析 在实际工作中有时需要对多个正态总体的均值进行比较，方差分析是解决此类问题的常用方法。 方差分析 1 单因素方差分析 例1 设有三个小麦品种，经试种得第公顷产量的数据如下（单位：） 品种1：4350 4650 4080 4275 品种2：4125 3720 3810 3960 3930 品种3：4695 4245 4620 现问不同品种的小麦产量之间有无显著的差异？ 方差分析 1 单因素方差分析 data var1; input kind$ yield@@; cards; 1 4350 1 4650 1 4080 1 4275 2 4125 2 3720 2 3810 2 3960 2 3930 3 4695 3 4245 3 4620 ; proc anova; class kind; model yield=kind; means kind/snk t alpha=0.05; /*多重比较*/ means kind; run; 方差分析 1 单因素方差分析 方差分析 1 单因素方差分析 多重比较的结果1：LSD比较结果 方差分析 1 单因素方差分析 多重比较的结果2：SNK比较结果 方差分析 1 单因素方差分析 例1 结果 多重比较的结果： 多重比较结果在表示上有如下约定： 标有相同字母的组，表示他们之间没有显著差异； 对标有不同字母的组，则表示有显著差异。 从二者比较的结果看所得的结论一致，也就是品种1和品种3的产量之间没有明显差异，而品种2与品种1和3的产量之间均有明显差异。 方差分析 1 单因素方差分析 means kind;语句输出的结果如下： 方差分析 1 单因素方差分析 PROC ANOVA（GLM）选择项1； CLASS 变量表； MODEL 因变量=自变量表； MEANS 效应/选择项2； BY 变量表； FREQ 变量表； TEST H=效应项 E=误差项； 方差分析 1 单因素方差分析 PROC ANOVA（GLM）的常用选择项为： DATA=SAS数据集 指明ANOVA（GLM）过程要处理的数据集，缺省值为SAS最近产生的数据集。 OUTSTA=SAS数据集 将结果输出到指定的数据集中。 方差分析 1 单因素方差分析 CLASS语句定义分组变量。 MODEL语句指定因变量和自变量，因变量为连续变量，自变量常常是分组变量。设A、B为分类变量，y为因变量，则以下几个MODEL语句的含义如下： ①MODEL y=A；单因素（A因素）方差分析。 ②MODEL y=A B 主效应模型，无交互作用的两因素方差分析。 ③MODEL y=A B A*B 析因模型（带交互作用项），有交互作用的两因素方差分析 其中CLASS语句必须出现在MODEL语句前面。 方差分析 1 单因素方差分析 MEANS 语句要求ANOVA（GLM）过程计算出该语句后列出的每个水平所对应的因变量的均值。选择项2指明进行多重比较的方法，常用的TUKEY、SCHEFFE、LSD、DUNCAN和SNK等，同时还可规定检验的显著性水平的值（由ALPHA=值，缺省值为0.05）。不用选择项2时，可计算各水平对应因变量的均值。 TEST语句要检验的效应和误差项。其中H=后指定的是要检验的效应，E=后指定作为误差项的效应，此项必需指定，缺省是用MSE作为误差项。 方差分析 方差分析 2 两因素方差分析 方差分析 2 两因素方差分析 方差分析 2 两因素方差分析 方差分析 2 两因素方差分析 方差分析 2 两因素方差分析 例2 不考虑交互作用的两因素方差分析 为了考察蒸馏水的pH值和硫酸铜溶液浓度对化验血清中白蛋白与球蛋白的影响。pH值（A）取四种水平其值为5.40 5.60 5.70 5.80，记为A1、A2、A3和A4，硫酸铜浓度(B)取三种水平其值分别为0.04 0.08 0.10，记为B1、B2和B3，采用两因素的全面试验，所得结果如下 方差分析 2 两因素方差分析 例2 方差分析 2 两因素方差分析 方差分析 2 两因素方差分析 例2 Source DF Anova SS F Value Pr F A 3 5 40.95 0.000