輸出匹配列

gooderpan

有一個文件file：（tab分隔）
Chrom     pos     ref      alt     T01     N01      T02    N02
1            123      A       B       0/1      ./.        1/1     0/0
2            1344    A       G       1/1      0/0      0/1      ./.
3            14553  T       C        1/1      ./.       1/1      0/0

在第一行中搜索并定位T02屬于哪一列，輸出：（事先不知道T02是哪一列）
Chrom     pos     ref      alt      T02
1            123      A       B        1/1
2            1344    A       G        0/1
3            14553  T       C         1/1

請問大神該怎么辦？

haooooaaa

1. awk 'NR==1{for(i=5;i<=NF;i++)if($i=="T02")a=i}{$5=$a;NF=5}1' file | column -t

復(fù)制代碼

gooderpan

謝謝大神��！回復(fù) 2# haooooaaa


   

gooderpan

請問，如果最后結(jié)果只要求輸出第2，4，5列，那么原命令怎么改呢？？回復(fù) 2# haooooaaa


   

jason680

回復(fù) 4# gooderpan

第2，4，5列  ==> $2,$4,$5
   

gooderpan

awk 'NR==1{for(i=5;i<=NF;i++)if($i=="T02")a=i}{$5=$a;NF=$2,$4,$5}1' file ??? 貌似不行噢回復(fù) 5# jason680


   

elu_ligao

1. {$5=$a;NF=$2,$4,$5}1
   
2. ->
   
3. {$5=$a;print $2,$4,$5}

復(fù)制代碼

gooderpan

請問大神，“{$5=$a;NF=$2,$4,$5}1” 中的“1”代表是啥意思呀？？回復(fù) 7# elu_ligao


   

elu_ligao

回復(fù) 8# gooderpan


    echo aaa | awk '{print}'
    echo aaa | awk '1'

    都是打印全行內(nèi)容

hz_oracle

大神好膩害