我用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