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