如何将文件的行排序为矩阵?

如何将文件的行排序为矩阵?

我想知道如何按以下方式对文件的行进行排序:

我的文件是 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

相关内容