『统计学』第五部分:方差分析和F检验

第四部分的卡方检验是研究类别变量之间的关系,而这一部分的方差分析则是研究类别型自变量与数值型因变量之间的关系,它在形式上是比较多个总体的均值是否相等

从形式上看,方差分析与之前的t检验z检验区别不大,都是检验均值是否相等,但在比较多个均值时,t检验需要做多次两两比较的假设检验,而方差分析只需要一次,并且方差分析中是将所有的样本信息结合在一起,增加了分析的可靠性。

下面我们仍旧从三个问题出发来研究一下方差分析:

  • 方差分析到底是啥?
  • 方差分析从哪里来?
  • 方差分析要到哪里去?(有什么用?)

这里并不打算按照这个顺序来分析这三个问题,因为很多理论都是为了解决某些实际问题而提出的,理论只是对解决方法的高度概括,从实际问题出发,反而可以更好地理解为什么要研究这些理论。下面将从以下几个方面来回答这些问题:

  • 1、为什么要做方差分析?(方差分析有什么用?)
  • 2、方差分析到底是啥?从哪里来的?
  • 3、方差分析的基本假定
  • 4、单因素方差分析
1、为什么要做方差分析?

这里举个实际的例子来引出方差分析,这样会对方差分析具体在做什么有更深的印象:

为了研究客户满意度是否与行业有关,消费者协会对各个行业抽样了数量不等的公司,分别统计了各个公司的被投诉次数:
在这里插入图片描述
那么基于以上的数据,你是否可以判断客户的投诉次数与行业有关。

问题其实比较清晰,就是比较这四个行业之间被投诉的均值是否相等,这里的行业就是类别型自变量,被投诉次数就是数值型因变量。所以方差分析要研究的就是行业对被投诉次数是否有显著影响。

2、方差分析到底是啥?从哪里来的?

方差分析的思想来自于误差的分解,对于来自同一个分布的数据,抽样带来的误差其实只有随机误差,即随机抽样抽到的不同值与均值之间的差距,以上面这个例子为例,我们可以把各个行业内部的抽样数据认为是来自同一个分布,即零售业不同企业有不同的被投诉次数,他们之间的不同可以认为只是随机误差,也被称为组内误差(SSE)

不同行业之间误差被称为组间误差(SSA),如果被投诉次数与行业无关,可以认为不同行业的抽样也是来自于同一总体分布的抽样,那么此时误差只有随机误差。而如果不同行业的抽样并不是来自同一总体分布,意味着被投诉次数与行业有关,那么此时组间误差包括随机误差和系统误差,这个系统误差就是由于行业这个因素带来。

综合以上,我们可以将总体误差(SST)分为组间误差和组内误差

当我们检验不同行业的均值是否相等时,其实就是检验组间误差是不是过大,包含了一定量的系统误差。

3、方差分析的基本假定

有三个基本假定:

每个总体都应服从正态分布;

每个总体的方差必须相同;

观测值是独立的;

其中第三点一般都满足,前两点需要借助样本值做一些检验的工作;

4、单因素方差分析(F检验)

如果只涉及一个类别的自变量的方差分析,则称其为单因素方差分析,最上面的例子就是一个单因素方差分析。

首先提出假设:

H 0 = μ 1 = μ 2 = . . . = μ k H_0 = \mu_1 = \mu_2 = ... = \mu_k H0=μ1=μ2=...=μk

H 1 = μ i ( i = 1... k ) H_1 = \mu_i (i = 1...k) H1=μi(i=1...k)不全相等

构造统计量:

1) 计算个样本的均值

即计算各个行业的均值:
在这里插入图片描述
2) 计算总均值

即计算总体均值:
在这里插入图片描述

这里:
X ‾ ‾ = 57 + 66 + . . . + 77 + 58 23 = 47.87 \overline {\overline X} = {57 + 66 + ... + 77 + 58 \over 23} = 47.87 X=2357+66+...+77+58=47.87
X ‾ 1 = 49 X ‾ 2 = 48 , X ‾ 3 = 35 , X ‾ 4 = 59 \overline X_1 = 49 \overline X_2 = 48, \overline X_3 = 35, \overline X_4 = 59 X1=49X2=48,X3=35,X4=59
求这四个平均值的平均值同样得到47.87。
所以,总均值可以看成所有数据点的均值,或每组均值的均值。

3) 计算各误差平方和

  • 总平方和:

S S T = ∑ i = 1 k ∑ j = 1 n i ( x i j − x ‾ ‾ ) 2 = ( 57 − 47.869 ) 2 + . . . + ( 58 − 47.47.869 ) 2 = 4164.608 SST = \sum^k_{i=1}\sum^{n_i}_{j=1}(x_{ij} - \overline {\overline x})^2 = (57 - 47.869)^2 + ... + (58 - 47.47.869)^2 = 4164.608 SST=i=1kj=1ni(xijx)2=(5747.869)2+...+(5847.47.869)2=4164.608

  • 组间平方和:

S S A = ∑ i = 1 k ( x ‾ i − x ‾ ‾ ) 2 = 7 ∗ ( 49 − 47.869 ) 2 + . . . + 5 ∗ ( 59 − 47.869 ) 2 = 1456.608 SSA = \sum^k_{i=1}(\overline x_{i} - \overline {\overline x})^2 = 7*(49 - 47.869)^2 + ... + 5*(59 - 47.869)^2 = 1456.608 SSA=i=1k(xix)2=7(4947.869)2+...+5(5947.869)2=1456.608

  • 组内平方差:

