统计百科 ·

偏η²(partial eta squared)怎么算怎么报?公式推导 + APA写法全解

统计百科:偏η²的核心公式 SSeffect/(SSeffect+SSerror)、手动计算步骤、单因素/双因素/重复测量三场景数值示例、与η²的数值差距、APA 7th可抄进论文的完整报告句式,配Cohen效应量阈值速查。

审稿意见写着"请报告偏 eta 平方",SPSS 输出了一列数字却不知道怎么来的,论文效应量一行写成 η² = .267 结果被要求改成偏η²……偏η² 的名字人人熟悉,但真正能从 SS 表格手算一遍、说清楚它和η² 在数值上差多少、在 APA 7th 里每种设计该怎么写的,并不多。本文从公式出发,逐步拆解偏η² 的计算逻辑,给出三类设计的完整 APA 报告句式,每一句都可以直接套进论文。


一、偏η²的核心公式:从哪里来,为什么这样设计

公式

$\text{partial }\eta^2 = \frac{SS_{\text{effect}}}{SS_{\text{effect}} + SS_{\text{error}}}$

三个符号逐一解释:

符号全称含义
SS_effectSum of Squares for Effect目标效应(某主效应或交互效应)的平方和
SS_errorSum of Squares for Error**误差(组内/残差)**的平方和
partial η²偏 eta 平方目标效应在"效应+误差"这个局部空间里占的比例

设计意图:分母里只保留"目标效应 + 误差",其他因素的 SS 全部排除在外。这样得到的比例反映的是:在控制所有其他变量后,该效应解释了多少"剩余可解释方差"

与η²的公式对比

$\eta^2 = \frac{SS_{\text{effect}}}{SS_{\text{total}}} \quad \text{vs.} \quad \text{partial }\eta^2 = \frac{SS_{\text{effect}}}{SS_{\text{effect}} + SS_{\text{error}}}$

两者分子相同,区别只在分母:

  • η² 的分母是 SS_total(全部方差,包含所有主效应、交互效应、误差)
  • 偏η² 的分母是 SS_effect + SS_error(只有目标效应和误差,其他效应被"排除")

结果:偏η² 的分母 ≤ SS_total,所以 偏η² ≥ η²,在多因素设计中两者差距会很显著。


二、手动计算:从ANOVA摘要表逐步算出偏η²

读懂ANOVA摘要表

拿到 SPSS / R / Stata 的方差分析输出后,找到这几列:

变异来源SSdfMSFp
因素ASS_Adf_AMS_AF_Ap_A
因素BSS_Bdf_BMS_BF_Bp_B
A × B 交互SS_AxBdf_AxBMS_AxBF_AxBp_AxB
误差(组内)SS_errordf_errorMS_error
总计SS_total

偏η² 的计算只需要 SS_effect(目标行)和 SS_error(误差行)两个数,SS_total 不用。

数值示例:双因素ANOVA

研究情境:性别(2水平)× 学习策略(4水平)的学业成绩分析,N = 200。

ANOVA 摘要表:

变异来源SSdf
性别(A)1201
学习策略(B)4803
性别 × 学习策略(A×B)903
误差(组内)810192
总计1500199

逐行计算偏η²:

性别主效应

$\text{partial }\eta^2_A = \frac{120}{120 + 810} = \frac{120}{930} = .129$

学习策略主效应

$\text{partial }\eta^2_B = \frac{480}{480 + 810} = \frac{480}{1290} = .372$

交互效应

$\text{partial }\eta^2_{A \times B} = \frac{90}{90 + 810} = \frac{90}{900} = .100$

同一数据的η²(对比)

$\eta^2_A = \frac{120}{1500} = .080 \quad \eta^2_B = \frac{480}{1500} = .320 \quad \eta^2_{A\times B} = \frac{90}{1500} = .060$

差距:学习策略的偏η² = .372,η² = .320,相差 0.052,不是小数点舍入误差。这就是为什么多因素设计必须报告偏η²,而不是η²。


三、三类设计的偏η²计算逻辑

单因素ANOVA(特殊情形:偏η² = η²)

$SS_{\text{total}} = SS_{\text{between}} + SS_{\text{within}}$

此时 SS_error = SS_within,且没有其他因素的 SS,因此:

$\text{partial }\eta^2 = \frac{SS_{\text{between}}}{SS_{\text{between}} + SS_{\text{within}}} = \frac{SS_{\text{between}}}{SS_{\text{total}}} = \eta^2$

结论:单因素 ANOVA 里偏η² 和η² 数值完全相同,报告哪个都正确,但必须写清楚用的是哪个。

双因素/多因素ANOVA

每个效应的偏η² 独立计算,分母永远是"该效应的 SS + 误差 SS":

$\text{partial }\eta^2_{\text{任意效应}} = \frac{SS_{\text{该效应}}}{SS_{\text{该效应}} + SS_{\text{error}}}$

注意:每个效应共用同一个 SS_error,这导致各效应的偏η² 之和可能超过 1——这是正常现象,不是计算错误。

重复测量ANOVA

重复测量设计中误差项是分层的,不同效应对应不同误差项:

