统计百科 ·

ANCOVA 协变量调整后均值(EMM)的 APA 报告怎么写?— 调整均值、F 检验、偏 eta 平方、协变量假设全解

统计百科:专攻 ANCOVA 协方差分析的 APA 7th 报告格式——协变量调整后均值(EMM/M_adj)、组别 F 检验、偏 eta 平方效应量、回归斜率齐性与线性假设检验,给出可直接套进论文的完整模板与常见错误对照。

ANCOVA 的 APA 报告难点不在跑软件,而在"调整后均值写哪儿、协变量假设检验怎么交代、偏 eta 平方和普通 eta 平方选哪个"这三件事。本文只讲 ANCOVA 报告这一具体格式,给出可抄进论文的完整 APA 7th 模板,并逐一列出审稿人最常挑的格式错误。


你的 ANCOVA 报告有这些问题吗?

审稿人或导师关于 ANCOVA 报告最常见的反馈:

  • "报告了原始均值,但 ANCOVA 应该报告协变量调整后的估计边际均值(EMM/M_adj)"
  • "协变量假设检验(回归斜率齐性、线性关系)没有交代"
  • "效应量写的是 eta 平方,ANCOVA 应报告 partial eta 平方"
  • "调整后均值没有给出标准误(SE)和置信区间"
  • "组别 F 检验自由度写错了,ANCOVA 的误差 df 需减去协变量个数"
  • "协变量本身是否显著没有报告"

这些问题集中在 APA Publication Manual 第七版关于协方差分析报告的具体规则上。本文专攻 ANCOVA 报告格式,默认你已了解 p 值不加前导零(.032 而非 0.032)等基础规则。

如果你还没看过姊妹篇 F 值的 APA 7th 报告写法,建议先读那篇,本文不再重复基础格式规则。


一、ANCOVA 核心概念回顾:为什么要报告调整后均值

什么是协变量(Covariate)

协变量是 ANCOVA 中用于"统计控制"的连续变量,通常是:

  • 前测成绩(控制入组基线差异)
  • 年龄、智力、焦虑程度等个体差异变量
  • 实验前的测量值(消除因随机分组不完美导致的组间偏差)

ANCOVA 的逻辑是:从因变量的总变异中"扣除"协变量所能解释的部分,再检验组别效应。这样,组别效应是在统计控制了协变量之后的净效应,比普通 ANOVA 更精准。

什么是调整后均值(EMM / M_adj)

估计边际均值(Estimated Marginal Means,EMM),又称调整后均值(Adjusted Means,M_adj),是将协变量固定在某一代表值(通常是全样本均值)后,各组因变量的预测均值。

$M_{adj} = \bar{Y}{group} - \hat{b} \cdot (\bar{X}{group} - \bar{X}_{grand})$

其中 $\hat{b}$ 是协变量的组内回归系数,$\bar{X}{group}$ 是该组协变量均值,$\bar{X}{grand}$ 是总体协变量均值。

为什么 ANCOVA 一定要报告调整后均值而非原始均值?

各组的原始均值已经包含了协变量差异的影响——如果实验组入组时成绩就更高,原始均值差异部分来自这个"起点优势",而非干预效果。调整后均值才真正代表"假设各组协变量相同时的组间差异",是 ANCOVA 的核心结论。报告原始均值而不报调整后均值,是 ANCOVA 报告的重大格式错误。


二、ANCOVA 的三大假设检验(先于主效应报告)

在报告组别主效应之前,APA 规范要求先报告 ANCOVA 的前提假设检验结果,通常包含三项:

假设一:回归斜率齐性(Homogeneity of Regression Slopes)

含义:协变量与因变量的回归关系在各组中斜率相同(不存在协变量 × 组别的交互效应)。若斜率不齐,说明协变量对各组的"调整量"不一样,ANCOVA 调整后均值的可比性存疑。

