第十章 方差分与正交试验设计.docVIP

第十章 方差分析与正交试验设计 方差分析与试验设计是英国统计学家和遗传学家费希尔进行农业试验发展起来的通过试验获取数据并进行分析的统计方法。方差分析讨论的是生产和科学试验中有哪些因素对试验结果有显著作用，哪些因素没有显著作用。讨论的是一个因素对试验结果是否有影响称为一元方差分析，讨论的是多个因素对试验结果是否有影响称为多元方差分析．对于因素多于两个的方差分析，公式变得相当复杂，试验次数较多，我们介绍一个试验次数少的试验设计方案，正交试验设计。 10．1 一元方差分析 人们常常通过试验来考察了解各种因素对产品或成品的性能，成本、产量等的影响，我们把性能、成本、产量等统称为试验指标。有些指标可以直接用数量表示，称为定量指标；不能直接用数量表示的，称为定性指标，可按评定结果打出分数或评出等级，这时就能用数量表示了。在试验中，影响试验指标的原因称为因素。因素在试验中所处的各种状态称为因素的水平，某个因素在试验中需要考察它的几种状态，就称它为几水平的因素。 在生产实践和科学试验中，人们经常要研究这样的问题：如果改变生产条件是否会对产 品（指标）产生显著影响?如果改变试验条件是否会对试验结果（指标）产生显著影响?方差分析的作用就在于通过对试验数据的统计分析，从而推断试验数据间的差异是由于生产条件的改变还是由于随机误差的影响，并分析出最佳的试验条件。为此弄清楚方差分析处理问题的基本思想，下面举例说明。 例10.1.1 某灯泡厂用四种不同配料方案制成的灯丝生产四批灯泡,在每批灯泡中取若干个做寿命试验，它们的寿命分别记为xij，其中下标i表示第i批灯泡，第二个下标j表示第j次试验。具体数据如下 表10.1.1 四批灯泡的寿命试验表 品种 寿命（小时） A1 A2 A3 A4 1600，1610，1650，1680，1700，1720，1800 1580，1640，1640，1700，1750 1460，1550，1600，1620，1660，1740，1820,1640 1510，1520，1530，1570，1600，1680 现在要研究的问题是灯丝的不同配料方案，即不同的品种对灯泡寿命有无显著影响。 在这里灯泡的寿命就是指标，灯泡品种就是因子，四种不同品种的灯泡就是四个水平,因此这是一个单因子四水平试验。我们将每一种配料制成的灯泡，其寿命看成同一总体，而不同品种的灯泡就是不同总体，因而出现四个不同总体。每一种的灯泡寿命都有一个理论上的平均值，即分布的数学期望，不同品种的灯泡的寿命的数学期望可能有显著差异，也可能没有显著差异，试验的目的就是通过假设检验对这个问题给出一个推断。一般可假定母体的方差相同。由于其他试验条件相同，如果灯泡品种对灯泡寿命无显著性影响，我们可认为四个总体的概率分布相同，换句话说，灯泡品种对灯泡寿命是否有显著性影响，就是要检验四个总体的均值是否相等．按参数估计的假设检验方法可以逐个地进行检验，但这个方法显得繁而复杂．特别当水平数较多时，需要做许多假设和检验，计算量也相当大．如果能导出一个可以用来检验所有这些假设的统计量，那么解决这样的问题就方便多了。方差分析就是解决这样的问题． 假设试验只考虑一个因素A，它有I个水平A1，A2，…，AI，总共有次试验，xij表示第i水平第j次试验，其数据如下表 表10.1.2 一元方差试验数据表 水平 试验结果 我们再作如下假设：为个子总体，且相互独立，，而为的样本。显然个水平对试验结果有无显著性影响，就是看是否为相同的总体，或它们的分布是否相同。由于它们都是正态总体，就只要看它们分布的参数是否相同，已知方差相同，这就只须判断数学期望是否相等。换句话说，只要在一定的显著性水平上检验统计假设 令 ， 分别表示第个子总体的样本均值（组平均值）和总体样本均值（总平均值）。总偏差平方和 它描述全部数据离散程度（总波动）的大小。容易证明 10.1.1 其中 ， 反映的是各子总体样本均值（组平均值）的不同而引起的误差，是各组平均值与总体样本平均值的离差平方和，它表示因试验水平差异带来的误差大小，称为组间偏差平方和，也称为系统误差。反映的是每一个子总体的（组内）数据不同而引起的误差，是每个观测值与其组内平均值的离差平方和，它表示试验误差的大小，称为组内偏差平方和，也称为误差平方和。因此，通过的大小可以反映原假设是否成立。若显著地大于，说明各子总体（水平）之间差异显著，那么可能不成立。这种比较方差大小来判断原假设是否成立的方法就是方差分析的由来。那么的值大到什么程度可以否定呢？在理论上已经证明 10.1.2 统计量可以作为判断是否成立的检验统计量。在给定显著水平的情况下，当时，则拒绝，认为因素对试验的指标是显著的，否则接受。在实际进行一元方差分析时，通常将有关的统计量连同分析结果列在一张