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,我也愿意探索其他包.
首先，请记住，您可能会犯一个巨大的错误。我看不出有什么理由不从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