检验方法:在模型中加入"协变量 × 组别"交互项,检验该交互项是否显著。

APA 报告格式

回归斜率齐性假设以协变量与组别的交互项进行检验,结果不显著,
F(k-1, N-k-c-(k-1)) = X.XX, p = .XXX,满足斜率齐性假设。
(注:含 k-1 个交互项的扩充模型误差 df = N - k - c - (k-1);
当 k=2、c=1 时化简为 N-4,即 N-k-c-1。)

或不满足时:

回归斜率齐性假设检验显示交互项显著,F([df]) = X.XX, p = .XXX,
斜率齐性假设不成立;各组斜率分别为...,建议改用 Johnson-Neyman
分析或分组回归报告结果(Wilcox,2017)。

假设二:协变量与组别的独立性

含义:协变量不应受到实验操纵的影响(即组别不应显著预测协变量)。若协变量是"后测"而非"前测",则可能已经被干预影响,不适合作协变量。

检验方法:以组别为自变量、协变量为因变量做单因素 ANOVA 或 t 检验。

APA 报告格式

协变量独立性检验结果显示,两组在前测成绩上无显著差异,
t([df]) = X.XX, p = .XXX,满足协变量独立性假设。

假设三:协变量与因变量的线性关系

含义:协变量对因变量应具有线性预测关系;若是曲线关系,线性协变量调整会不准确。

检验方法:绘制散点图,或在模型中加入协变量的平方项检验。

APA 报告格式

散点图显示协变量(前测成绩)与因变量(后测成绩)之间呈线性关系,
满足 ANCOVA 线性假设。

三、自由度怎么算:ANCOVA 的 df 与 ANOVA 有何不同

ANCOVA 中,每加入一个协变量,误差自由度就减少 1,因为模型需要估计额外的回归系数。

设:

  • $k$ = 组数
  • $c$ = 协变量个数
  • $N$ = 总样本量
来源df
组别(Between Groups)$k - 1$
协变量(Covariate)$c$(每个协变量 1 个 df)
误差(Error)$N - k - c$
总计$N - 1$

单协变量的典型示例(2组,1个协变量,N = 74):

来源df
组别1
协变量1
误差71(= 74 - 2 - 1)
总计73

常见错误:把误差 df 写成 N - k(忽略协变量消耗的 df),导致 F 值自由度报告有误。


四、APA 7th 格式规范:ANCOVA 报告必须包含的要素

必报要素清单

要素格式示例
假设检验结果斜率齐性、协变量独立性分别报告 F 或 t,p 值
协变量效应F(1, N-k-c), p, partial eta^2F(1, 71) = 38.72, p < .001, partial eta^2 = .353
组别主效应F(k-1, N-k-c), p, partial eta^2F(1, 71) = 28.46, p < .001, partial eta^2 = .286
各组调整后均值M_adj = XX, SE = X.XXM_adj = 78.91, SE = 1.47
置信区间(可选但推荐)95% CI [LL, UL]95% CI [76.00, 81.82]
原始均值与 SD(对照用)M = XX, SD = XXM = 77.3, SD = 8.5

p 值与 F 值的精度要求

与所有 F 检验相同:

  • F 值精确到两位小数:F = 28.46
  • p 值精确到三位小数,不加前导零:p = .032
  • p = .000 改为 p < .001
  • 调整后均值报告两位小数:M_adj = 78.91

五、APA 7th ANCOVA 完整报告模板

模板一:单协变量、两组设计(最常见)

情形:实验组(n = 37)与对照组(n = 37),前测为协变量,后测为因变量。

完整报告段落模板

分析前,以协变量与组别的交互项检验回归斜率齐性假设,结果不显著,
F(1, [N-k-c-1]) = X.XX, p = .XXX,满足斜率齐性假设。
独立样本 t 检验显示两组前测成绩无显著差异,
t([N-k]) = X.XX, p = .XXX,满足协变量独立性假设。

