Question

我想知道是否可以访问（以某种形式）-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

Answer 1

怎么样：

> 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

然后yi，ci.lb和ci.ub提供相同的信息。

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

1 个答案: