统计百科 ·

Meta分析森林图与漏斗图的APA报告怎么写?合并效应量、I²异质性、发表偏倚完整模板

统计百科:Meta分析APA报告全解——森林图合并效应量与95%CI写法、漏斗图发表偏倚评估、Egger检验、I²异质性报告规范、随机/固定效应模型选择,给出可直接套用的APA 7th完整报告模板。

Meta分析的论文格式难点不在于跑出结果,而在于如何把森林图、漏斗图、I²检验、Egger检验这些结果按APA 7th规范组织成一段连贯的结果报告。本文专攻这一具体格式点:给出每个Meta分析核心指标的规范写法,以及可直接套用的完整报告段落模板。


你的Meta分析报告有这些问题吗?

审稿人或导师关于Meta分析结果报告最常见的反馈:

  • "合并效应量必须报告95%置信区间,不能只写SMD或OR值"
  • "异质性检验Q统计量和I²都要报告,只写I² = 68%不够"
  • "随机效应模型和固定效应模型的选择依据要在方法章节说明"
  • "漏斗图只呈现图片不够,Egger检验p值要写进结果正文"
  • "森林图中合并效应量的写法应与图注格式一致"
  • "I²的解读标准要引用Higgins et al.(2003),不能凭感觉描述"

这些问题集中在结果报告的规范性上。本文的目标只有一个:Meta分析每一个核心结果,到底怎么规范写进APA论文

如果你还没读过姊妹篇 Meta分析用AI一句话完成,建议先看那篇,了解从文献检索到效应量合并的全流程;本文默认你已有分析结果,专注讲报告格式。


一、Meta分析结果报告的整体框架

APA 7th对Meta分析报告没有单独章节,但结合PRISMA 2020声明与APA Publication Manual第七版的统计报告要求,结果章节通常按以下顺序组织:

顺序内容对应指标
1纳入研究概况研究数k、总样本量N
2模型选择依据固定/随机效应,Q检验p值
3合并效应量SMD/OR/RR/r + 95%CI + p值
4异质性评估Q统计量、I²、tau²
5发表偏倚漏斗图视觉检查 + Egger/Begg检验
6敏感性分析(若有)去掉影响力最大研究后的结果变化

每个部分都有固定的APA写法,下文逐一拆解。


二、森林图核心指标的APA写法

合并效应量(Pooled Effect Size)

合并效应量是Meta分析的核心结论,必须包含:效应量统计量类型 + 数值 + 95%CI + p值。

不同效应量指标的APA格式

效应量类型适用场景APA格式示例
SMD(标准化均值差)连续型结局,不同量表SMD = 0.52, 95% CI [0.31, 0.73], p < .001
MD(均值差)连续型结局,相同量表MD = 3.41, 95% CI [1.87, 4.95], p < .001
OR(比值比)二分类结局OR = 1.74, 95% CI [1.35, 2.24], p < .001
RR(相对危险度)队列研究二分类结局RR = 1.48, 95% CI [1.22, 1.80], p < .001
r(相关系数)相关性Metar = .34, 95% CI [.22, .45], p < .001

关键格式细节

  • 效应量数值:SMD/MD保留两位小数,OR/RR保留两位小数
  • 95% CI:方括号内两个数值,用逗号+空格分隔:[0.31, 0.73](不写"0.31至0.73")
  • p值:精确三位小数,不加前导零(.032而非0.032),p = .000改写为p < .001
  • SMD的Cohen(1988)效应大小参考:小效应 0.20,中效应 0.50,大效应 0.80

森林图的图注写法(APA Figure Caption)

APA 7th要求图注独立、自含。森林图图注标准格式:

Figure X
[数字]. [研究名称]Meta分析的森林图