以[协变量名]为协变量的单因素协方差分析(ANCOVA)结果如下:
协变量效应显著,F(1, [N-k-c]) = X.XX, p [值], partial eta^2 = .XX,
表明[协变量名]对[因变量名]具有显著预测作用。
控制[协变量名]后,组别主效应显著,
F([k-1], [N-k-c]) = X.XX, p [值], partial eta^2 = .XX,为[大小]效应
(Cohen,1988:.01 = 小效应,.06 = 中效应,.14 = 大效应)。

调整后均值比较显示,实验组(M_adj = XX, SE = X.XX, 95% CI [LL, UL])
显著[高于/低于]对照组(M_adj = XX, SE = X.XX, 95% CI [LL, UL])。

填入数值的示例

分析前,以协变量(前测成绩)与组别的交互项检验回归斜率齐性假设,结果不显著,F(1, 70) = 0.48,p = .490,满足斜率齐性假设。独立样本 t 检验显示,两组前测成绩无显著差异,t(72) = 0.37,p = .714(实验组:M = 63.2,SD = 8.1;对照组:M = 62.8,SD = 7.9),满足协变量独立性假设。

以前测成绩为协变量的单因素协方差分析(ANCOVA)结果显示,协变量效应显著,F(1, 71) = 38.72,p < .001,partial eta^2 = .353,表明前测成绩对后测成绩具有显著预测作用。控制前测成绩后,组别主效应显著,F(1, 71) = 28.46,p < .001,partial eta^2 = .286,为大效应(Cohen,1988)。调整后均值比较显示,实验组(M_adj = 78.91,SE = 1.47,95% CI [76.00, 81.82])显著高于对照组(M_adj = 67.23,SE = 1.53,95% CI [64.19, 70.27]),调整后均值差为 11.68(p < .001)。


模板二:单协变量、多组设计(三组及以上)

情形:三种干预方案(A/B/C),前测为协变量,后测为因变量,各组 n = 30,N = 90。

完整报告段落模板

回归斜率齐性假设检验不显著,F([k-1], [N-k-c-(k-1)]) = X.XX, p = .XXX,
(含 k-1 个交互项后,误差 df 在基础误差 N-k-c 基础上再减 k-1,即 N-k-c-(k-1))
协变量独立性假设通过单因素 ANOVA 确认,F([k-1], [N-k]) = X.XX, p = .XXX。

单因素协方差分析(One-Way ANCOVA)结果显示,
协变量效应显著,F(1, [N-k-1]) = X.XX, p [值], partial eta^2 = .XX;
控制协变量后,组别主效应显著,
F([k-1], [N-k-1]) = X.XX, p [值], partial eta^2 = .XX。

事后多重比较(Bonferroni 校正)显示:
[组A](M_adj = XX, SE = X.XX)显著[高于/低于][组B](M_adj = XX, SE = X.XX),
均值差 = X.XX(95% CI [XX, XX], p [值]);
[其余显著对比];[不显著对比]未达显著(p > .05)。

填入数值的示例

回归斜率齐性假设以交互项检验,未达显著,F(2, 84) = 0.62,p = .540。前测成绩的组间差异不显著,F(2, 87) = 0.14,p = .870,满足协变量独立性假设。

单因素协方差分析结果显示,协变量(前测成绩)效应显著,F(1, 86) = 42.31,p < .001,partial eta^2 = .330;控制前测后,干预方案主效应显著,F(2, 86) = 16.74,p < .001,partial eta^2 = .280,为大效应。Bonferroni 事后多重比较显示,方案 A(M_adj = 82.4,SE = 1.6)和方案 B(M_adj = 79.1,SE = 1.5)均显著高于方案 C(M_adj = 70.3,SE = 1.7),均值差分别为 12.1(p < .001)和 8.8(p = .002);方案 A 与方案 B 之间的差异未达显著(均值差 = 3.3,p = .180)。


