r - Metafor：如何在固定效应荟萃分析中获得 I^2 的置信区间

Question

我试图从我使用 R 包 metafor 运行的固定效应荟萃分析中获得 I^2 异质性统计的置信区间。

似乎 confint() 函数不适用于固定效应荟萃分析？

感谢任何指针！

Answer 1

固定效应荟萃分析的假设是效应中没有异质性。观察到的效果的所有差异都假定是由于抽样方差。由于 I^2 是可归因于异质性而非抽样方差的观察方差的比例，因此当您使用固定效应模型时，metafor 不会估计 I^2 值（有关使用固定效应和随机效应的样本输出，请参见下文）。这就是您无法获得 I^2 值周围的置信区间的原因。

library(metafor)
#> Warning: package 'metafor' was built under R version 3.5.2
#> Loading required package: Matrix
#> Loading 'metafor' package (version 2.0-0). For an overview 
#> and introduction to the package please type: help(metafor).

rma(yi = yi, vi = vi, data = dat.bangertdrowns2004, method = "FE")
#> 
#> Fixed-Effects Model (k = 48)
#> 
#> Test for Heterogeneity: 
#> Q(df = 47) = 107.1061, p-val < .0001
#> 
#> Model Results:
#> 
#> estimate      se    zval    pval   ci.lb   ci.ub     
#>   0.1656  0.0269  6.1499  <.0001  0.1128  0.2184  ***
#> 
#> ---
#> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

rma(yi = yi, vi = vi, data = dat.bangertdrowns2004)
#> 
#> Random-Effects Model (k = 48; tau^2 estimator: REML)
#> 
#> tau^2 (estimated amount of total heterogeneity): 0.0499 (SE = 0.0197)
#> tau (square root of estimated tau^2 value):      0.2235
#> I^2 (total heterogeneity / total variability):   58.37%
#> H^2 (total variability / sampling variability):  2.40
#> 
#> Test for Heterogeneity: 
#> Q(df = 47) = 107.1061, p-val < .0001
#> 
#> Model Results:
#> 
#> estimate      se    zval    pval   ci.lb   ci.ub     
#>   0.2219  0.0460  4.8209  <.0001  0.1317  0.3122  ***
#> 
#> ---
#> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

由reprex 包(v0.2.0) 于 2019 年 2 月 16 日创建。