注. 图中每个正方形代表单个研究的效应量估计,正方形大小与权重成正比;
横线代表95%置信区间;菱形代表合并效应量及其95%CI。
效应量为标准化均值差(SMD),正值表示[实验组优于对照组/自变量正向效应]。
采用随机效应模型(DerSimonian-Laird法),
合并SMD = 0.52, 95% CI [0.31, 0.73], p < .001。

三、异质性I²的APA写法

异质性检验三件套:Q + I² + tau²

APA格式要求同时报告Cochran's Q检验、I²和tau²(或tau),三者缺一不可。

格式模板

异质性检验结果显示,Cochran's Q([k-1]) = X.XX, p [值],
I² = XX%(XX%置信区间 [XX%, XX%]),tau² = X.XX(tau = X.XX),
表明各研究间存在[低/中/高]度异质性。

I²解读标准(引用Higgins et al., 2003)

I²范围异质性程度文字描述
0%–25%低度各研究结果较为一致,异质性低
25%–50%中度存在一定程度异质性
50%–75%高度异质性较高,建议开展亚组分析或Meta回归
>75%极高度异质性极高,合并效应量的解释需谨慎

报告I²时务必引用:Higgins, J. P. T., Thompson, S. G., Deeks, J. J., & Altman, D. G. (2003). Measuring inconsistency in meta-analyses. BMJ, 327(7414), 557–560.

模型选择依据的写法

鉴于纳入研究在被试特征、干预方案和测量工具上存在差异,
且Q检验结果显示存在显著异质性,Q([k-1]) = XX.XX, p < .001,
本研究采用随机效应模型(DerSimonian-Laird法)进行效应量合并。

若Q检验p > .05(无显著异质性),仍推荐报告I²,并可注明"虽Q检验未达显著,为保守起见仍采用随机效应模型"。


四、漏斗图与发表偏倚的APA写法

漏斗图图注(APA Figure Caption)

Figure X
[数字]. [研究名称]Meta分析的漏斗图

注. 横轴为效应量(SMD),纵轴为标准误(SE);
每个点代表一项纳入研究;虚线表示合并效应量;
两侧斜线为95%伪置信区间边界。
漏斗图呈[基本对称/轻度不对称/明显不对称],
提示发表偏倚[可能较小/待检验/可能存在]。

Egger检验结果的APA写法

Egger检验是定量评估发表偏倚的最常用方法,格式:

Egger's test: b = X.XX, SE = X.XX, t([k-2]) = X.XX, p = .XXX

完整发表偏倚段落模板

为评估发表偏倚,本研究绘制漏斗图并进行Egger线性回归检验。
漏斗图(见Figure X)呈[基本对称/轻度不对称],
Egger检验结果[显示截距不显著,b = X.XX, SE = X.XX, t([k-2]) = X.XX, p = .XXX,
提示发表偏倚可能较小/截距显著,b = X.XX, SE = X.XX, t([k-2]) = X.XX, p = .XXX,
提示可能存在一定程度的发表偏倚]。
[若Begg秩相关检验:Begg's test: Kendall's tau = X.XX, p = .XXX,
与Egger检验结果一致/结果不一致,综合判断发表偏倚...]

发表偏倚结论性语言

Egger检验p值漏斗图结论表述
p > .10,对称漏斗图对称未见显著发表偏倚证据(Egger, p = .XXX)
p < .10,轻度不对称漏斗图轻度不对称可能存在一定发表偏倚,结论需谨慎解读
p < .05,明显不对称漏斗图明显不对称存在显著发表偏倚,可结合Trim-and-Fill法校正

五、完整Meta分析结果报告段落模板

情形设定

纳入18项随机对照试验(k = 18,总N = 2,340),研究正念训练对焦虑症状的干预效果,效应量为SMD。

可直接套用的完整结果段落

最终纳入18项随机对照试验(k = 18),共涉及2,340名参与者
(实验组1,187人,对照组1,153人)。各研究样本量范围为42至312人,
干预周期为4至12周。

