2017-03-24 90 views
0

我創造了R的森林圖既然不能顯示我自己的數據,我使用示例代碼,並且我發現這裏 https://cran.r-project.org/web/packages/forestplot/vignettes/forestplot.html森林圖中的R標在tabletext矩陣

cochrane_from_rmeta <- 
structure(list(
mean = c(NA, NA, 0.578, 0.165, 0.246, 0.700, 0.348, 0.139, 1.017, NA, 0.531), 
lower = c(NA, NA, 0.372, 0.018, 0.072, 0.333, 0.083, 0.016, 0.365, NA, 0.386), 
upper = c(NA, NA, 0.898, 1.517, 0.833, 1.474, 1.455, 1.209, 2.831, NA, 0.731)), 
.Names = c("mean", "lower", "upper"), 
row.names = c(NA, -11L), 
class = "data.frame") 


tabletext<-cbind(
c("", "Study", "Auckland", "Block", 
"Doran", "Gamsu", "Morrison", "Papageorgiou", 
"Tauesch", NA, "Summary"), 
c("Deaths", "(steroid)", "36", "1", 
"4", "14", "3", "1", 
"8", NA, NA), 
c("Deaths", "(placebo)", "60", "5", 
"11", "20", "7", "7", 
"10", NA, NA), 
c("", "OR", "0.58", "0.16", 
"0.25", "0.70", "0.35", "0.14", 
"1.02", NA, "0.53"), 
c("",NA,NA,NA,NA,NA,NA,NA,NA,NA,"Heterogeniety I^2 = 20%"))[enter image description here][1] 


library(forestplot) 
forestplot(tabletext, 
     cochrane_from_rmeta,new_page = TRUE, 
     is.summary=c(TRUE,TRUE,rep(FALSE,8),TRUE), 
     clip=c(0.1,2.5), 
     xlog=TRUE, 
     col=fpColors(box="royalblue",line="darkblue", summary="royalblue")) 

我想上標數據Heterogeniety的最後一個標籤。我試過表達式(Heterogeniety〜I^2),但我得到以下錯誤信息:

cbind中的錯誤(c(「結果」,「死亡或BPD」,NA,NA,NA,「交互作用p值= .XX」:無法創建類型‘表達’

矩陣如何包括在我的森林圖上標

回答

0

R沒有允許您創建表達的矩陣的文檔中還提到了? forestplot

您還可以提供IDE矩陣雖然這不能被設計有表達

幸運的是,對於forestplot文件還提到,它可以在列表labeltext

該列表應被包裹在m×n個編號以類似於矩陣:列表(列表( 「rowname 1第1欄」, 「rowname 2欄1」),則列表( 「R1C2」,表達式(測試版) )

因此,您可能希望將tabletext轉換成嵌套列表如下:

tabletext <- list(
    list("", "Study", "Auckland", "Block", 
     "Doran", "Gamsu", "Morrison", "Papageorgiou", 
     "Tauesch", NA, "Summary"), 
    list("Deaths", "(steroid)", "36", "1", 
     "4", "14", "3", "1", 
     "8", NA, NA), 
    list("Deaths", "(placebo)", "60", "5", 
     "11", "20", "7", "7", 
     "10", NA, NA), 
    list("", "OR", "0.58", "0.16", 
     "0.25", "0.70", "0.35", "0.14", 
     "1.02", NA, "0.53"), 
    list("",NA,NA,NA,NA,NA,NA,NA,NA,NA, expression(Heterogeniety~I^2==20~'%'))) 


forestplot(tabletext, 
    cochrane_from_rmeta,new_page = TRUE, 
    is.summary=c(TRUE,TRUE,rep(FALSE,8),TRUE), 
    clip=c(0.1,2.5), 
    xlog=TRUE, 
    col=fpColors(box="royalblue",line="darkblue", summary="royalblue")) 

HTH

+0

這很好,但現在我有一個新問題。我無法使用graph.pos = 4(而之前這是可能的)在案例號和OR之間定位圖。我現在得到以下錯誤,如果我想改變forestplot的位置:forestplot中的錯誤(tabletext, cochrane_from_rmeta,graph.pos = 4,: 圖表位置必須介於1和2之間。您提供的值爲' 4'有沒有改變劇情位置的方法? – lenn