[英]Access / save information from metafor forest plot in meta-analysis
I'm wondering if it's possible to access (in some form) the information that is presented in the -forest- command in the -metafor- package.我想知道是否可以访问(以某种形式)在 -metafor- 包中的 -forest- 命令中显示的信息。
I am checking / verifying results, and I'd like to have the output of values produced.我正在检查/验证结果,并且我希望生成值的输出。 Thus far, the calculations all check, but I'd like to have them available for printing, saving, etc. instead of having to type them out by hand.
到目前为止,所有计算都经过检查,但我希望它们可用于打印、保存等,而不必手动输入它们。
Sample code is below :示例代码如下:
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)
I'd like to access the values (effect size and ci for each study, as well as the overall estimate, and ci) that are presented on the right of the graphic.我想访问图形右侧显示的值(每项研究的效应大小和 ci,以及总体估计和 ci)。
Thanks so much,非常感谢,
-Jon -乔恩
How about:怎么样:
> 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
Then yi
, ci.lb
, and ci.ub
provides the same info.然后
yi
、 ci.lb
和ci.ub
提供相同的信息。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.