1
我使用rmeta包裝製作森林地塊。我的問題是節省製作的情節,因爲它沒有顯示整個數字。保存與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