2

我试图让它以逗号分隔的形式输出。当前版本根本不起作用(我得到一个空白文件作为输出),而以前的版本(我保留 awk BEGIN 语句但没有排序分隔符)只会以制表符分隔而不是逗号分隔输出. 在以前的版本中,在不尝试获取逗号分隔符的情况下,我确实得到了预期的答案(使用复杂的过滤器等),所以我不寻求这部分的帮助。我意识到这是一种非常丑陋的过滤方式,而且数字也很丑陋/非常大。

问题背景:在文件 lamina.bed 中找到与区域 chr12:5000000-6000000 重叠的区域,并按第 4 列降序排序,输出为逗号分隔。染色体是第一列,区域的起始位置是第 2 列,结束位置是第 3 列,值是第 4 列。我们应该使用 awk(在 Unix bash shell 中)。预先感谢您的帮助!

awk 'BEGIN{FS="\t"; OFS=","} ($2 <= 5000000 && $3 >= 5000000) || ($2 >= 5000000 && $3 <= 6000000) || ($2 <= 6000000 && $3 >= 6000000) || ($2 <= 5000000 && $3 >= 6000000)' /vol1/opt/data/lamina.bed | awk 'BEGIN{FS=","; OFS=","} ($1 == "chr12") ' | sort -t$"," -k4rn > ~/MOLB7621/PS_2/results/2015_02_05/PS2_p3_n1.csv
cat ~/MOLB7621/PS_2/results/2015_02_05/PS2_p3_n1.csv

输入示例行(制表符分隔,包括 chr12 上应该工作的行):

#chrom  start   end value
chr1    11323785    11617177    0.86217008797654
chr1    12645605    13926923    0.934891485809683
chr1    14750216    15119039    0.945945945945946
chr12   3306736 5048326 0.913561847988077
chr12   5294045 5393088 0.923076923076923
chr12   5505370 6006665 0.791318864774624
chr12   7214638 7827375 0.8562874251497
chr12   8139885 10173149    0.884353741496599
4

1 回答 1

2

要获得逗号分隔的输出,请使用以下命令:

$ awk 'BEGIN{FS="\t"; OFS=","} ($2 <= 5000000 && $3 >= 5000000) || ($2 >= 5000000 && $3 <= 6000000) || ($2 <= 6000000 && $3 >= 6000000) || ($2 <= 5000000 && $3 >= 6000000) {$1=$1;print}' file | awk 'BEGIN{FS=","; OFS=","} ($1 == "chr12") ' | sort -t$"," -k4rn 
chr12,5294045,5393088,0.923076923076923
chr12,3306736,5048326,0.913561847988077
chr12,5505370,6006665,0.791318864774624

上面唯一的变化是对动作的添加:

{$1=$1;print}

awk如果行上的一个或多个字段已以某种方式更改,则只会使用新的字段分隔符重新格式化该行。 $1=$1足以表明字段 1 已更改。因此,插入了新的字段分隔符。

此外,这两个调用awk可以合并为一个调用:

awk 'BEGIN{FS="\t"; OFS=","} ($2 <= 5000000 && $3 >= 5000000) || ($2 >= 5000000 && $3 <= 6000000) || ($2 <= 6000000 && $3 >= 6000000) || ($2 <= 5000000 && $3 >= 6000000) {$1=$1; if($1 == "chr12") print}' file | sort -t$"," -k4rn

更简单的例子

在下文中,输入以制表符分隔,输出字段分隔符 ,OFS设置为逗号。在第一个示例中,使用了以下awk命令print

$ echo $'a\tb\tc' | awk -v OFS=, '{print}'
a       b       c

尽管如此OFS=,,输出仍保留制表符分隔符。

现在,我们添加简单的语句$1=$1并观察输出:

$ echo $'a\tb\tc' | awk -v OFS=, '{$1=$1;print}'
a,b,c

输出现在以逗号分隔。同样,这是因为只有在它认为行上的某个字段已以某种方式更改awk时才使用新行重新格式化。OFSto 本身的分配$1足以触发重新格式化。

请注意,仅进行影响整个生产线的更改是不够的。例如,以下内容不会触发重新格式化:

$ echo $'a\tb\tc' | awk -v OFS=, '{$0=$0;print}'
a       b       c

需要单独更改该行的一个或多个字段。在下面,作为一个整体sub运行$0,因此不会触发重新格式化:

$ echo $'a\tb\tc' | awk -v OFS=, '{sub($1,"NEW");print}'
NEW     b       c

但是,在下面的示例中,sub专门针对字段进行操作$1,因此会触发重新格式化:

$ echo $'a\tb\tc' | awk -v OFS=, '{sub($1,"NEW", $1);print}'
NEW,b,c
于 2015-02-08T23:42:38.010 回答