第五章方差分析.ppt

第五章 方 差 分 析 5.1 方差分析中的有关概念 5.2 单因素方差分析 5.3 双因素方差分析 5.4 均值估计与多重比较 5.1 方差分析中的有关概念 5.1.1 单因素方差分析问题与模型 5.1.2 双因素方差分析问题与模型 5.1.3 方差分析中的基本假定 5.1.1 单因素方差分析问题与模型 1. 数学模型 进行单因素方差分析时，需要得到如图5-1所示的数据结构。 设xij表示第i个总体的第j个观测值（ i = 1，2，…，m， j = 1，2，…，ni），希望由此对不同水平下总体的均值进行比较。 对此，观察到的xij常用以下的模型表示： xij = ?i + ?ij ，1≤j≤ni，1≤i≤m 其中?i表示第i个总体的均值，?ij为随机误差，在方差分析中为了得到有效的检验法还常假定?ij满足： ● ?ij为相互独立的； ● ?ij都服从正态分布，且?ij的均值都为0，方差都相同。 2. 方差分析的过程 为了方便起见，可将?i记为： ?i = ? + ?i 其中 称为总均值，?i = ?i – ?，i = 1，2，…，m称为因素A的第i个水平的附加效应，这样比较不同水平下均值是否相同。问题的检验假设： H0：?1 = ?2 = … = ?m，H1：?1，?2，…，?m不全相等； 就可以表示为： H0：?1 = ?2 = … = ?m= 0，H1：?1，?2，…，?m不全为零。 在H0成立下检验用统计量： 其中 、 称为组间、组内（变差）平方和；这里 称为组内平均； 称为总平均，n = n1 + n2 + … + nm；另外 称为全部（变差）平方和；可以证明 SST = SSMA + SSE。 当原假设成立时，各总体均值相等，各样本均值间的差异应该较小，模型平方和也应较小，F统计量取很大值应该是稀有的情形。 所以对给定显著性水平α?(0, 1)，若p = P{F ? F0} α，则拒绝原假设H0（F0为F统计量的观测值），可以认为所考虑的因素对响应变量有显著影响；否则不能拒绝H0，认为所考虑的因素对响应变量无显著影响。 3. 方差分析表 通常将上述计算结果表示为表5-1所示的方差分析表。 表5-1 单因素方差分析表 其中，MSA = SSMA/(m – 1)，MSE = SSE/(n – m)。利用方差分析表中的信息，就可以对因素各水平间的差异是否显著做出判断。 例. 设有三台机器，用来生产规格相同的铝合金薄板. 取样，测量薄板的厚度精确至千分之一厘米. 得结果如表所示. 问不同机器对生产的铝合金板的厚度有无影响? Solution→Analysis→Analyst(分析员系统) (出现空白数据表)→ File→Open By Sas Name…　 (在Make one selection窗口中)work 选中数据名E411→ (OK) →Statistics→ANOVA → One-Way ANOVA… Independent填分类变量c→ Dependent因变量y Plots—可选择分水平的盒形图(Box-＆-Whisker Plot) ,条形图(Bar Chart)及均值、标准差图 Means—Comparisons Methods给出了10种多重比较的方法 α为选择的显著性水平, Breakdown可按水平分组出描述性统计量 → OK (点击运行后的结果树标签则会打开相应图. 5.1.2 双因素方差分析问题与模型 1. 无交互作用的双因素方差分析 对于多因素问题，通常考虑有重复观测的情形，其数据结构如图5-2所示。 图5-2 双因素方差分析中数据结构 若第一个因素A有l个水平，第二个因素B有m个水平。在因素A的第i个水平和因素B的第j个水平下进行了多次观测，记为{xijk，1≤k≤n}。 对xijk考虑以下模型： xijk= ? + ?i + ?j + ?ijk， 1≤i≤l，1≤j≤m，1≤k≤n 其中?表示平均的效应，?i和?j分别表示因素A的第i个水平和因素B的第j个水平的附加效应，