R语言生物分析：修改物种分布图函数以支持批量处理

你的批量处理思路是可行的，以下是针对代码问题的修正方案及说明：

核心问题排查与修正点 #

1. 统一输入类型为字符向量
   如果传入的是数据框的因子列（factor），循环迭代时会出现匹配错误，需先转换为字符向量。

2. 替换subset为更稳定的筛选函数
   subset在循环场景下对特殊字符的兼容性较差，改用dplyr::filter能提升数据筛选的稳定性。

3. 增加空数据判断
   若物种在birds_maps中无匹配数据，直接生成图表会报错，需跳过这类物种并给出提示。

4. 优化文件路径生成
   使用file.path替代paste0，自动适配不同系统的路径分隔符，避免路径错误。

修正后的完整代码 #

library(dplyr)
library(sf)
library(ggplot2)
library(stringr)

species.distrib.map <- function(species_names, distribution_color = "#f16d3d", file_type = ".pdf") {
  # 转换输入为字符向量，兼容因子列输入
  species_names <- as.character(species_names)
  
  for (species_name in species_names) {
    # 用dplyr::filter替代subset，提升筛选稳定性
    spec_pol <- birds_maps %>% filter(Species == species_name)
    
    # 空数据判断，跳过无匹配的物种
    if(nrow(spec_pol) == 0){
      message(paste("Warning: No data found for species", species_name, ", skipping..."))
      next
    }
    
    spec_sf <- st_as_sf(spec_pol)
    
    spec_map <- ggplot() +
      geom_sf(data = Depart_sf %>%
                filter(!DEPTO == "SAN ANDRES,PROVIDENCIA Y SANTA CATALINA"),
              fill = "transparent", color = "darkgray", linewidth = 0.1) +
      geom_sf(data = spec_sf, fill = distribution_color, alpha = 0.5, color = NA) +
      geom_sf(data = Colb_sf, fill = "transparent", color = "darkgray", linewidth = 0.3) +
      geom_point(data= CIAT, aes(x = Longitud, y = Latitud), color = "#265f66", size = 1) +
      geom_rect(data = rect_df_colb, aes(xmin = xmin, xmax = xmax, ymin = ymin, ymax = ymax), color = "black", fill = NA, linewidth= 0.57) +
      theme_void()
    
    # 用file.path生成兼容跨系统的路径
    filename <- file.path("D:/R/birdmaps/Species maps", paste0(species_name, file_type))
    message(paste("Saving map for species:", species_name, "to file:", filename))
    
    ggsave(filename = filename, plot = spec_map, width = 4.91, height = 8, units = "in", dpi = 700)
    print(spec_map)
  }
}

使用示例 #

• 传入字符列表：

  species_list <- c("Hirundo rustica", "Passer domesticus")
  species.distrib.map(species_list)
  

• 传入数据框的物种列：

  # 假设your_data包含Species列
  species.distrib.map(your_data$Species)
  

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