模板三:多协变量设计

情形:两组,年龄(连续)和前测(连续)两个协变量,N = 120。

关键差异:误差 df = N - k - c = 120 - 2 - 2 = 116;每个协变量分别报告 F。

完整报告段落模板

以[协变量1]和[协变量2]为协变量的单因素协方差分析(ANCOVA)结果:
[协变量1]效应显著,F(1, [N-k-c]) = X.XX, p [值], partial eta^2 = .XX;
[协变量2]效应[显著/不显著],F(1, [N-k-c]) = X.XX, p = .XXX, partial eta^2 = .XX。
控制上述两个协变量后,组别主效应显著,
F([k-1], [N-k-c]) = X.XX, p [值], partial eta^2 = .XX。
实验组调整后均值(M_adj = XX,SE = X.XX)显著[高于/低于]
对照组(M_adj = XX,SE = X.XX)。

填入数值的示例

以年龄和前测成绩为协变量的单因素协方差分析显示,前测成绩效应显著,F(1, 116) = 44.18,p < .001,partial eta^2 = .276;年龄效应不显著,F(1, 116) = 1.94,p = .167,partial eta^2 = .016。控制上述两个协变量后,组别主效应显著,F(1, 116) = 22.37,p < .001,partial eta^2 = .162,为大效应。实验组调整后均值(M_adj = 76.4,SE = 1.31,95% CI [73.8, 79.0])显著高于对照组(M_adj = 67.8,SE = 1.38,95% CI [65.1, 70.5])。


六、方法章节如何描述 ANCOVA

标准方法章节描述(单协变量、两组)

采用单因素协方差分析(One-Way ANCOVA)检验[干预/组别]对[因变量]的效应,
以[协变量名]为协变量以控制入组基线差异。
分析前检验回归斜率齐性(协变量 × 组别交互项)和协变量独立性假设;
协变量与因变量之间的线性关系通过散点图确认。
效应量报告偏 eta 平方(partial eta^2),参照 Cohen(1988)标准
(.01 = 小效应,.06 = 中效应,.14 = 大效应)。
事后多重比较(若三组及以上)采用 Bonferroni 方法校正。
统计分析使用 ChatSRS(R 引擎),显著性水平设为 alpha = .05。

七、在 ChatSRS 一句话获得规范 ANCOVA 报告

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

"对实验组和对照组的后测成绩做单因素 ANCOVA,以前测成绩为协变量;先检验回归斜率齐性假设和协变量独立性假设,再做主分析;输出各组调整后均值(EMM)、SE 和 95% CI、协变量 F 值、组别 F 值(含两个自由度)、p 值、partial eta 平方效应量,并给出符合 APA 7th 格式的完整报告段落。"

ChatSRS 返回的段落即为上文示例格式,自动填入数值、交代假设检验、给出调整后均值和置信区间,可直接复制进论文。

如需三组以上 ANCOVA 或多个协变量,在指令中说明设计类型,ChatSRS 自动切换报告模板并处理 Bonferroni 事后比较。


八、ANCOVA 报告常见错误对照表

常见错误正确写法违反的规则
报告原始均值而非调整后均值必须报告 M_adj + SE + 95% CIANCOVA 核心输出是 EMM
不检验回归斜率齐性假设先报告斜率齐性交互项 F,p 值ANCOVA 前提假设必须交代
忽略协变量本身的 F 检验分别报告协变量 F 和组别 F协变量是否显著需单独报告
误差 df 未减去协变量 dfdf_error = N - k - c(而非 N - k)每个协变量消耗 1 个 df
效应量报告 eta^2(非偏)ANCOVA 必须报告 partial eta^2多因素/协变量设计用偏值
只说"调整后"不给数值M_adj = 78.91, SE = 1.47APA 要求报告具体数值
p = .000p < .001.000 是舍入假值
p = 0.032p = .032APA 不加前导零
忽略协变量独立性检验报告前测组间 t 检验或 ANOVA 结果协变量需先于干预存在

