课程笔记
9.1单因素方差分析
模型与假设
设有 个总体,分别对应因素的 个水平。从第 个总体抽取 个观测,单因素固定效应模型写作
其中 是第 组的总体均值, 是随机误差。总样本量 ,并要求 、。
经典精确 检验要求:
- 各组误差相互独立;
- ,即各组都服从正态分布;
- 各组方差相等,共用同一个 。
要检验的假设是
也可以把组均值写成
其中 是加权总均值, 是第 组的效应,原假设等价于所有 。后文记组均值和总均值分别为
平方和分解与检验
平方和分解
总平方和、组间平方和与组内平方和分别定义为
三者满足平方和分解
证明
把每个观测相对总均值的离差拆成“相对组均值”与“组均值相对总均值”两部分:
平方后对 求和,交叉项为
因为组内离差之和为零。于是对 求和得到
即 。
实际计算可用
若已经算出各组样本方差 ,则
各组样本量不等时,总均值要按 加权,不能把组均值简单平均。自由度也同样分解:
方差分析表与F检验
把平方和除以对应自由度得到均方,列成方差分析表:
| 变异来源 | 平方和 | 自由度 | 均方 | 统计量 |
|---|---|---|---|---|
| 组间(因素) | ||||
| 组内(误差) | ||||
| 总计 |
在 下,
若
则拒绝各组均值相等。 是上 分位点,记号与上侧分位点一致。
尽管均值差异可正可负,方差分析始终用右尾:差异被平方后反映为组间均方增大, 过大才说明组间变异相对组内误差过大。
是共同误差方差 的无偏估计; 不是。 接近 表示组间均方与组内均方相近;是否拒绝原假设仍要结合自由度与显著性水平判断。
使用与判断
- 多次做两样本 检验可能使整体误拒绝风险增加,不能替代整体 检验。
- 方差分析的目标是比较均值,名称来自用方差比检验均值差异。
- 是共同误差方差 的无偏估计,不是 。
- 残差图可辅助检查方差不齐与异常值;独立性还依赖试验设计,不能只凭正态概率图判断。
- 两组时,经典等方差单因素 检验与合并方差双侧 检验满足 。
组间和组内平方和的自由度分别为 和 。先将平方和除以自由度,再比较两个均方;与两个总体的方差比检验使用相同的分布族,但自由度来自不同的模型。
拒绝总体均值都相等的假设,只能说明至少有两个总体均值不同。确定具体哪些组有差异需要进一步比较,并控制多重比较带来的误拒绝风险。