R:如何在包meta的metaprop中指定置信区间

Question

library(meta)
event <- c(81, 15, 0, 1)
n <- c(263, 148, 20, 29)
#
m1 <- metaprop(event, n, sm="PLOGIT", method.ci="SA")

我有兴趣在荟萃分析中结合比例。在上面的例子中，我有 4 个研究，每个研究报告一个比例。看来 metaprop计算 4 项研究中每一项的 CI。然而，由于 4 项研究已经报告了 CI(在原始论文中)，有没有办法将实际报告的 CI 纳入元分析计算中，而不是使用 metaprop自己计算？

如果其他软件包允许我指定 CI，我也愿意探索它们。

Answer 1

首先，请记住，您可能犯了一个巨大的错误。我认为没有理由不从 metaprop 中获取 CI功能。您应该从包中获取计算出的 CI，特别是在您使用随机效应模型时。平均个体效应和 CI 的计算方式与固定效应不同。无论如何，这里是:

您可以在以下位置找到置信下限:

> m1$lower
[1] 0.25219013 0.05273001 0.00000000 0.00000000

以及置信上限:

> m1$upper
[1] 0.3637795 0.1499727 0.0000000 0.1008922

您还可以访问 TE 中的单个平均效应:

> m1$TE
[1] -0.8095575 -2.1822989 -3.7135721 -3.3322045

然后，您可以在这些向量中手动添加您自己的计算值，它们将出现在您的森林图中。

假设您自己的较低值是 myLCL :

myLCL <- seq(-1, -.25, .25) # define lower confidence limits

m1$lower <- myLCL # assign them to the respective location
m1$lower
[1] -1.00 -0.75 -0.50 -0.25

你的置信上限是 myUCL :

myUCL <- seq(.25, 1, .25) # define upper confidence limits

m1$upper <- myUCL # assign them to the respective location
m1$upper
[1] 0.25 0.50 0.75 1.00