效应SS_effectSS_error
时间(被试内因素)SS_timeSS_time×subjects(时间×被试交互)
组别(被试间因素)SS_groupSS_subjects-within-group(组内被试间误差)
时间 × 组别交互SS_time×groupSS_time×group×subjects(三阶交互误差)

SPSS / R 会自动配对正确的误差项,直接读输出的 Partial Eta Squared 即可,不需要手动计算。


四、Cohen效应量阈值速查

Cohen(1988)基于单因素 ANOVA 的η² 给出阈值,实践中同等适用于偏η²:

效应大小偏η² 值直觉参照
小效应.01统计可见,肉眼几乎无法察觉
中效应.06仔细观察才能发现的差异
大效应.14肉眼可见的明显差异

实操提示:偏η² 因为分母比η² 小,数值往往略高,用 Cohen 阈值解读时会略微高估效应大小,但目前学界普遍接受这种做法。报告时说明"参照 Cohen(1988)标准"即可。


五、APA 7th 偏η² 报告句式——可直接抄进论文

格式规则

  1. 偏η² 写法:partial η² = .XX(两位小数,不加前导零);行文中可写"偏η²"、"partial η²"或"偏 eta 平方"
  2. 必须与 F 值和 p 值同行报告:F(df1, df2) = X.XX, p [值], partial η² = .XX
  3. 效应大小解读:加括号注明 (Cohen, 1988)
  4. 方法章节需说明使用了偏η² 而非η²

单因素ANOVA(η² = 偏η²,报告η²或说明等价)

单因素方差分析结果显示,[自变量]对[因变量]的主效应显著,
F([k-1], [N-k]) = X.XX, p [= .XXX / < .001], η² = .XX,
为[小/中/大]效应(Cohen, 1988)。

填入数值示例

单因素方差分析结果显示,教学方法对学生期末成绩的主效应显著,F(3, 196) = 23.78, p < .001, η² = .267,为大效应(Cohen, 1988)。Tukey HSD 事后多重比较显示,项目式教学(M = 88.9, SD = 8.4)显著优于传统讲授(M = 76.4, SD = 8.2),均值差 = 12.5(95% CI [9.2, 15.8], p < .001)。


双因素ANOVA(必须报告偏η²)

双因素方差分析结果显示:
[因素A]主效应显著,F([a-1], [N-ab]) = X.XX, p [值], 偏η² = .XX,
为[效应大小](Cohen, 1988);
[因素B]主效应显著,F([b-1], [N-ab]) = X.XX, p [值], 偏η² = .XX,
为[效应大小];
[A × B]交互效应[显著/不显著],F([(a-1)(b-1)], [N-ab]) = X.XX,
p [值], 偏η² = .XX。

填入数值示例(使用上方示例数据)

双因素方差分析结果显示,性别主效应显著,F(1, 192) = 28.44, p < .001, 偏η² = .129,为中效应(Cohen, 1988);学习策略主效应显著,F(3, 192) = 37.93, p < .001, 偏η² = .372,为大效应;性别 × 学习策略交互效应显著,F(3, 192) = 7.11, p < .001, 偏η² = .100,为中效应。进一步的简单效应分析显示,女生的学习策略效应(偏η² = .421)大于男生(偏η² = .318),表明女生成绩对学习策略变化更为敏感。


重复测量ANOVA(报告偏η²,需附球形检验)

重复测量方差分析前进行 Mauchly 球形假设检验,
结果[显著/不显著(W = X.XX, χ²(df) = X.XX, p [值])],
[违反球形假设,采用 Greenhouse-Geisser 校正(ε = X.XX)/
球形假设成立,自由度未作调整]。
时间主效应显著,F([校正df1], [校正df2]) = X.XX, p [值],
偏η² = .XX,为[效应大小]效应(Cohen, 1988)。

填入数值示例

重复测量方差分析前进行 Mauchly 球形假设检验,结果显著(W = 0.487, χ²(5) = 56.32, p < .001),违反球形假设,采用 Greenhouse-Geisser 方法对自由度进行校正(ε = 0.685)。时间主效应显著,F(2.05, 162.15) = 142.7, p < .001, 偏η² = .643,为大效应(Cohen, 1988),表明被试焦虑水平在 12 周干预期间发生了显著变化。


ANCOVA(协变量分析,报告控制后偏η²)

以[协变量]为协变量进行单因素协方差分析(ANCOVA)。
协变量对[因变量]的效应显著,F(1, [N-k-1]) = X.XX, p [值], 偏η² = .XX。
控制[协变量]后,组别主效应显著,F([k-1], [N-k-1]) = X.XX, p [值],
偏η² = .XX,为[效应大小]效应(Cohen, 1988)。

填入数值示例

以前测成绩为协变量进行单因素协方差分析(ANCOVA)。前测成绩对后测成绩的效应显著,F(1, 70) = 38.72, p < .001, 偏η² = .356。控制前测成绩后,教学方法的主效应显著,F(1, 70) = 28.46, p < .001, 偏η² = .289,为大效应(Cohen, 1988),实验组调整后均值(M_adj = 78.91, SE = 1.47)显著高于对照组(M_adj = 67.23, SE = 1.53)。


