方差分析(2)

之前介绍了单因子方差分析的基本思想，从假设检验、给出检验统计量到最后的判断，以及使用Excel进行分析示例。对于多因子方差分析的基本思想与单因子方差分析类似，同样是利用方差和 F F F检验统计量，这里介绍的多因素方差分析主要是给出检验统计量并做出判断，以及用Excel示例。

双因子方差分析方法与单因子有相似，但有差异。考虑两个因素对结果的影响时，不仅有两个因素各自对结果的影响，还有两个因素交互作用产生的影响。下面列出无重复双因子方差分析和可重复上因子方差分析的数据格式和偏差平方和，以及方差分析表。

无重复双因子方差分析是考虑两个因子( A A A和 B B B)的变化对数据结果的影响，且两个因子不进行重复实验， 建立假设与单因素相似，即 H 0 H_0 H0​: r ∗ k r*k r∗k个总体均值相等， H 1 H_1 H1​： r ∗ k r*k r∗k个总体均值不相等。其数据结构如下表：

其中， x ˉ ⋅ 1 = ∑ i x i 1 r , x ˉ ⋅ 2 = ∑ i x i 2 r , . . . , x ˉ ⋅ k = ∑ i x i k r , ( i = 1 , 2 , . . . , r ) \bar{x}_{·1}=\frac{\sum_{i}{x_{i1}}}{r}, \bar{x}_{·2}=\frac{\sum_{i}{x_{i2}}}{r}, ..., \bar{x}_{·k}=\frac{\sum_{i}{x_{ik}}}{r}, (i=1,2,...,r) xˉ⋅1​=r∑i​xi1​​,xˉ⋅2​=r∑i​xi2​​,...,xˉ⋅k​=r∑i​xik​​,(i=1,2,...,r)

x ˉ 1 ⋅ = ∑ j x 1 j k , x ˉ 2 ⋅ = ∑ j x 2 j k , . . . , x ˉ r ⋅ = ∑ j x r j k , ( j = 1 , 2 , . . . , k ) \bar{x}_{1·}=\frac{\sum_{j}{x_{1j}}}{k}, \bar{x}_{2·}=\frac{\sum_{j}{x_{2j}}}{k}, ..., \bar{x}_{r·}=\frac{\sum_{j}{x_{rj}}}{k}, (j=1,2,...,k) xˉ1⋅​=k∑j​x1j​​,xˉ2⋅​=k∑j​x2j​​,...,xˉr⋅​=k∑j​xrj​​,(j=1,2,...,k)

总体均值 x ˉ ˉ \bar{\bar{x}} xˉˉ： x ˉ ˉ = ∑ i = 1 i = k x ˉ i ⋅ k = ∑ j = 1 j = r x ˉ ⋅ j r = ∑ i ∑ j x i j r k \bar{\bar{x}} = \frac{\sum_{i=1}^{i=k}\bar{x}_{i·}}{k} =\frac{\sum_{j=1}^{j=r}\bar{x}_{·j}}{r} = \frac{\sum_{i}{\sum_{j}{x_{ij}}}}{rk} xˉˉ=k∑i=1i=k​xˉi⋅​​=r∑j=1j=r​xˉ⋅j​​=rk∑i​∑j​xij​​

下面先给出对应因子的偏差平方和：

总体偏差平方和 S S T = ∑ j = 1 k ∑ i = 1 r ( x i j − x ˉ ˉ ) 2 , SS_T = \sum_{j=1}^{k}{\sum_{i=1}^{r}{(x_{ij}-\bar{\bar{x}})^2}}, SST​=j=1∑k​i=1∑r​(xij​−xˉˉ)2,自由度为 d f T = k r − 1 df_T=kr-1 dfT​=kr−1.

