我想知道如何按以下方式对文件的行进行排序:
我的文件是 file.txt (制表符分隔):
g1 00A98_01563 00554_01552 CCUG38_01373
g2 00554_01444
g3 00A98_04566 CCUG38_05322
我想要得到这个(制表符分隔):
g 00A98 00554 CCUG38
g1 1 1 1
g2 0 1 0
g3 1 0 1
或采用以下格式(制表符分隔):
g 00A98 00554 CCUG38
g1 00A98_01563 00554_01552 CCUG38_01373
g2 00554_01444
g3 00A98_04566 CCUG38_05322
我如何使用带有 sort、awk、grep 或其他命令行的命令行来执行此操作?
祝一切顺利,问候
答案1
使用磨坊主,将输入视为以 TAB 作为输入字段分隔符、以下划线作为输入对分隔符的分隔键值对,并将输出设置为 TSV。然后您可以不稀疏化您的数据
$ mlr --idkvp --ifs tab --ips '_' --otsv unsparsify file.txt
1 00A98 00554 CCUG38
g1 01563 01552 01373
g2 01444
g3 04566 05322
然后您可以添加各种转换,例如
$ mlr --idkvp --ifs tab --ips '_' --otsv unsparsify --fill-with 0 then put '
for(k,v in mapexcept($*,"1")){if(v != 0){$[k] = 1}}
' then rename "1","g" file.txt
g 00A98 00554 CCUG38
g1 1 1 1
g2 0 1 0
g3 1 0 1
或者
$ mlr --idkvp --ifs tab --ips '_' --otsv unsparsify then put -S '
for(k,v in mapexcept($*,"1")){if(v != ""){$[k] = k ."_". v}}
' then rename "1","g" file.txt
g 00A98 00554 CCUG38
g1 00A98_01563 00554_01552 CCUG38_01373
g2 00554_01444
g3 00A98_04566 CCUG38_05322
在最后一种情况下,对齐看起来“不正确”,但用--ocsv
代替输出--tsv
应该可以确认它是正确的。
答案2
这是toMatrix.awk
#!/usr/bin/env gawk -f
BEGIN { FS = OFS = "\t" }
{
for (i=2; i<=NF; i++) {
x=$i
sub(/_.*/, "", x)
if (!(x in values)) {
values[x] = 1
ordered[++value] = x
}
g[NR] = $1
data[NR][x]=1
}
}
END {
printf "%s", "g"
for (i = 1; i <= value; i++)
printf "%s%s", OFS, ordered[i]
print ""
for (nr = 1; nr <= NR; nr++) {
printf "%s", g[nr]
for (i = 1; i <= value; i++)
printf "%s%s", OFS, 0 + data[nr][ordered[i]]
print ""
}
}
$ gawk -f toMatrix.awk file.txt
g 00A98 00554 CCUG38
g1 1 1 1
g2 0 1 0
g3 1 0 1