单个率Meta分析结果解读
一、方法概述
Meta分析(Meta-analysis)是一种对多项独立研究结果进行定量合并的统计分析方法,由Glass于1976年首次提出并命名。单个率Meta分析是Meta分析的基本形式之一,其目标是将多项独立研究中报告的同一事件发生率(或患病率、阳性率等)进行加权合并,从而获得更为精确的总体率估计及其置信区间。该方法在流行病学、循证医学及公共卫生领域具有广泛应用,尤其适用于整合不同地区、不同时间或不同样本来源的率数据,以克服单项研究样本量有限的局限。SPSSAU软件提供了完整的单个率Meta分析流程,涵盖效应量合并、异质性检验、发表偏倚检测、敏感性分析及累积Meta分析等核心模块。
在方法选择上,单个率Meta分析通常采用随机效应模型(Random-effects Model),以考虑不同研究间可能存在的真实效应差异。本研究采用DerSimonian-Laird(DL)法估计研究间方差τ²,并基于倒方差法(Inverse Variance Method)计算各研究权重。同时采用Hartung-Knapp方法对合并效应量的置信区间进行校正,以提高小样本合并时的覆盖概率。合并指标为率Ratio,即各研究的事件发生率。
在SPSSAU【Meta荟萃分析】模块选择【单个率】,按提示粘贴或修改数据,操作如下图:
二、分析基本信息
表1报告了本次Meta分析的基本设定。本研究共纳入5项独立研究,采用随机效应模型进行效应量合并,tau²估计方法为DL法,效应量计算采用倒方差法。合并指标为率Ratio,并启用Hartung-Knapp校正以提高置信区间的准确性。同时执行了发表偏倚检验(Egger检验和Begg检验)、敏感性检验(逐一剔除法)和累积Meta分析,以全面评估合并效应的稳健性。
三、效应量合并结果
表2报告了各单项研究及合并效应量的结果。五项研究的样本量分别为82、181、46、142和104,合计555例。各研究的率估计值介于0.413~0.486之间,波动范围较小,表明不同研究间的事件发生率具有较好的一致性。权重分配方面,AHS 2018的权重最高(32.924%),其次为CARET 2010(25.187%)和Weinstein 2015(19.001%),ATBC 2011的权重最低(8.215%)。权重分布与各研究的样本量基本对应,样本量最大的AHS 2018(n=181)获得最高权重,样本量最小的ATBC 2011(n=46)权重最低,符合倒方差法的权重分配逻辑。
合并效应量结果显示,总体合并率为0.443(95% CI: 0.403~0.483),z=30.739,p<0.001,达到极其显著的统计学水平。置信区间宽度为0.080个单位,区间较窄且不跨越无效值,表明合并效应量的估计精度较高。从合并率的数值来看,0.443意味着在综合五项研究后,目标事件的总体发生率约为44.3%。
四、异质性检验
表3报告了异质性检验结果。Cochran's Q统计量为1.874(p=0.759),未达到显著水平,表明各研究间的效应量差异可以归因于抽样误差,不存在显著的统计学异质性。I²=0.000%,远低于通常所采用的25%低异质性阈值,进一步证实了研究间一致性良好。tau²=0.000,表明研究间方差的估计值为零,随机效应模型的权重分配实际上等同于固定效应模型。H=1.000(95% CI: 1.000~2.193),H²=1.000,均处于理论最低值,表明研究间不存在过度离散。综合各项异质性指标,当前Meta分析的证据基础具有高度一致性,合并效应量的可靠性较强。
五、发表偏倚检验
表4报告了两种发表偏倚检验的结果。基于回归的Egger检验中,截距beta₁=0.517,t=0.327(df=3),p=0.765,远大于0.05的显著性阈值,表明漏斗图不存在显著的不对称性。非参数秩相关的Begg检验中,Kendall's score=4.000,z=0.980,p=0.327,同样未达到显著水平。两种方法的结果一致表明,当前纳入的五项研究不存在可检测的发表偏倚,即小样本研究缺失或选择性发表的风险较低,合并效应量可以较为放心地被解释。
六、Trim and Fill剪补法
表5报告了Trim and Fill剪补法的结果。分析显示,在当前漏斗图对称性良好的情况下,算法未识别出任何需要补充的缺失研究(补充0项),修正后的合并效应量与修正前完全一致(率=0.443, 95% CI: 0.403~0.483)。这一结果进一步印证了发表偏倚检验的结论,即当前Meta分析的合并效应量未受到发表偏倚的影响,总体估计结果是稳健的。
七、敏感性分析
表6报告了逐一剔除法的敏感性分析结果。剔除Hartman 2008后,合并率为0.441(95% CI: 0.385~0.498);剔除AHS 2018后为0.454(95% CI: 0.401~0.508);剔除ATBC 2011后为0.442(95% CI: 0.387~0.496);剔除CARET 2010后为0.428(95% CI: 0.398~0.459);剔除Weinstein 2015后为0.450(95% CI: 0.399~0.501)。五项剔除后的合并率波动范围为0.428~0.454,极差仅为0.026,波动幅度非常有限。所有剔除后的合并效应量均达到极显著水平(p<0.001),且I²均为0.000%,表明无论剔除哪项研究,合并结论均保持稳定。其中,剔除CARET 2010后合并率略降至0.428,剔除AHS 2018后略升至0.454,但变化幅度均不超过0.015,说明不存在对合并结论具有过度影响力的单项研究。
八、累积Meta分析
表7展示了按发表时间依次累积纳入研究的Meta分析结果。仅纳入Hartman 2008时,合并率为0.451(95% CI: 0.344~0.559),置信区间较宽,反映了单项研究估计精度有限。加入AHS 2018后,合并率降至0.430,但置信区间大幅拓宽至0.246~0.613,这与随机效应模型在研究数量较少时的估计特性有关。继续纳入ATBC 2011后,合并率稳定在0.434,置信区间收窄至0.384~0.483。纳入CARET 2010后,合并率升至0.450(95% CI: 0.399~0.501),区间进一步收窄。最终纳入Weinstein 2015后,合并率收敛至0.443(95% CI: 0.403~0.483),置信区间达到最窄宽度。从累积轨迹来看,合并率在0.430~0.454之间波动后逐步收敛,I²在纳入第三项研究后稳定在0.000%,表明证据基础已趋于一致和稳定。
九、图形结果解读
图1 森林图
图1展示了SPSSAU软件输出的森林图。从图中可以直观看出,五项研究的点估计值(方块)均集中在合并效应线(率=0.443)附近,各研究的95%置信区间(横线)存在较大重叠。合并效应量的菱形完全位于无效线右侧,且未跨越无效线,表明总体合并率具有统计学显著性。各研究方块的宽度(代表权重)与表2中的权重分配一致,AHS 2018的方块最宽,ATBC 2011最窄。
图2 漏斗图
图2展示了漏斗图。各研究点大致沿中心对称轴均匀分布,左右两侧近似对称,未出现明显的偏斜或缺失区域。这一视觉判断与Egger检验(p=0.765)和Begg检验(p=0.327)的统计结论一致,均表明不存在显著的发表偏倚。
图3 Trim and Fill剪补图
图3展示了Trim and Fill剪补法的结果图。图中原始研究以空心方块表示,算法未补充任何缺失研究(填充0项),合并效应量的菱形位置未发生改变。这进一步从图形角度证实了发表偏倚对当前合并结论无实质性影响。
图4 敏感性分析图
图4展示了逐一剔除法的敏感性分析图。从图中可以看出,剔除任一研究后的合并效应量(点估计及置信区间)与总体合并效应量高度重叠,波动范围极为有限,直观地展示了合并结论的稳健性。
图5 累积Meta分析图
图5展示了累积Meta分析的结果图。随着研究按时间顺序逐步纳入,合并效应量的点估计呈现先波动后收敛的趋势,置信区间逐步收窄,最终稳定在0.403~0.483范围内。这一收敛轨迹表明,当前领域的研究证据已逐步趋于一致。
十、结论
本研究基于随机效应模型对5项独立研究(合计n=555)进行了单个率Meta分析。主要发现如下:第一,总体合并率为0.443(95% CI: 0.403~0.483, z=30.739, p<0.001),估计精度较高。第二,异质性检验显示研究间一致性良好(Q=1.874, p=0.759, I²=0.000%),不存在显著的统计学异质性。第三,发表偏倚检验(Egger p=0.765; Begg p=0.327)和Trim and Fill剪补法均表明不存在可检测的发表偏倚。第四,敏感性分析显示逐一剔除任一研究后合并率在0.428~0.454之间波动,极差仅0.026,合并结论稳健。第五,累积Meta分析显示合并效应随研究纳入逐步收敛至0.443,证据基础趋于稳定。综上,当前Meta分析的合并效应方向明确、结论稳健,可为后续循证决策提供可靠的定量参考。