如何将多个BED文件的指定列合并为单个表格文件？

先排查bedops "Unable to find file"报错 #

• 确认文件路径：用pwd查看当前工作目录，ls检查目标文件是否存在；若文件不在当前目录，使用绝对路径（如/home/user/data/1.bed）或先通过cd切换到文件所在目录
• 核对文件名大小写：Linux/macOS系统区分大小写，确保命令中的文件名（如1.bed）和实际文件名完全一致
• 处理特殊文件名：如果文件名包含空格或特殊字符，用双引号包裹（如"sample 1.bed"）

正确合并BED文件的方法 #

方法1：用bedops + bedmap（推荐，适合大规模数据） #

1. 生成统一的基础坐标文件
   提取所有BED文件的前3列，去重并按BED格式排序（bedops要求输入为排序后的文件）：

   cat *.bed | cut -f1-3 | sort-bed - | uniq > all_coords.bed
   

2. 生成合并表格
   先写入表头，再用bedmap映射每个文件的第4列数值：

   # 写入表头
   echo -e "chrom\tstart\tend\t$(ls *.bed | tr '\n' '\t')" | sed 's/\t$//' > merged_table.txt
   # 映射数值并追加到表格
   bedmap --echo --echo-map-id-uniq --delim '\t' all_coords.bed $(printf "%s " *.bed | sed 's/ / --map-id-file /g') >> merged_table.txt
   

   说明：--echo输出基础坐标，--echo-map-id-uniq提取每个文件的第4列值，--delim '\t'指定制表符作为分隔符，--map-id-file指定要提取数值的目标BED文件。

方法2：用awk脚本（无需额外工具，灵活通用） #

如果没有安装bedops/bedtools，可直接用awk处理（确保所有BED文件已按chrom/start/end排序；未排序的话先执行for f in *.bed; do sort-bed $f > sorted_$f; done，再用sorted_*.bed作为输入）：

awk '
BEGIN {
    OFS="\t"
    for (i=1; i<ARGC; i++) {
        headers = headers (headers == "" ? "" : OFS) ARGV[i]
    }
    print "chrom", "start", "end", headers
    ARGC = 1
}
FNR == 1 { file_idx++ }
{
    key = $1 OFS $2 OFS $3
    values[key][file_idx] = $4
    keys[key] = 1
}
END {
    for (key in keys) {
        printf "%s", key
        for (i=1; i<=file_idx; i++) {
            printf "%s%s", OFS, (values[key][i] == "" ? "." : values[key][i])
        }
        print ""
    }
}
' *.bed > merged_table.txt

说明：脚本自动读取所有输入BED文件，以文件名作为表头列名，将同一坐标的数值对应填充，缺失值用.表示。

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