【发布时间】:2023-03-15 05:29:01
【问题描述】:
如何在回收数据时制作多个相同数据但因不同因素(列)而不同颜色的图?这是gridExtra 与cowplot 所做的不同吗?
目标: 我的目标是直观地比较有效聚类相同数据的不同结果。 我目前认为,直观地比较 2-4 种聚类算法的最简单方法是将它们彼此相邻绘制。
因此,如何将相同的数据并排绘制不同的颜色?
挑战/规范:性能非常重要。我要制作大约 30,000 张图表,每张图表有 450 - 480 分。数据被“回收”是至关重要的。
我可以使用包cowplot 和gridExtra 并排绘制它们。我今天才开始使用gridExtra,但它似乎可以回收数据,并且就我的目的而言比cowplot 更好。 更新: u/eipi10 证明如果我在绘图前收集列,facet_wrap 可以工作。
设置
#Packages
library(ggplot2)
library(cowplot)
library(gridExtra)
library(pryr) #memory profile
#Data creation
x.points <- c(1, 1, 1, 3, 3, 3, 5, 5, 5)
y.points <- c(1, 3, 5, 1, 3, 5, 1, 3, 5)
cl_vert <- c("A", "A", "A", "B", "B", "B", "C", "C", "C")
cl_hoz <- c("A", "B", "C", "A", "B", "C", "A", "B", "C")
cl_cent <- c("A","A","A","A", "B", "A","A","A","A")
df <- data.frame(x.points, y.points, cl_vert, cl_hoz, cl_cent)
绘制图表
#Graph function and individual plots
graph <- function(data = df, Title = "", color.by, legend.position = "none"){
ggplot(data, aes(x = `x.points`, y = `y.points`)) +
geom_point(aes(color = as.factor(color.by))) + scale_color_brewer(palette = "Set1") +
labs(subtitle = Title, x = "log(X)", y = "log(Y)", color = "Color" ) +
theme_bw() + theme(legend.position = legend.position)
}
g1 <- graph(Title = "Vertical", color.by = cl_vert)
g2 <- graph(Title = "Horizontal", color.by = cl_hoz)
g3 <- graph(Title = "Center", color.by = cl_cent)
#Cowplot
legend <- get_legend(graph(color.by = cl_vert, legend.position = "right")) #Not a memory waste
plot <- plot_grid(g1, g2, g3, labels = c("A", "B", "C"))
title <- ggdraw() + draw_label(paste0("Data Ex ", "1"), fontface = 'bold')
plot2 <- plot_grid(title, plot, ncol=1, rel_heights=c(0.1, 1)) # rel_heights values control title margins
plot3 <- plot_grid(plot2, legend, rel_widths = c(1, 0.3))
plot3
#gridExtra
plot_grid.ex <- grid.arrange(g1, g2, g3, ncol = 2, top = paste0("Data Ex ", "1"))
plot_grid.ex
pryr 的内存使用情况
#Comparison
object_size(plot_grid) #315 kB
object_size(plot3) #1.45 MB
#Individual objects
object_size(g1) #756 kB
object_size(g2) #756 kB
object_size(g3) #756 kB
object_size(g1, g2, g3) #888 kB
object_size(legend) #43.6 kB
其他问题:
写完这个问题并提供了样本数据,我才想起gridExtra,试了一下,似乎比它的组件图的组合数据占用内存少。我认为 g1、g2 和 g3 共享相同的数据,但着色分配除外,这就是为什么各个组件与总对象大小之间存在大约 130 kB 差异的原因。为什么 plot_grid 占用的空间比这还少? ls.str(plot_grid) 似乎没有显示 g1、g2 和 g3 的任何合并。我最好的选择是使用lineprof() 并逐行比较吗?
我浏览/阅读/咨询过的来源:
- http://adv-r.had.co.nz/memory.html #不完全理解
- Add a common Legend for combined ggplots #稍后修复 gridExtra
请多多包涵,因为我是一名新程序员(12 月才真正开始编写脚本);我还不了解所有技术细节,但我想了解。
【问题讨论】:
-
你打算如何查看 30,000 个地块?也许有更好的方法来做你想做的事?
-
我不会一次处理所有 30,000 组数据。我也不打算保存图表本身,而是保存用于制作它们的数据。实际上,我将相同格式的数据分布在约 500 个文件夹中(按簇数和 HDBSCAN 给出的 p_value 分类)。我一次加载一个 csv 文件。那个 csv 可以有 2 组数据到 200。我使用 tidyverse
group_by()和nest()来形成一个嵌套列表,然后运行 lapply 来制作所有的图表。然后我一次调用一个图表(我有一个闪亮的应用程序来导航目录并显示正确的图表)
标签: r performance memory ggplot2 gridextra