方法章节说明效应量指标(直接复制)

多因素/重复测量/ANCOVA 设计

效应量以偏 eta 平方(partial η²)报告,以反映每个效应在
控制其他因素后的独立贡献;参照 Cohen(1988)标准评估效应大小
(小效应 .01,中效应 .06,大效应 .14)。统计分析使用
ChatSRS(SPSS+R+Stata 三引擎),显著性水平设为 α = .05。

六、在chatsrs.com一句话获得偏η²及APA报告

打开 chatsrs.com,上传数据后输入:

"对性别(2组)× 学习策略(4组)做双因素ANOVA,输出各主效应和交互效应的F值(含两个自由度)、p值、偏eta平方;按Cohen(1988)标准标注每个效应的大小;给出完整APA 7th报告段落,方法章节注明使用了偏eta平方而非eta平方,同时输出各效应对应的SS和误差SS以供核验。"

ChatSRS 会自动识别设计类型,单因素时输出η²,多因素时切换为偏η²,并在输出段落中精确标注使用了哪个指标。报告段落包含完整的统计数字,可直接复制进论文。

如需同时输出无偏估计,加一句"同时输出偏omega平方作为偏η²的保守对照"即可。


七、常见计算陷阱与检查清单

陷阱1:把η²当偏η²报告(最常见错误)

SPSS GLM 输出的 Partial Eta Squared 是偏η²,数值比η² 偏高。若将 SPSS 给的 .372 写成"η² = .372",实际上η² = .320,差距 0.052,审稿人核算后会要求修改。

检查方法:手算 SS_effect / SS_total,与 SPSS 给的值比较;不相等则 SPSS 给的是偏η²,应标注为偏η²。

陷阱2:多个效应的偏η²相加超过1

这不是错误。每个效应的偏η² 分母不同(共用同一 SS_error,但分子各异),不能把多个效应的偏η² 简单相加作概率解读。

陷阱3:重复测量设计用了错误的误差项

重复测量 ANOVA 的不同效应(被试内 vs. 被试间)对应不同误差项。手算时需确认用的是哪个 SS_error。直接读 SPSS / R 输出的 Partial Eta Squared,避免手算错配。

陷阱4:p < .001 仍然写成 p = .000

偏η² 与 p 值必须同行。如果软件显示 p = .000,报告时改写为 p < .001。

偏η² 报告自检清单

[  ] 分母确认为 SS_effect + SS_error(非 SS_total)
[  ] 精确到两位小数,不加前导零:偏η² = .267
[  ] F 值行同时报告:F(df1, df2) = X.XX, p [值], 偏η² = .XX
[  ] 注明效应大小并引用 Cohen(1988)
[  ] 方法章节说明"效应量以偏η²报告"
[  ] 多因素设计未将 SPSS 输出误写为 η²

常见FAQ

Q:SPSS 输出了 Partial Eta Squared,直接报告就行吗,还需要手算验证?

A:直接报告即可,SPSS 的 Partial Eta Squared 正是公式 SS_effect / (SS_effect + SS_error) 的计算结果。如果想核验,勾选 SPSS 输出中的"SS"列,手动验算一行即可。多数情况下 SPSS 数值是正确的,手算验证主要是帮助理解公式来源,不是必须步骤。

Q:偏η² = .05 是"小效应"还是"中效应"?介于两个阈值之间怎么写?

A:.05 在 Cohen 的小效应(.01)和中效应(.06)之间,更接近中效应。报告时可写"接近中等效应"(approaching medium effect)或直接给出数值不做定性,如 偏η² = .05,让读者自行对照 Cohen 阈值。不必强行归入某个档次。

Q:eta generalized(广义η²)和偏η² 有什么区别?

A:广义η²(eta_G²)是 Olejnik & Algina(2003)提出的改进版效应量,分母包含所有测量误差但排除设计误差,旨在提高跨研究可比性。偏η² 是最主流的期刊报告指标,广义η² 主要出现在心理学方法学论文和元分析比较中。除非目标期刊或审稿人明确要求广义η²,否则报告偏η² 即可。

Q:样本量很小(N = 30)时,偏η² 还可靠吗?

A:小样本下偏η² 对总体效应量有较明显的高估。N < 100 时推荐同时报告偏ω²(partial omega squared)作为无偏对照。偏ω² 的计算比偏η² 复杂,ChatSRS 可以一键输出两者并给出对照解读。

Q:不显著的效应(p > .05)也需要报告偏η² 吗?

A:必须报告。p > .05 只说明无法拒绝零假设,不等于效应量为零。报告偏η²(即使 .01 或更小),结合样本量和统计功效(power),读者才能判断是真的无效应,还是因样本不足导致功效低下。APA 7th 明确要求所有统计检验均须报告效应量,不区分显著与否。


相关阅读


本文首发于 ChatSRS 官方博客。如有问题或反馈,欢迎在 ChatSRS 站内 [用户中心 -> 帮助与反馈] 联系我们。