Shiny应用中system()调用vsearch生成空文件的解决方法

以下是针对你遇到的问题的具体排查和修复步骤：

1. 修复文件路径问题（核心原因） #

Shiny运行时的工作目录与你单独执行R脚本的目录不同，且你调用input$seq$name仅获取了文件名，没有完整路径，导致vsearch找不到输入文件，最终生成空输出。

修复代码：
修改pl函数，使用上传文件的完整临时路径（input$seq$datapath），并为所有输出文件指定绝对路径：

pl <- function(file_path) {
  # 指定输出目录（可换成你想要的固定目录）
  output_dir <- getwd()
  # 生成输出文件的绝对路径
  sort_out <- file.path(output_dir, "CLADE_sort.fa")
  centroids_out <- file.path(output_dir, "CLADE.clustered.fa")
  uc_out <- file.path(output_dir, "CLADE.cluster")
  
  # 构造带绝对路径的vsearch命令
  sort_cmd <- paste("vsearch --sortbylength", file_path, 
                    "--output", sort_out, 
                    "--minseqlength 500 -notrunclabels")
  system(sort_cmd, intern = TRUE, wait = TRUE)
  
  cluster_cmd <- paste("vsearch --cluster_smallmem", sort_out, 
                       "--id 0.97 --centroids", centroids_out, 
                       "-uc", uc_out)
  system(cluster_cmd, intern = TRUE, wait = TRUE)
}

然后修改触发逻辑，用input$seq$datapath替代input$seq$name，并绑定按钮事件而非自动触发：

# 替换原来的output$pipeline，改用按钮触发
observeEvent(input$seqbut, {
  req(input$seq)  # 确保文件已上传
  pl(input$seq$datapath)
  showNotification("vsearch处理完成！", type = "success")
})

2. 捕获vsearch执行的错误信息 #

添加错误捕获，查看vsearch实际运行时的报错，快速定位问题：

# 修改system调用，捕获输出和错误信息
sort_result <- system(sort_cmd, intern = TRUE, wait = TRUE, ignore.stderr = FALSE)
cluster_result <- system(cluster_cmd, intern = TRUE, wait = TRUE, ignore.stderr = FALSE)
# 在Shiny控制台打印结果，便于调试
print(sort_result)
print(cluster_result)

3. 确保vsearch可被Shiny进程找到 #

如果Shiny的PATH环境变量未包含vsearch路径，直接使用vsearch的绝对路径调用：

# 先在终端执行`which vsearch`获取绝对路径，替换下方路径
vsearch_path <- "/usr/local/bin/vsearch"
sort_cmd <- paste(vsearch_path, "--sortbylength", file_path, ...)

4. 检查文件读写权限 #

确保Shiny进程对输出目录有读写权限，可指定一个明确的权限目录：

output_dir <- path.expand("~/Shiny_vsearch_results")
# 创建目录（如果不存在）
dir.create(output_dir, recursive = TRUE, showWarnings = FALSE)
# 后续输出文件都放到这个目录
sort_out <- file.path(output_dir, "CLADE_sort.fa")

内容的提问来源于stack exchange，提问作者Marta Noguera Cardona