Metafor 森林图显示优势比的置信区间不正确

Question

我正在尝试绘制优势比和相应的 90% 置信区间，这是我之前使用包 brms 从比例优势模型中获得的。我制作了包含从模型中获得的优势比的向量（后验样本的中位数；x），以及下（ci.lb）和上（ci.ub）置信区间界限，并使用以下代码制作森林图：

x <- c(6.587028, 10.67589, 1.578881, 1.396755, 5.447785, 1.852427, 1.828179, 1.725313, 1.526206, 1.993046, 1.191804, 0.7648945)

ci.lb <- c(2.682959, 4.196124, 0.6783311, 0.547960, 2.011936, 0.7649611, 0.7492622, 0.6538183, 0.5299715, 0.8348141, 0.4814904, 0.2945799)

ci.ub <- c(17.843931, 29.200081, 3.7121095, 3.612463, 15.466248, 4.6075681, 4.3162345, 4.3404646, 4.3962568, 4.7397496, 2.9739579, 1.9560032)

forest(x, ci.lb, ci.ub, annotate=TRUE, showweights=FALSE, 
       header=headers, top=3, steps=5, 
       refline=0, digits=2L, xlab="Odds ratios", 
       slab=labels, efac=1, pch=15, col) 

但是，我的森林图上的置信区间与我的 ci.lb 和 ci.ub 向量中给出的置信区间不对应。这是我得到的情节：

谁能告诉我我做错了什么？

Answer 1

OP 显示的森林图来自 metafor 包，而不是 meta 包。后者加载metafor，运行forest(x, ci.lb, ci.ub)时使用的函数是metafor::forest.default()。如果您查看此函数的文档 (https://wviechtb.github.io/metafor/reference/forest.default.html)，您会看到前五个参数是 forest(x, vi, sei, ci.lb, ci.ub, ...)。因此，您将ci.lb 和ci.ub 向量传递给参数vi 和sei，这是没有意义的。所以使用forest(x, ci.lb=ci.lb, ci.ub=ci.ub) 来确保这些向量被传递给正确的参数。这里的总体教训是不要依赖参数的位置匹配。

Answer 2

您很可能正在使用 meta 包中的 forest。如果我加载那个库，我会得到和你一样的图：

library(meta)
forest(x, ci.lb, ci.ub)

这个函数需要一个meta 类的对象。对于你现在所拥有的，也许这会起作用：

library(forestplot)
studynames = paste0("study",1:length(x))

naming = list(
       "study"= studynames ,
        "odds ratio[95% CI]" = paste(round(x,digits=2),
        "[",
        round(ci.lb,digits=2),"-",
        round(ci.ub,digits=2),"]")
          )


forestplot(labeltext=naming,mean=x,lower=ci.lb,upper=ci.ub)