我将ggplot脚本包装成一个函数,如下所示。但是,运行这个会给我一些错误消息,例如
xy.coords(x,y,xlabel,ylabel,log)中的错误:
'x'是一个列表,但没有组件'x'和'y'
没有将这些脚本包装到函数中,它们只是工作正常。因此,如何编写这种函数,使用ggplot生成树形图并将其保存为pdf。
nicedendro<-function(inputdat, outputfile){
library(ggdendro)
library(ggplot2)
x <- read.table(inputdat, head=TRUE)
y <- 1-x
d <- as.dist(y,diag=FALSE,upper=FALSE)
hc <- hclust(d,"ave")
dhc <- as.dendrogram(hc)
ddata <- dendro_data(dhc,type="rectangle")
ddata$labels$text <- gsub("\\."," ",ddata$labels$text)
pdf(outputfile, width=30,height=35)
plot(ggplot(segment(ddata)) +
geom_segment(aes(x=x0,y=y0,xend=x1,yend=y1)) +
xlab(NULL) +
ylab(NULL) +
scale_x_discrete(limits=ddata$labels$text) +
opts(panel.grid.major = theme_blank()) +
opts(panel.grid.minor=theme_blank()) +
coord_flip())
dev.off()
}
答案 0 :(得分:5)
这是一个非常常见的问题:您必须使用print()
才能打印图表。
d <- function(){
g1 <- qplot(...)
print(g1)
}