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

Question

我正在尝试做几乎完全相同的事情： 从 R 中的输出中提取元分析估计和 CI

我有 metaprop 输出，我想提取置信区间估计值以放入数据帧，以便我可以将它们添加到我的图表中： 森林图，排序和汇总多个变量作为汇总数据的一行。

我使用了双反正弦变换，因此如果我想获得整个图的估计和 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)))

获取整个事件的信息。 即上面的代码给了我 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 没有丢失。”

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

Answer 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)

有关更多信息，请参阅https://www.researchgate.net/publication/325486099_How_to_Conduct_a_Meta-Analysis_of_Proportions_in_R_A_Comprehensive_Tutorial