最新方法請(qǐng)看：http://www.itdecent.cn/p/cc31cb7c614d

此方法能達(dá)到目的但不是最優(yōu)方法

說(shuō)明

需要兩個(gè)文件，f1為需要提取的列名，自上而下為提取順序；f2為待提取的矩陣，第一列默認(rèn)輸出。

f1示例

image.png

f2示例

image.png

腳本

awk 'BEGIN{FS=OFS="\t"}NR==FNR{a[NR]=$1}NR>FNR&&FNR==1{for(i=1;i<length(a)+1;i++){for(j=1;j<NF+1;j++){if(a[i]==$j){b[j]=$j}}}}NR>FNR{printf $1;for(i=1;i<length(a)+1;i++){for(j=1;j<NF+1;j++){if(a[i]==b[j]){printf "\t"$j}}}print ""}'  pheno.file  matrix.02

即

awk '
BEGIN{
  FS=OFS="\t"
}NR==FNR{
  a[NR]=$1
}NR>FNR&&FNR==1{
  for(i=1;i<length(a)+1;i++){
    for(j=1;j<NF+1;j++){
      if(a[i]==$j){
        b[j]=$j
       }
    }
  }
}NR>FNR{
  printf $1
  for(i=1;i<length(a)+1;i++){
    for(j=1;j<NF+1;j++){
      if(a[i]==b[j]){
        printf "\t"$j
      }
    }
  }
  print ""
}' pheno.file  matrix.02

結(jié)果

image.png