文件中对应列的处理

gooderpan

有一个文件file：（tab分隔）
Chr    Pos    Ref     01N    01T    02N   02T    03N   03T
1       123    A         a       c         b       c        a      d
1       223    G         c       c         a       c        d      c
1       3234   T         b       d        c       d        a      c
2       123    A         a       d         a       d        b      c
3       123    A         b       c         b       d        c      c

我们把01N和01T，02N和02T，03N和03T看作是一对的对照。每一行读取时如果同一对中N出现了“c”或“d”，那么将这一对N和T的值都替换成“e”
Chr    Pos    Ref     01N    01T    02N   02T    03N   03T
1       123    A         a       c         b       c        a      d
1       223    G         e       e         a       c        e      e
1       3234   T         b       d        e       e        a      c
2       123    A         a       d         a       d        b      c
3       123    A         b       c         b       d        e      e

请问大神们，这个该怎么办呀？？？

jason680

回复 1# gooderpan

http://bbs.chinaunix.net/thread-4241164-1-1.html
   

gooderpan

大神，这个例子中首先需要对于每一对中N那一列取值的判断，该怎么加入呀？？回复 2# jason680


   

zxy877298415

1. awk 'NR>1{for(i=4;i<=NF;i++) {a[i]=$i} for(i=4;i<NF;i+=2)
   
2. {if((a[i]=="c"&&a[i+1]=="d")||(a[i]=="c"&&a[i+1]=="c")||(a[i]=="d"&&a[i+1]=="c")||(a[i]=="d"&&a[i+1]=="d")) {sub(a[i],"e",$i);sub(a[i+1],"e",$(i+1))}}}1' file |column -t

复制代码

回复 1# gooderpan


   

elu_ligao

1. awk 'NR==1{print; next}{for(i=4;i<NF;i+=2){if($i~/[cd]/)$i=$(i+1)="e"}}1' file | column -t

复制代码

jcdiy0601

awk '{if($4=="c"||$4=="d")$4=$5="e";if($6=="c"||$6=="d")$6=$7="e";if($8=="c"||$8=="d")$8=$9="e";print}' file |column -t
Chr  Pos   Ref  01N  01T  02N  02T  03N  03T
1    123   A    a    c    b    c    a    d
1    223   G    e    e    a    c    e    e
1    3234  T    b    d    e    e    a    c
2    123   A    a    d    a    d    b    c
3    123   A    b    c    b    d    e    e

我是一隻羊

1. awk   -vNum='4|6|8' 'BEGIN{split(Num,array,"|")}{for(i in array){$array[i]~/c|d/?$array[i]=$(array[i]+1)="e":1}}1'|column -t

复制代码

moperyblue

1. sed -r ':a;s/^((\S+\s+){2}\S+((\s+\S+){2})*)(\s+)[cd](\s+)\S+/\1\5e\6e/;ta'

复制代码