在Meta分析中从metafor森林图中访问/保存信息

时间:2017-01-25 18:36:40

标签: r

我想知道是否可以访问(以某种形式)-metafor-package中-forest-命令中显示的信息。

我正在检查/验证结果,我希望生成值的输出。到目前为止,计算都是检查,但我想让它们可用于打印,保存等,而不必手动输入。

示例代码如下:

es <- read.table(header=TRUE, text = "
b se_b
0.083 0.011
0.114 0.011
0.081 0.013
0.527 0.017
" )


library(metafor)

es.est <- rma(yi=b, sei=se_b, dat=es, method="DL")


studies <- as.vector(  c("Larry (2011)" , "Curly (2011)", "Moe (2015)" , "Shemp (2010)" )   )


forest(es.est , transf=exp  , slab = studies , refline = 1 ,  xlim=c(0,3), at = c(1, 1.5, 2, 2.5, 3, 3.5, 4) , showweights=TRUE)

我想访问图表右侧显示的值(每项研究的效果大小和c.i.,以及总体估计值和c.i.)。

非常感谢,

-Jon

1 个答案:

答案 0 :(得分:0)

怎么样:

> summary(escalc(measure="GEN", yi=b, sei=se_b, data=es), transf=exp)
      b  se_b     yi     vi    sei      zi  ci.lb  ci.ub
1 0.083 0.011 1.0865 0.0001 0.0110  7.5455 1.0634 1.1102
2 0.114 0.011 1.1208 0.0001 0.0110 10.3636 1.0968 1.1452
3 0.081 0.013 1.0844 0.0002 0.0130  6.2308 1.0571 1.1124
4 0.527 0.017 1.6938 0.0003 0.0170 31.0000 1.6383 1.7512

然后yici.lbci.ub提供相同的信息。