好吧,我有几个文件,但我无法使用 awk 将它们连接起来。以下是三个示例文件:
文件A
2 2 31771 15 5327353 4.73E-04 1 1
2 2 40642 23 27428869 -8.29E-04 1 1
2 2 21517 7 89002990 -2.41E-04 1 1
2 2 33787 16 44955568 2.86E-05 1 1
Run Code Online (Sandbox Code Playgroud)
文件B
2 2 31771 15 5327353 5.07E-04 1 1
2 2 40642 23 27428869 5.45E-04 1 1
2 2 21517 7 89002990 1.85E-04 1 1
2 2 33787 16 44955568 3.73E-04 1 1
Run Code Online (Sandbox Code Playgroud)
文件C
2 2 31771 15 5327353 4.28E-04 1 1
2 2 40642 23 27428869 -7.55E-04 1 1
2 2 21517 7 89002990 -2.01E-04 1 1
2 2 33787 16 44955568 3.09E-05 1 1
Run Code Online (Sandbox Code Playgroud)
每个文件有 8 列,但我不需要第 1、2、7 和 8 列。第 3、4 和 5 列是所有文件共有的,这些可能是用于连接文件的参考列,第 6 列将是每个文件的不同信息。一般来说,我的最终文件如下所示:
最终文件
31771 15 5327353 4.73E-04 5.07E-04 4.28E-04
40642 23 27428869 -8.29E-04 5.45E-04 -7.55E-04
21517 7 89002990 -2.41E-04 1.85E-04 -2.01E-04
33787 16 44955568 2.86E-05 3.73E-04 3.09E-05
Run Code Online (Sandbox Code Playgroud)
我尝试了以下命令:
awk 'NR==FNR{a[$3]=$6;next}{print $3"\t"$4"\t"$5"\t"$6"\t"a[$3]}' FileA FileB FileC > Finalfile
Run Code Online (Sandbox Code Playgroud)
但不幸的是它只适用于两个文件,我得到这样的结果:
2 2 31771 15 5327353 4.73E-04 5.07E-04
2 2 40642 23 27428869 -8.29E-04 5.45E-04
2 2 21517 7 89002990 -2.41E-04 1.85E-04
2 2 33787 16 44955568 2.86E-05 3.73E-04
Run Code Online (Sandbox Code Playgroud)
有人可以帮忙吗?请记住,有多个文件 (10),而不仅仅是三个。非常感谢!
尝试一下:
awk '{a[$3FS$4FS$5]=a[$3FS$4FS$5]""$6FS}END{for (i in a){print i, a[i]}}' file*
Run Code Online (Sandbox Code Playgroud)
更干净的版本(感谢@james-brown):
awk '{ a[$3 OFS $4 OFS $5 FS] = a[$3 OFS $4 OFS $5 FS] ( a[$3 OFS $4 OFS $5 FS] == "" ? "" : OFS) $6 }
END{ for (i in a){print i,a[i]} }' OFS="\t" file*
Run Code Online (Sandbox Code Playgroud)
输出
33787 16 44955568 2.86E-05 3.73E-04 3.09E-05
21517 7 89002990 -2.41E-04 1.85E-04 -2.01E-04
40642 23 27428869 -8.29E-04 5.45E-04 -7.55E-04
31771 15 5327353 4.73E-04 5.07E-04 4.28E-04
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
2712 次 |
| 最近记录: |