九、EMM 表格怎么呈现

APA 7th 建议在文字报告之外,用表格呈现各组的原始均值、调整后均值和相关信息。标准三线表格式:

表 X. 各组后测成绩的原始均值与协变量调整后估计边际均值

组别nM(原始)SDM_adjSE95% CI
实验组3777.38.578.911.47[76.00, 81.82]
对照组3768.19.267.231.53[64.19, 70.27]

注. M_adj = 协变量(前测成绩)固定在总样本均值(M = 63.0)后的估计边际均值。


常见 FAQ

Q:ANCOVA 的协变量一定要连续变量吗?可以用二分类变量作协变量吗?

A:ANCOVA 的协变量通常是连续变量(如前测、年龄、智力分数)。二分类变量(如性别)也可以加入模型,但此时 ANCOVA 实质上变成了包含组别 × 二分类变量的双因素模型,统计意义更接近于控制变量的双因素 ANOVA。如果确实需要控制性别影响,建议在方法章节明确说明"以性别为哑变量协变量",并分别报告性别效应的 F 值;不建议将性别称为"协变量"而不加说明,容易让读者误解为连续协变量。

Q:斜率齐性假设不成立时,ANCOVA 还能用吗?

A:若斜率齐性假设不成立(即协变量 × 组别交互项显著),标准 ANCOVA 的调整后均值会产生误导——调整量在各组不同,"控制协变量后的组别差异"本身因协变量值不同而不同。此时有几种处理方案:(1) 报告 Johnson-Neyman 分析,说明在协变量哪些取值范围内组别差异显著;(2) 在各协变量典型值处分别估计组别差异(pick-a-point 法);(3) 改用分层回归,将组别 × 协变量交互效应直接解释为调节效应。无论选哪种方案,均需在方法/结果章节明确说明假设不成立及其处理方式。

Q:ANCOVA 的 partial eta^2 和 ANOVA 的 partial eta^2 标准一样吗?

A:判断效应大小的 Cohen(1988)标准相同:.01 = 小效应,.06 = 中效应,.14 = 大效应。但值得注意的是,ANCOVA 通过控制协变量减小了误差方差(SS_error 变小),因此在实验设计不变的情况下,ANCOVA 得到的 partial eta^2 通常大于未控制协变量的 ANOVA 的 partial eta^2。这是 ANCOVA 提升统计功效的体现,不是计算错误。报告时应明确说明效应量是"控制[协变量]后的组别效应 partial eta^2",与纯 ANOVA 结果不可直接对比。

Q:调整后均值比原始均值高(或低)很多,正常吗?

A:正常,且差距大小取决于两个因素:(1) 各组在协变量上的差异大小;(2) 协变量与因变量的回归系数强度。当某组协变量均值明显低于总样本均值时,该组的调整后均值会高于原始均值("向上调整"),反之亦然。这正是 ANCOVA 的目的——消除因各组协变量基线不等带来的均值偏差。若调整前后均值差距极大,通常意味着各组基线差异较大,此时更应该在报告中同时呈现原始均值和调整后均值,让读者看到调整幅度,并在讨论中解释其含义。

Q:ANCOVA 结果需要做功效分析(power analysis)吗?

A:建议在方法章节简要报告。对于 ANCOVA,样本量计算通常比 ANOVA 更宽松——协变量的引入降低了误差方差,相当于提升了统计功效,因此达到同等功效所需的样本量更少。事后功效分析(post hoc power)可以在 G*Power 中选"ANCOVA: Fixed effects, main effects and interactions",输入实际 partial eta^2、alpha、样本量和协变量个数来计算观察功效(observed power),并在结果章节报告:如"基于 partial eta^2 = .286 的观察统计功效为 .94(alpha = .05)"。


相关阅读


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