因子A偏差平方和 因子A在不同水平间的效应差异(主效应) S S A = ∑ j = 1 k ∑ i = 1 r ( x ˉ i ⋅ − x ˉ ˉ ) 2 = k ∑ i = 1 r ( x ˉ i ⋅ − x ˉ ˉ ) 2 , SS_A = \sum_{j=1}^{k}{\sum_{i=1}^{r}{(\bar{x}_{i·}-\bar{\bar{x}})^2}} =k{\sum_{i=1}^{r}{(\bar{x}_{i·}-\bar{\bar{x}})^2}}, SSA​=j=1∑k​i=1∑r​(xˉi⋅​−xˉˉ)2=ki=1∑r​(xˉi⋅​−xˉˉ)2,自由度为 d f A = r − 1 df_A=r-1 dfA​=r−1.

因子B偏差平方和 因子B在不同水平间的效应差异(主效应) S S B = ∑ j = 1 k ∑ i = 1 r ( x ˉ ⋅ j − x ˉ ˉ ) 2 = ∑ j = 1 k r ( x ˉ ⋅ j − x ˉ ˉ ) 2 , SS_B = \sum_{j=1}^{k}{\sum_{i=1}^{r}{(\bar{x}_{·j}-\bar{\bar{x}})^2}} ={\sum_{j=1}^{k}{r(\bar{x}_{·j}-\bar{\bar{x}})^2}}, SSB​=j=1∑k​i=1∑r​(xˉ⋅j​−xˉˉ)2=j=1∑k​r(xˉ⋅j​−xˉˉ)2,自由度为 k − 1 k-1 k−1.

随机误差的偏差平方和 S S e = ∑ j = 1 k ∑ i = 1 r ( x i j − x ˉ i ⋅ − x ˉ ⋅ j + x ˉ ˉ ) 2 , SS_e = \sum_{j=1}^{k}{\sum_{i=1}^{r}{(x_{ij}-\bar{x}_{i·}-\bar{x}_{·j} +\bar{\bar{x}})^2}}, SSe​=j=1∑k​i=1∑r​(xij​−xˉi⋅​−xˉ⋅j​+xˉˉ)2,自由度为 d f e = ( k − 1 ) ( r − 1 ) df_e=(k-1)(r-1) dfe​=(k−1)(r−1).

偏差平方和之间有如下等式: S S T = S S A + S S B + S S e SS_T= SS_A+SS_B+SS_e SST​=SSA​+SSB​+SSe​

可重复双因子方差分析是对两个因子的实验，且有同一水平下的重复实验，其数据格式如下(重复次数相同)：

其中，

所有 r ∗ k ∗ p r*k*p r∗k∗p个数据的平均： x ˉ ˉ = 1 r k p ∑ i = 1 r ∑ j = 1 k ∑ m = 1 p x i j m \bar{\bar{x}} = \frac{1}{rkp}{\sum_{i=1}^{r}\sum_{j=1}^{k}\sum_{m=1}^{p}{x_{ijm}}} xˉˉ=rkp1​i=1∑r​j=1∑k​m=1∑p​xijm​

A i A_i Ai​, B j B_j Bj​条件下, p p p次重复的平均： x ˉ i j ⋅ = 1 p ∑ m = 1 p x i j m , e g : x ˉ 11 ⋅ = 1 p ∑ m = 1 p x 11 m ， x ˉ 12 ⋅ = 1 p ∑ m = 1 p x 12 m \bar{x}_{ij·} = \frac{1}{p}{\sum_{m=1}^{p}{x_{ijm}}} , \\ eg: \bar{x}_{11·} = \frac{1}{p}{\sum_{m=1}^{p}{x_{11m}}} ， \bar{x}_{12·} = \frac{1}{p}{\sum_{m=1}^{p}{x_{12m}}} xˉij⋅​=p1​m=1∑p​xijm​,eg:xˉ11⋅​=p1​m=1∑p​x11m​，xˉ12⋅​=p1​m=1∑p​x12m​

