从子组的输出中提取元分析估计和置信区间答案

【问题标题】：Extract Meta-Analysis Estimate and Confidence Interval from output for subgroups从子组的输出中提取元分析估计和置信区间
【发布时间】：2020-09-22 04:02:43
【问题描述】：

我正在尝试做这件事：Extracting meta analysis estimate and CI from output in R

我有 metaprop 输出，我想提取置信区间估计值以放入数据框中，以便将它们添加到我的图表中：Forest Plot, ordering and summarizing multiple variables 作为汇总数据的一行。

我使用了双反正弦变换，因此如果我想获得整个绘图的估计值和 CI，那么第一个答案中列出的代码可以完美运行，但我有子组。那么如何提取下图中绿色的数字，而不是红色的数字呢？

编辑：更新更多信息，因此根据上面引用的 SE 问题，您可以：

random.est2 <- c(pes.summary$TE.random,pes.summary$lower.random,pes.summary$upper.random)
unlist(lapply(random.est2, meta:::backtransf,  sm="PFT",  n=1/mean(1/pes.summary$n)))

为了得到整个事情的信息。 IE。上面的代码给了我 0.988492 [0.976040, 0.997163]。这些数字与照片中红色方块中的数字完全一致。因此，它是您为此提取的“TE.Random”对象。嗯...根据 metaprop 的文档，“Te.Random.w”应该是子组的估计治疗效果。所以我想我会为子组执行与上面相同的代码并尝试：

random.est2 <- c(pes.summary$TE.random.w,pes.summary$lower.random.w,pes.summary$upper.random.w)
unlist(lapply(random.est2, meta:::backtransf,  sm="PFT",  n=1/mean(1/pes.summary$n)))

但这给了我大约 0.01 的数字。所以我假设我可能错误地反向转换它们？

最后，似乎有一个 n.harmonic.mean.w，它是：“子组中观察数的谐波平均值（用于 Freeman-Tukey 双反正弦变换的反变换） - 如果不缺少 byvar。”

这听起来更像是我应该使用的，但我也无法让它发挥作用。

【问题讨论】：

标签： r object

【解决方案1】：

试试这个。

library(metafor)

Estimate<-transf.ipft(pes.summary$TE.random.w,ni=pes.summary$n.harmonic.mean.w)
Upper<-transf.ipft(pes.summary$lower.random.w,ni=pes.summary$n.harmonic.mean.w)
Lower<-transf.ipft(pes.summary$upper.random.w,ni=pes.summary$n.harmonic.mean.w)


new.df<-data.frame(Estimate,Upper,Lower)

【讨论】：