gpt4 book ai didi

linux - 在管道 awk 输出上打印循环输入文件名

转载 作者:太空狗 更新时间:2023-10-29 12:40:23 28 4
gpt4 key购买 nike

我希望在 awk 的每一行输出以下内容:

for fname in $(ls *.bam | cut -c -9 | uniq)
do
echo $fname
samtools depth -a "$fname".bam | awk '{sum+=$3} END {print "Average = ", sum/NR}' >> mean.depth
done

当前输出:

$ less mean.depth  
Average = 5.40089
Average = 4.68387

期望的输出:

$ less mean.depth  
FILENAME1 Average = 5.40089
FILENAME2 Average = 4.68387

这样输出的平均深度是可追踪的。我试过 FILENAME 和各种各样的东西但没有成功。谢谢!

最佳答案

解决方案

改编自this answer

awk -v var=定义shell变量$fname,然后用END {print var打印它

samtools depth  -a "$fname".bam | awk -v var=$fname '{sum+=$3} END {print var, " = ", sum/NR}' >> mean.depth

关于linux - 在管道 awk 输出上打印循环输入文件名,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/47091464/

28 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com