我正试图从我使用R包metafor进行的固定效应荟萃分析中获得I ^ 2异质性统计量的置信区间。
似乎confint（）函数不适用于固定效果的荟萃分析？
赞赏任何指针！
答案 0 :(得分：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 package（v0.2.0）于2019-02-16创建。