异质性检验结果显示,Cochran's Q(17) = 48.63, p < .001,
I² = 65.06%(95% CI [41.3%, 79.4%]),tau² = 0.073(tau = 0.270),
表明各研究间存在高度异质性(Higgins et al., 2003)。
鉴于此,本研究采用随机效应模型(DerSimonian-Laird法)合并效应量。

随机效应模型结果显示,正念训练对焦虑症状的合并效应量显著,
SMD = 0.52, 95% CI [0.31, 0.73], z = 4.89, p < .001,
为中等效应(Cohen, 1988),即正念训练组相比对照组,
焦虑症状改善幅度约为对照组标准差的0.52倍(见Figure 1)。

为评估发表偏倚,绘制漏斗图(见Figure 2)并进行Egger线性回归检验。
漏斗图视觉检查呈基本对称分布,Egger检验截距不显著,
b = 0.47, SE = 0.61, t(16) = 0.77, p = .451,
未见显著发表偏倚证据。

六、方法章节的规范描述模板

APA 7th要求方法章节说明Meta分析的统计方案,包括:

采用随机效应模型(DerSimonian-Laird法)合并各研究效应量,
以标准化均值差(SMD)为主要效应量指标,并计算95%置信区间。
采用Cochran's Q检验评估研究间异质性,并以I²统计量量化异质性程度;
I²解读参照Higgins et al.(2003)标准(0–25%为低度,25–50%为中度,50–75%为高度,>75%为极高度)。
发表偏倚通过漏斗图视觉检查及Egger线性回归检验(alpha = .10)评估。
所有Meta分析计算采用ChatSRS(R引擎,meta程序包)完成。

七、在ChatSRS一句话获得规范Meta分析报告

打开 chatsrs.com,上传效应量数据(含每项研究的均值、SD、样本量或OR/CI)后输入:

"对这18项随机对照试验的焦虑干预效果做Meta分析,效应量用SMD,采用随机效应模型(DerSimonian-Laird),输出:①合并SMD及95%CI和z/p值;②森林图;③Cochran's Q、I²(含95%CI)和tau²异质性检验;④漏斗图及Egger检验;⑤符合APA 7th格式的完整结果报告段落,含图注写法。"

ChatSRS将自动:选择DerSimonian-Laird随机效应模型 → 生成森林图(含各研究权重方块)→ 生成漏斗图(含95%伪置信区间)→ 输出APA格式完整段落(合并效应量+异质性+发表偏倚)→ 自动判断I²水平并给出Higgins等(2003)对应的文字描述。

如需按亚组(如干预类型、随访时长)分层分析,在指令中补充"按[变量名]做亚组分析"即可。


八、Meta分析APA报告常见错误对照表

常见错误正确写法违反的规范
SMD = 0.52(无CI)SMD = 0.52, 95% CI [0.31, 0.73]CI必须报告
95% CI (0.31–0.73)95% CI [0.31, 0.73]APA用方括号,不用圆括号
只报告I² = 65%,无Q检验Q(17) = 48.63, p < .001, I² = 65.1%Q与I²必须同时报告
I² = 65%"高异质性"(无引用)I² = 65%,高度异质性(Higgins et al., 2003)解读标准须引用文献
漏斗图只贴图,无Egger结果补充Egger t(k-2)值和p值定性图示须配定量检验
p = .000p < .001.000是舍入假值
合并效应为"固定效应"无依据说明选择依据(Q检验结果或研究设计异质性)模型选择须有方法论依据
tau²未报告tau² = 0.073(tau = 0.270)随机效应须报告组间方差
95% CI写成"0.31 to 0.73"95% CI [0.31, 0.73]APA格式为方括号
OR = 1.74(无CI无z/p)OR = 1.74, 95% CI [1.35, 2.24], p < .001效应量必须配齐CI和p值

常见 FAQ

Q:随机效应模型和固定效应模型,论文里怎么选、怎么写?