S S E = ∑ i = 1 k ∑ j = 1 n i ( x i j − x ‾ i ) 2 = ( 57 − 49 ) 2 + . . . + ( 58 − 59 ) 2 = 2708 SSE = \sum^k_{i=1}\sum^{n_i}_{j=1}(x_{ij} - \overline x_i)^2 = (57 - 49)^2 + ... + (58 - 59)^2 = 2708 SSE=i=1kj=1ni(xijxi)2=(5749)2+...+(5859)2=2708

并且: S S T = S S A + S S E SST = SSA + SSE SST=SSA+SSE

4) 计算统计量

SST的自由度为n-1,n为观测值的个数;

SSA的自由度为k-1,k为因素水平的个数;

SSE的自由度为n-k。

平方和自由度自由度计算公式
SST= 4164.60822n - 1
SSA = 1456.6083k - 1
SSE = 270819n - k

故:
总平方和 = 组内平方和 + 组间平方和

我们其实是要比较SSA与SSE,因为我们确认SSE,即组内误差来自随机误差,那我们就以SSE为基准,去比较SSA,若二者相差不大,则认为组间误差也仅仅是随机误差而没有系统误差,而由于二者之间的数量级不一样(自由度),难以直接比较,所以我们可以比较二者的均方误差:

组间均方误差: M S A = S S A k − 1 MSA = {SSA \over k-1} MSA=k1SSA 服从 χ 2 ( k − 1 ) \chi^2(k-1) χ2(k1)

组内均方误差: M S E = S S E n − k MSE = {SSE \over n-k} MSE=nkSSE 服从 χ 2 ( n − k ) \chi^2(n-k) χ2(nk)

于是,我们构造 F F F统计量: F = M S A M S E ∼ F ( k − 1 , n − k ) F = {MSA \over MSE} \sim F(k-1, n-k) F=MSEMSAF(k1,nk)

计算可得统计量的值为:

F = 3.406 F = 3.406 F=3.406

这里定义的统计量叫F统计量,它服从F分布,这里可以认为它是两个 χ 2 \chi^2 χ2分布之比,两者的自由度可以相等,也可以不等。
   \;
F统计量是组间平方和除以其自由度k-1,也被称为组间均方误差(MSA),然后除以组内平方和除以其自由度n-k:
F = M S A M S E = S S A k − 1 S S E n − k F = {MSA \over MSE} = {{SSA \over k-1} \over {SSE \over n-k}} F=MSEMSA=nkSSEk1SSA

如果分子比分母大很多,那就说明波动大多数来自于各组之间,较少来自于各组之内,这时我们应该相信,总体均值之间存在差异,如果这个数字很大,那就意味着零假设成立的概率较低;
如果这个数字很小,分母更大,那就意味着组内波动比组间波动在总波动中占比更多,这意味着差异可能只是随机产生的,这就更难拒绝零假设。

对于这个问题, α = 0.05 , n 1 = 3 , n 2 = 19 \alpha = 0.05, n_1 = 3, n_2 = 19 α=0.05,n1=3,n2=19,从下面的F分布表中可以读出临界F值:
在这里插入图片描述
这里从F分布统计表中读出的临界F值是:
F c = F 0.05 ( 3 , 19 ) = 3.13 F_c = F_{0.05}(3, 19) = 3.13 Fc=F0.05(3,19)=3.13

因为 F = 3.406 > F c = 3.13 F = 3.406 > F_c = 3.13 F=3.406>Fc=3.13,所以,在零假设前提下,得到观测值的概率非常低,因此我们拒绝零假设,相信很有可能总体均值存在差异。

参考文献:
[1] 方差分析
[2] 统计 可汗学院

  • 32
    点赞
  • 143
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
在医学统计学中,析因设计方差分析可以用来研究一个或多个自变量对因变量的影响,以及它们之间的交互作用。以下是析因设计方差分析的结论模板: 1. 对于主效应: 在α水平下,我们发现(自变量名称)显著影响(因变量名称)(F(自由度, F值) = F统计量, p < α),表明(自变量名称)的不同水平与(因变量名称)之间存在显著差异。 2. 对于交互作用: 在α水平下,我们发现(自变量1名称)和(自变量2名称)之间存在显著交互作用(F(自由度, F值) = F统计量, p < α),表明(自变量1名称)和(自变量2名称)在影响(因变量名称)方面不是独立的。 3. 对于主效应和交互作用的结合: 在α水平下,我们发现(自变量1名称)显著影响(因变量名称)(F(自由度, F值) = F统计量, p < α),而(自变量2名称)与(因变量名称)之间不存在显著差异(F(自由度, F值) = F统计量, p > α)。此外,我们还发现(自变量1名称)和(自变量2名称)之间存在显著交互作用(F(自由度, F值) = F统计量, p < α)。这表明(自变量1名称)和(自变量2名称)对(因变量名称)的影响不是独立的,它们之间的关系需要进一步探讨。 需要注意的是,对于每个结论,需要提供统计量(如F值)和显著性水平(如p值)的值,并解释它们的含义。此外,还需要根据实验的设计和问题来选择合适的结论模板。
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值