Meta和Metafor R软件包-

编程问答 2022-08-14

问题描述

我目前正在使用“ Metafor”软件包在R中进行元分析。在进行研究时，我遇到了另一套用于R中的荟萃分析的软件包，即“元”。我更喜欢在后者的包中（设计方面）创建的森林图，但是不幸的是，某些数据与我使用Metafor创建的图不同。

具体来说，只有I ^ 2和汇总估算值的数据不同。

Meta_1 <- rma(yi=yi,vi=vi,measure="SMD",method="ML",slab=Citation,data=dat)
forest(Meta_1)

Meta_2 <- Metagen(yi,vi^.5,data = dat,studlab = paste(Citation),comb.fixed = FALSE,comb.random = TRUE,hakn = TRUE,method.tau = "ML",sm = "SMD")
forest(Meta_2)

有人知道为什么会出现这些差异吗？

解决方法

因此，我能够获得预测间隔以匹配各个函数，但不能匹配I ^ 2值（即使相差仅2％）。与其他软件包相比，一个软件包可能在做一些统计上的修正，或者与RE / FE类型的建模方法有关。

无论如何，我希望这段代码能为您指明正确的方向。为了使配置项匹配，还必须使用method.tau.ci中的参数metagen()。

library(meta)
library(metafor)
study<- c(1:10) 
yi<- c( -0.48965031,0.64970214,0.11201680,0.07945655,-0.70874645 -0.54922759,0.66768916,-0.45523574 ) 
vi <- c(0.10299697,0.14036855,0.05137812,0.03255550,0.34913525,0.34971466,0.07539957,0.08428983) 
dat <- cbind(study,yi,vi) 
dat <- as.data.frame(dat)

meta_1 <- rma(yi=dat$yi,vi=dat$vi,measure="SMD",method="REML",slab=paste(study),data=dat) 
forest(meta_1) 

meta_2 <- meta::metagen(TE =dat$yi,seTE = dat$vi^.5,method.tau = 'REML',method.tau.ci = 'BJ',comb.random = TRUE,comb.fixed = TRUE,sm = 'SMD') 

forest(meta_2)
metaforrr