北京林业大学《数理统计A》ch5方差分析.docVIP

第五章 方差分析 §5.2单因素方差分析 引例：灯泡厂用四种不同材料的灯丝制成四批灯泡。现在，从四批灯泡中分别随机抽取几只进行使用寿命试验。结果如下： 材料 灯泡寿命（小时） 平均寿命 1600，1610，1650，1680，1700， 1720，1800 1680 1580，1640，1640，1700，1750 1662 1460，1550，1600，1620，1640， 1660，1740，1820 1636.25 1510，1520，1530，1570，1600，1680 1568.33 分别表示四种灯丝材料制成的灯泡寿命。记，. 需检验统计假设 . 灯丝材料：因素， 四种灯丝材料：因素的四个不同水平。 一．问题的一般化： 独立，，. 其中未知。 有样本；样本均值，样本方差； 有样本样本均值，样本方差； … … 有样本；样本均值，样本方差. 问题：利用联合样本 检验统计假设 ， 中至少有两个不相等。 二．平方和的分解 总的离差平方和：， 其中，。 的取值反映联合样本取值的分散程度。 联合样本取值分散的原因： （1）不全相同； （2）可能完全相同，但由于是随机变量，其样本观察值也会大小不一。 为区分这两种情况，将做一个分解： 记 ，， 则 其中，. 记 ， 则 分析： 其中第一项是（第一个水平）的样本离差平方和，反映的样本分散程度；第二项是（第二个水平）的样本离差平方和，反映的样本分散程度；……这些分散完全是随机性造成的。所以，从分解出来的完全是随机性造成的，称为组内离差平方和（误差平方和）。 这样，如果不全相同，那么由此带来的联合样本的分散会在这部分反映出来。称为组间离差平方和（因素的效应平方和）。 三．假设检验问题的拒绝域 如果和比较起来不太大，只能认为也是随机性造成的，即认为是随机性造成的，这时应接受；反之，如果和比较，太大了，就有理由怀疑不完全是随机性造成的，各总体均值不全相同，从而拒绝。 取检验的拒绝域的形式：. 当真时，，独立、同分布，，， 于是，独立、同分布， ，. . 又由，独立，知 ，. 且这个变量独立。所以， ， 即 . 可以证明，当真时， ， 即 ，且与独立。 这样，当真时，记，，则 . 给定显著水平，查表12(P262)得， . 检验的拒绝域：。 检验方法：若，则拒绝； 若，则接受. 四．方差分析表 1．平方和、的计算 方法一：， 其中， 是联合样本的样本方差； ， 其中是第个水平样本的样本方差，于是， . 此法适合数据量较小的情形。 方法二：， 记， 则 （*） 即 （**） 其中，（）， 所以， 输入数据，给出下表： 水平 1 2 求和 按（*）（**）计算平方和： ；。 方法三：使用统计软件Excel、Matlab2．P139）： 变差来源 平方和 自由度 平均平方和 F值 显著性 组间 **或* 或 - 组内 总和 标**，标*，标-。 注：时，检验 ， 将 代入 得 例1．四种不同材料的灯丝制成的灯泡寿命是独立随机变量，，。从四批灯泡中分别随机抽取几只做寿命试验。结果如下： 材料 灯泡寿命（小时） 平均寿命 1600，1610，1650，1680，1700， 1720，1800 1680 1580，1640，1640，1700，1750 1662 1460，1550，1600，1620，1640， 1660，1740，1820 1636.25 1510，1520，1530，1570，1600，1680 1568.33 . 中至少两个不同。 利用方法二：由联合样本观察值得 材料 7 =11760 5 =8310 8 =130905 6 =94107 求和 42570 =2 = 于是， ， . =44360.7. 若利用方法一： =195711.5， =28600+16880+85187.5+20683.3=151350.8 =44360.7. 方差分析表： 变差来源 平方和 自由度 平均平方和 值 显著性 组间 =44360.7 =3 =14786.9 =2.15 - 组内 =151350.8 =22 =6879.6 总和 =195711.5 =25 ， ， 接受，认为四种材料的灯丝制成的灯泡寿命无显著差异。 方法三：使用Excel:工具→数据分析→单因素方差分析 差异源 SS df MS F P-value F crit