组内观察值数目不等的单向分类资料.docVIP

二、组内观察值数目不等的单向分类资料 如果资料中含有k组数据，其中第i组含ni个观察值，全部共有个观察值。此类资料观察值的数学模型为： （i＝1，2，…，k，j ＝1，2，…，ni） 因此总变异分解为两大部分：组间变异和组内变异。方差分析表如表7.19所示。 表7.19 组内观察值数目不等的单向分类资料的方差分析表变异 来源自由度平方和均方F值期望均方(EMS)固定模型随机模型组间dft＝k－1MStF＝MSt/ MSe误差dfe＝－kMSe总变异dfT＝－1SST＝其中为各组样本含量的调和平均数，即。但是如果利用随机模型对进行估计，即)时，可能会出现方差估计值为负的奇怪现象。这时可以使用代替进行计算。关于在样本容量不等时，利用或代替n进行计算的讨论。 各组内观察值数目相等的试验资料称为平衡资料；各组内观察值数目不等的试验资料称为不平衡资料。利用平衡资料所得到的参数估计值是无偏的，利用不平衡资料所得到的参数估计值不是无偏的。一般地说，如果是经过精心设计来安排试验的话，应尽可能安排平衡试验，避免采用不平衡试验设计。如果是进行一些野外调查，由于特殊条件的要求或限制，不得不采用不平衡设计时，仍应特别注意估计结果的有偏性。 例7.8 调查四个水库中的氯离子浓度，按水库的水面面积大小取不同数目的样本。数据如表7.20所示，如果各水库的效应为固定模型，试测验各水库间氯离子的浓度之间是否有显著差异。 ⑴ 数据处理： 表7.20 四个水库氯离子浓度的数据资料及其整理水库观 察 值ni Ⅰ121314151416147981382137214 Ⅱ9108891097635715679 Ⅲ9210111213121188088480010 Ⅳ14101112131267287486412　总和2831337113603⑵ 平方和及自由度的分解： 总自由度dfT = 观察值总数 – 1 = 28 – 1 = 27 处理间自由度dft = 处理数 – 1 = 4 – 1 = 3 误差自由度dfe = dfT –dft = 27 – 3 = 24 矫正数C.T. = 总和平方/观察值总数 = 3132/28 = 3498.8929 总平方和 = 3711－3498.8929 = 212.1071 处理间平方和 ＝3603 － 3498.8929 ＝104.1071 误差平方和 = 212.1071 － 104.1071 = 108 * 建议记住文字表述 ⑶ 列方差分析表： 表7.21 例7.8的方差分析表变异来源自由度平方和均方FF0.05F0.01处理间3104.107134.70247.712**3.014.72误差24108.00004.5000　总变异27212.1071　　⑷ 多重比较（以Duncan法为例）： 1. 先计算样本容量的调和平均数或： 2. 计算标准误：0.806 3. 查表并计算各种LSR值（误差自由度dfe = 24）如表7.22所示： 表7.22查得的SSR值及各种LSR值的计算gSSR0.05SSR0.01LSR0.05LSR0.0122.923.962.3533.19233.074.112.4743.33743.154.242.5393.4174. 列梯形表进行比较： 表7.23 例7.8的多重比较梯形表处理平均数Ⅰ145**4**2Ⅳ123*2Ⅲ101Ⅱ9⑸ 作统计推断： 1. 第Ⅰ号水库与第Ⅱ、Ⅲ号水库的氯离子浓度有极显著差异； 2. 第Ⅱ号水库与第Ⅳ号水库的氯离子浓度有极显著差异； 3. 其余各水库的氯离子浓度之间没有显著差异。