A i A_i Ai​条件下， k ∗ p k*p k∗p个数据的平均: x ˉ i ⋅ ⋅ = 1 k p ∑ j = 1 k ∑ m = 1 p x i j m , e g : x ˉ 1 ⋅ ⋅ = 1 k p ∑ j = 1 k ∑ m = 1 p x 1 j m ,   x ˉ 2 ⋅ ⋅ = 1 k p ∑ j = 1 k ∑ m = 1 p x 2 j m \bar{x}_{i··} = \frac{1}{kp}{\sum_{j=1}^{k}\sum_{m=1}^{p}{x_{ijm}}}, \\ eg: \bar{x}_{1··} = \frac{1}{kp}{\sum_{j=1}^{k}\sum_{m=1}^{p}{x_{1jm}}}, \ \bar{x}_{2··} = \frac{1}{kp}{\sum_{j=1}^{k}\sum_{m=1}^{p}{x_{2jm}}} xˉi⋅⋅​=kp1​j=1∑k​m=1∑p​xijm​,eg:xˉ1⋅⋅​=kp1​j=1∑k​m=1∑p​x1jm​, xˉ2⋅⋅​=kp1​j=1∑k​m=1∑p​x2jm​

B j B_j Bj​条件下， r ∗ p r*p r∗p个数据的平均： x ˉ ⋅ j ⋅ = 1 r p ∑ i = 1 r ∑ m = 1 p x i j m , e g : x ˉ ⋅ 1 ⋅ = 1 r p ∑ i = 1 r ∑ m = 1 p x i 1 m ,   x ˉ ⋅ 2 ⋅ = 1 r p ∑ i = 1 r ∑ m = 1 p x i 2 m \bar{x}_{·j·} = \frac{1}{rp}{\sum_{i=1}^{r}\sum_{m=1}^{p}{x_{ijm}}}, \\ eg: \bar{x}_{·1·} = \frac{1}{rp}{\sum_{i=1}^{r}\sum_{m=1}^{p}{x_{i1m}}}, \ \bar{x}_{·2·} = \frac{1}{rp}{\sum_{i=1}^{r}\sum_{m=1}^{p}{x_{i2m}}} xˉ⋅j⋅​=rp1​i=1∑r​m=1∑p​xijm​,eg:xˉ⋅1⋅​=rp1​i=1∑r​m=1∑p​xi1m​, xˉ⋅2⋅​=rp1​i=1∑r​m=1∑p​xi2m​

注：上述表中，列出数据是以等重复次数列出，对于重复次数不等的情形，计算均值方法类似。

总体偏差平方和 S S T = ∑ i = 1 r ∑ j = 1 k ∑ m = 1 p ( x i j m − x ˉ ˉ ) 2 ,   自 由 度 d f T = r k p − 1 SS_T = {\sum_{i=1}^{r}\sum_{j=1}^{k}\sum_{m=1}^{p}{(x_{ijm} - \bar{\bar{x}})^2}}, \ 自由度df_T = rkp -1 SST​=i=1∑r​j=1∑k​m=1∑p​(xijm​−xˉˉ)2, 自由度dfT​=rkp−1

因子 A A A偏差平方和 S S A = ∑ i = 1 r ∑ j = 1 k ∑ m = 1 p ( x i ⋅ ⋅ − x ˉ ˉ ) 2 = k p ∑ i = 1 r ( x i ⋅ ⋅ − x ˉ ˉ ) 2 SS_A = {\sum_{i=1}^{r}\sum_{j=1}^{k}\sum_{m=1}^{p}{(x_{i··} - \bar{\bar{x}})^2}} =kp\sum_{i=1}^{r}(x_{i··} - \bar{\bar{x}})^2 SSA​=i=1∑r​j=1∑k​m=1∑p​(xi⋅⋅​−xˉˉ)2=kpi=1∑r​(xi⋅⋅​−xˉˉ)2 自由度为 d f B = r − 1 df_B = r-1 dfB​=r−1.

因子 B B B偏差平方和 S S B = ∑ i = 1 r ∑ j = 1 k ∑ m = 1 p ( x ⋅ j ⋅ − x ˉ ˉ ) 2 = r p ∑ j = 1 k ( x ⋅ j ⋅ − x ˉ ˉ ) 2 , SS_B = {\sum_{i=1}^{r}\sum_{j=1}^{k}\sum_{m=1}^{p}{(x_{·j·} - \bar{\bar{x}})^2}} =rp\sum_{j=1}^{k}{(x_{·j·} - \bar{\bar{x}})^2}, SSB​=i=1∑r​j=1∑k​m=1∑p​(x⋅j⋅​−xˉˉ)2=rpj=1∑k​(x⋅j⋅​−xˉˉ)2, 自由度为 d f B = k − 1 df_B = k-1 dfB​=k−1.

