2016-03-10 17 views
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

回答

0

這是我創建forestplot包的部分原因。它基於rmeta :: foresplot,因此您不必在函數調用中更改任何內容。