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)
謝謝你這麼多,@rawr!第一個代碼做了訣竅:-) – Christina