随机误差的偏差平方和 S S e = ∑ i = 1 r ∑ j = 1 k ∑ m = 1 p ( x i j m − x ˉ i j ⋅ ) 2 , SS_e = {\sum_{i=1}^{r}\sum_{j=1}^{k}\sum_{m=1}^{p}{(x_{ijm} - \bar{x}_{ij·})^2}}, SSe​=i=1∑r​j=1∑k​m=1∑p​(xijm​−xˉij⋅​)2, 自由度为 d f e = r k ( p − 1 ) df_e=rk(p-1) dfe​=rk(p−1).

因子 A B AB AB的交互影响偏差平方和 S S A B = S S T − S S A − S S B = ∑ i = 1 r ∑ j = 1 k ∑ m = 1 p [ ( x ˉ ˉ + x ˉ i j ⋅ ) − ( x i ⋅ ⋅ + x ⋅ j ⋅ ) ] 2 , SS_{AB} = SS_T - SS_A -SS_B = \sum_{i=1}^{r}\sum_{j=1}^{k}\sum_{m=1}^{p}[(\bar{\bar{x}} + \bar{x}_{ij·}) - (x_{i··} + x_{·j·}) ]^2, SSAB​=SST​−SSA​−SSB​=i=1∑r​j=1∑k​m=1∑p​[(xˉˉ+xˉij⋅​)−(xi⋅⋅​+x⋅j⋅​)]2, 自由度为 d f A B = d f T − d f A − d f B = ( r − 1 ) ( k − 1 ) df_{AB}=df_T - df_A -df_B = (r-1)(k-1) dfAB​=dfT​−dfA​−dfB​=(r−1)(k−1).

对应偏差平方和的均方（ M S = S S / d f MS = SS/df MS=SS/df），分别记为 M S A , M S B , M S e , M S A B MS_A, MS_B, MS_e, MS_{AB} MSA​,MSB​,MSe​,MSAB​。 与单因子方差分析相似，仍使用 F F F检验统计量。

当交互作用不显著时，可用随机误差平方和与交互作用偏差平方和作为新误差平方和： S S E = S S A B + S S e , SS_E = SS_{AB} + SS_e, SSE​=SSAB​+SSe​, 自由度 d f E = d f A B + d f e df_E =df_{AB} + df_e dfE​=dfAB​+dfe​, 均方为 M S E = S S E / d f E MS_E = SS_E/ df_E MSE​=SSE​/dfE​. 则对应 A A A, B B B的 F F F值和 p p p的计算则为： F A = M S A / M S E , p = F α ( d f A , d f E ) ; F B = M S B / M S E , p = F α ( d f B , d f E ) . F_A= MS_A / MS_E, p = F_\alpha(df_A, df_E); \\ F_B = MS_B / MS_E, p = F_\alpha(df_B, df_E). FA​=MSA​/MSE​,p=Fα​(dfA​,dfE​);FB​=MSB​/MSE​,p=Fα​(dfB​,dfE​).

Excel版本： Microsoft Office 2016

类似无重复双因子方差分析，选择"可重复上因素分析"工具：

选择数据，填写"样本行数"（示例中对因子A进行3次重复）、置信水平和输出区域：

上一步确定后，生成如下结果： 上面结果中，差异源中样本、列、交互和内部分别对应的是因子 A A A，因子 B B B、 A B AB AB交互影响和随机误差 e e e的信息。各列中SS即为对应偏差平方和，df为对应自由度，MS为对应均方，F为对应 F F F值，P-value 为对应 p p p值，F-crit为 F F F临界值。

与无重复双因子方差分析相同，其 p p p值和 F F F临界值也可通过Excel中函数得到： =F.DIST.RT({F_value},df1,df2) 和 =F.INV(1-{alpha},df1,df2)。

2019-07- 26