你能帮我请按 $2 列中的坐标合并行吗?有一系列坐标增加一个。我想输出 f.e. : 第一行合并到第四行 9079811-9079814 之后没有系列所以将它合并到另一行等。对于输入中的第三列我想计算平均值。
我写了一些脚本,但是这个脚本合并了从第一个坐标到最后一个坐标的所有行。没有条件系列。
awk -F'\t' -v OFS="\t" '{print $2,$4,$3,$1}' input | awk '!x[$2]{x[$2]=$1}y[$2]<$1{y[$2]=$1;}x[$2]>$1{x[$2]=$1} {sum+=$3} END{for(i in y)print $1,x[i],y[i],sum/NR,i}' | sort -V -k1,1 > output
输入:
chr12 9079811 29 A2M
chr12 9079812 29 A2M
chr12 9079813 29 A2M
chr12 9079814 28 A2M
chr12 9091202 5 A2M
chr12 9091203 5 A2M
chr12 9091204 5 A2M
chr12 9091390 15 A2M
chr12 9091391 15 A2M
chr12 9091392 13 A2M
输出:
chr12 9079811 9079814 28.75 A2M
chr12 9091202 9091204 5 A2M
chr12 9091390 9091392 14.3 A2M
最佳答案
Awk
解决方案:
awk 's{
if ($2-prev == 1) { sum += $3; c++; prev=$2 }
else { print $1, s, prev, sum/c, $4; s=sum=c=0 }
}
!s{ s=prev=$2; sum=$3; c++ }
END{ print $1, s, prev, sum/c, $4 }' file
s
- 指向事件处理系列的变量;包含系列的起始元素(例如 9079811
) prev
- 保存已处理系列的每个前一项 输出:
chr12 9079811 9079814 28.75 A2M
chr12 9091202 9091204 5 A2M
chr12 9091390 9091392 14.3333 A2M
关于awk - 合并系列并识别它的结尾 - AWK,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/48458008/