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自己进行计算？

如果允许我指定配置项，我也愿意探索其他软件包。

Answer 1

首先，请记住，您可能犯了一个巨大的错误。我没有理由不从metaprop函数中获取配置项。您应该从包装中获取计算出的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

R：如何在包元的元属性中指定置信区间

1 个答案: