我正在尝试将ggplot2地块添加到列表中,以便可以使用ggpubr的ggarrange函数将它们组织在多个页面上。
由于有数百个图,因此我正在使用一个函数来生成和保存图,但是我无法获得将图返回到环境或将名称写入列表的功能。
我敢肯定,这很简单,我很想念但找不到它。
我正在使用的绘图功能是:
histFacet.plot <- function(x, results, info, ...) {
md<- names(x) %in% c("rn","Taxa","year","rep","block","column",
"range", "entity_id")
traits <- names(x[ , !md])
for (i in traits) {
i <-ggplot(data = x, aes_string(x = i)) +
geom_histogram(colour="black", fill="white") +
#facet_grid(x$year ~ .) +
theme_bw() +
xlab(paste0(i)) +
ylab("Frequency") +
theme(panel.grid.major = element_blank()) +
theme(panel.grid.minor = element_blank()) +
theme(axis.text = element_text(size = 15)) +
theme(axis.title = element_text(size = 15)) +
theme(strip.text = element_text(size = 15))
#ggsave(paste0(i,"_",info,".pdf"),path=paste(results, sep=''))
plotList<- list(plotList, list(i))
print(i)
}
return(i)
}
histFacet.plot(pd,'~/Dropbox/Research_Poland_Lab/AM Panel/Figures/Hist/',
"_raw_2018")
答案 0 :(得分:0)
您的大问题是您return(i)
而不是return(plotList)
。在迭代器中使用i
在for
循环内重新分配i
很奇怪。尤其是当情节使用i
作为字符串时,我会尝试这样做:
histFacet.plot <- function(x, results, info, ...) {
md <- names(x) %in% c("rn","Taxa","year","rep","block","column",
"range", "entity_id")
traits <- names(x[ , !md])
plotList = list()
for (i in traits) {
thisPlot <- ggplot(data = x, aes_string(x = i)) +
geom_histogram(colour="black", fill="white") +
#facet_grid(x$year ~ .) +
theme_bw() +
xlab(i) +
ylab("Frequency") +
theme(panel.grid.major = element_blank()) +
theme(panel.grid.minor = element_blank()) +
theme(axis.text = element_text(size = 15)) +
theme(axis.title = element_text(size = 15)) +
theme(strip.text = element_text(size = 15))
plotList[[i]] = thisPlot
print(i)
}
return(plotList)
}
由于您没有共享示例数据,因此未经测试,当然。如果您提供了一个小的,可重现的示例数据集,则很乐意进行测试/调试。
我不确定是否要print(i)
将绘图打印到图形设备上(如果是,请更改为print(thisPlot)
)还是将当前特征打印到控制台上以更新进度循环(如果是这样,请更改为message(i)
以使其易于禁用)。
其他一些注意事项:如果您是多方面的,请使用year ~ .
作为公式,而不要使用x$year ~ .
。如果i
已经成为书迷,则paste0(i)
与i
(在您的xlab
中)相同。