如何绘制嵌套列表中数据框的分组对比图？R语言栅格可视化问题

先明确数据结构假设 #

假设你的嵌套列表rast_df结构如下（每个子数据框仅包含栅格值一列）：

rast_df <- list(
  absorbance = list(
    farm1 = data.frame(value = rnorm(1000)),
    farm2 = data.frame(value = rnorm(800, mean = 1))
  ),
  resistance = list(
    farm1 = data.frame(value = rnorm(1200, mean = 2)),
    farm2 = data.frame(value = rnorm(900, mean = 3))
  )
)

步骤1：重构数据为长格式 #

要实现同一数据类型下两个农场的对比，首先需要把嵌套列表转换成包含数据类型、农场、栅格值的长数据框，这是ggplot做分组对比的基础：

library(tidyverse)

# 转换为长格式数据框
long_df <- rast_df %>%
  enframe(name = "data_type") %>%  # 提取外层的absorbance/resistance
  unnest(cols = value) %>%         # 展开内层的farm列表
  unnest(cols = value) %>%         # 展开每个农场的数据框
  rename(farm = name, value = value)  # 重命名列（如果你的列名不是value，自行调整）

步骤2：绘制对比小提琴图 #

分面对比（同一页面展示两种数据类型的农场对比） #

ggplot(long_df, aes(x = farm, y = value)) +
  geom_violin(na.rm = TRUE, scale = "count") +
  facet_wrap(~data_type, scales = "free_y") +  # 按数据类型分面，y轴自适应
  labs(x = "农场", y = "栅格值") +
  theme_bw()

单独绘制每个数据类型的对比图 #

如果需要分开保存或查看两种数据类型的图：

# 拆分数据并生成绘图列表
plot_list <- long_df %>%
  group_split(data_type) %>%
  map(function(sub_df) {
    ggplot(sub_df, aes(x = farm, y = value)) +
      geom_violin(na.rm = TRUE, scale = "count") +
      labs(x = "农场", y = "栅格值", title = unique(sub_df$data_type)) +
      theme_bw()
  })

# 查看absorbance的对比图
plot_list[[1]]
# 查看resistance的对比图
plot_list[[2]]

原代码的问题根源 #

1. aes参数写法错误：aes(x='', y=df[,1])不符合ggplot的语法逻辑——ggplot会自动从传入的df中读取列，应该直接使用列名（比如y=value），而非df[,1]这种外部引用方式，这会导致ggplot无法正确识别数据来源，触发类型错误。
2. rapply使用不当：rapply的how="unlist"和deflt=NA_integer_试图将ggplot对象强制转换为整数类型，而ggplot对象是复杂的S3对象，完全不兼容这种转换，这是报错的直接原因。
3. 缺少分组标识：原代码没有保留数据类型、农场的分组信息，无法实现同一数据类型下两个农场的对比，仅能单独绘制每个农场的孤立小提琴图，不符合需求。

内容的提问来源于stack exchange，提问作者Enrico Gabrielli