用rmeta保存森林地块

时间:2016-03-10 08:48:40

标签: r forestplot

我用rmeta包制作森林情节。我的问题是保存生产的情节,因为它没有显示整个数字。

我的示例数据:

Results
metabolite   expb   lower95 upper95  pvalue   se
a            1,02    0,98    1,04     0,05    0.01
b            1,03    0,99    1,05     0,06    0.02
c            1,05    1,02    1,10     0,08    0.03
...
number of metabolites is 70

我按如下方式制作森林情节

library(rmeta)
d <- meta.summaries(results$expb, results$se, names = results$metabolite, method = c("fixed"), data = results)

tabletext <- cbind(c("95%", "Lower", results$lower95, NA, NA), c("C.I.", "Upper", results$upper95, NA, NA), c("", "P-value", results$pvalue, NA, round(d$summary, 2)))

m <- c(NA, NA, results$expb, NA, d$summary)
m <- round(m, digits = 2)
l <- c(NA, NA, results$lower95, NA, d$summary - d$se.summary * 1.96)
l <- round(l, digits = 2)
u <- c(NA, NA, results$upper95, NA, d$summary + d$se.summary * 1.96)
u <- round(u, digits = 2)

forestplot(tabletext, m, l, u, zero = 1, summary = c(TRUE, TRUE, rep(FALSE, 
8), TRUE), boxsize = 0.2,cex.axis=0.2)

然而,当我保存时,它会削减产生的情节。你能帮我解决一个问题吗?谢谢你的任何提示。 PS。您可以在以下链接中找到示例图。 Example figure

1 个答案:

答案 0 :(得分:0)

这是我创建forestplot包的部分原因。它基于rmeta :: foresplot,因此您不必在函数调用中更改任何内容。