A:选择依据写在方法章节,结果章节只写使用了什么模型。固定效应模型假设所有研究估计的是同一个真实效应量,适用于研究设计高度同质、研究对象和干预几乎相同的场景;随机效应模型假设各研究的真实效应量来自某一分布,适用于纳入研究存在方法学差异的常见情形。实操判断:若Q检验显著(p < .10)或I² > 25%,应选随机效应模型;许多方法学家建议在系统综述中默认使用随机效应模型以获得更保守的估计。方法章节写法:"鉴于纳入研究在[被试/干预/测量]上存在差异,采用随机效应模型(DerSimonian-Laird法)"。

Q:I²报告时要不要给置信区间?

A:推荐给出,尤其是纳入研究数k较少(k < 20)时。I²的抽样方差在k较小时很大,点估计的不确定性高。APA格式:I² = 65.1%,95% CI [41.3%, 79.4%]。R的meta包和metafor包默认输出I²的置信区间,ChatSRS也自动输出。若期刊要求精简,至少保留点估计和Q检验。

Q:Egger检验显著了(p < .05),怎么处理,怎么写?

A:发现显著发表偏倚有两种处理方案,需在讨论章节说明局限性。方案一:使用Trim-and-Fill法估计"填充"缺失研究后的校正效应量,报告格式为"Trim-and-Fill法调整后的合并SMD = 0.38, 95% CI [0.19, 0.57],结论[保持不变/有所改变]";方案二:在讨论章节说明"可能存在发表偏倚"作为研究局限,结合灰色文献检索和注册研究情况进行综合判断。无论哪种方案,结果正文中Egger检验的写法不变:b = X.XX, SE = X.XX, t(k-2) = X.XX, p = .XXX。

Q:Meta分析只有8项研究,I²和Egger检验还有效吗?

A:k < 10时统计检验效力不足,需特别说明。Q检验在k较小时功效低,可能无法检出真实异质性(type II error);Egger检验在k < 10时也缺乏足够功效,结果可靠性有限。报告写法建议:"本研究纳入研究数k = 8,Q检验和Egger检验在样本量较小时功效有限(Egger et al., 1997),结果解读需谨慎"。可补充漏斗图视觉判断,并建议未来研究扩充文献量重新验证。

Q:森林图里的菱形(合并效应量)和每项研究的方块大小代表什么,图注怎么说明?

A:方块大小与该研究的权重成正比(随机效应模型中权重取决于研究内方差和研究间方差tau²),方块越大说明该研究对合并结果的贡献越大;横线是该研究效应量的95%CI;菱形的中心是合并效应量,左右顶点是合并效应量的95%CI边界。图注中建议统一说明:"方块大小与权重成正比;菱形表示合并效应量及其95%CI",避免读者误解。


快速参考:Meta分析APA报告格式速查卡

[纳入研究概况]
k = XX项研究,总N = XXXX;单项样本量范围 XX–XX

[模型选择(方法章节)]
随机效应模型(DerSimonian-Laird法)/ 固定效应模型
选择依据:Q([k-1]) = XX.XX, p [值],I² = XX%

[合并效应量(结果章节)]
SMD = X.XX, 95% CI [X.XX, X.XX], z = X.XX, p [值]
OR = X.XX, 95% CI [X.XX, X.XX], z = X.XX, p [值]
r = .XX, 95% CI [.XX, .XX], z = X.XX, p [值]

[异质性三件套]
Q([k-1]) = XX.XX, p [值]; I² = XX%(95% CI [XX%, XX%]); tau² = X.XX

[发表偏倚]
Egger: b = X.XX, SE = X.XX, t([k-2]) = X.XX, p = .XXX
Trim-and-Fill(若有):调整后 SMD = X.XX, 95% CI [X.XX, X.XX]

[I²解读标准(须引用Higgins et al., 2003)]
0-25%: 低度  25-50%: 中度  50-75%: 高度  >75%: 极高度

相关阅读


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