2016-10-01 25 views
3

私はRに慣れていないので、素朴な質問を許してください。私はRブロガーに見て適応したコードを使って森林プロットを生成しました。私は脚注を追加したいと思いますが、どうすればよいか分かりません。私はさまざまなコードを見てきましたが、どのコードが適切かわかりません。どんなアドバイスも大歓迎です!ありがとう!forestplot Rパッケージを使用している場合、どのように脚注を追加しますか?

フォレストプロットコードは以下である:

library(ggplot2) 
library(grid) 
library(gridExtra) 
library(xtable) 
library(plyr)   
library(MASS)  
library(reshape2) 

d1 <- read.csv("/view/name_view/vob/folder/subfolder/forestTAF.csv", stringsAsFactors=FALSE) 

#Format data so column labels and columns match required graphical output 
subgps <- c(4,5,8,9,12,13,16,17,20,21,24,25,28,29) 

#Combine the rest of the columns in the table 
tabletext <- cbind(c("Change due to", d1$Variable), 
       c("PK Measures", d1$PK.Measures), 
       c("Recommendation", d1$Recommendation)) 

#Load forestplot package 
library(forestplot) 

#Create forest plot 
png(filename="forestTAF.png",width=1600, height=640) 
forestplot(labeltext=tabletext, graph.pos=3, align=c("l", "c", "l"), 
     mean=c(NA,d1$Point.Estimate), 
     lower=c(NA,d1$Low), upper=c(NA,d1$High), 
     title="Drugname \nFold Change and 90% CI", 
     xlab="Change Relative to Reference (i.e., without interaction)", 
      txt_gp=fpTxtGp(label=gpar(cex=1.25), 
          ticks=gpar(cex=1.1), 
          xlab=gpar(cex = 1.2), 
          title=gpar(cex = 1.2)), 
     col=fpColors(box="black", lines="black"), 
     zero=1, cex=0.9, lineheight = "auto", boxsize=0.5, colgap=unit(6,"mm"), 
     lwd.ci=2, ci.vertices=TRUE, ci.vertices.height = 0.4) 

答えて

2

grid::gridtextは一つの選択肢grid::unit"npc"が正規0,1そう0.05は全体の5%であろう座標

library('forestplot') 

## ?forestplot 

row_names <- list(list("test = 1", expression(test >= 2))) 
test_data <- data.frame(coef=c(1.59, 1.24), 
         low=c(1.4, 0.78), 
         high=c(1.8, 1.55)) 
forestplot(row_names, test_data$coef, test_data$low, test_data$high, 
      zero = 1, cex = 2, lineheight = "auto", xlab = "Lab axis txt") 

grid::grid.textあります画像

x <- unit(.05, 'npc') 
y <- unit(.05, 'npc') 
grid.text('A footnote', x, y, gp = gpar(fontsize=10, font = 3)) 

または、grid::grid.locatorを使用して、希望の位置をポイントしてクリックすることができます。

g <- grid.locator('npc') 
grid.text('A footnote', g$x, g$y, gp = gpar(fontsize=10, font = 3)) 

enter image description here

+0

、そんなに@rawrをありがとうございます!最初のコードはトリックでした:-) – Christina

関連する問題