多谢大家了,另外一个问题awk 提取匹配

liaoboda

比如边边内容(标红),我要提取第一行和后边匹配KXXX
Ame_R016766        "PPP2R3A; protein phosphatase 2, regulatory subunit B'', alpha; K11583 protein phosphatase 2 (formerly 2A), regulatory subunit B''"
Ame_R013656        "CDH13; cadherin 13, H-cadherin (heart); K06808 cadherin 13, H-cadherin"
Ame_R018461        NOTCH4; notch 4; K02599 Notch
Ame_R018610        rho guanine nucleotide exchange factor 11-like; K12331 Rho guanine nucleotide exchange factor 11
Ame_R015728        GL19043 gene product from transcript GL19043-RA; K07295 nuclear receptor subfamily 2 group A

liaoboda

有可能有些Ame 后边没有这个K值

liaoboda

Ame_R019376     glutathione K00799
Ame_R019673     similar K04257
Ame_R020529     ADAM21 K08610
Ame_R020683     "CDH23, K06813
Ame_R018879     "SCN10A K04842
Ame_R019377     similar K00799
Ame_R020632     "IGLL1, K06554
Ame_R020201     threonine K01733
Ame_R018946     NRXN3 K07377
Ame_R020727     similar K06553
Ame_R018801     "ABCA1 K05641
Ame_R017225     "similar K04700
Ame_R020267     forkhead K09401
这个前边的去不掉,我尝试了K* K0效果都不是太好,只要K06553这样的

liaoboda

我把" ; ,这些符号都去掉试试

liaoboda

awk  '{for(i=1;i<=NF;i++){if($i~/\<K0/)print $1,$i}}' Negative.above-Fst-result.txt
这个是可以了,但是有一个问题,会有空值存在,输入的有1051行,这个的结果只有715行.所以空的可能要输出$1 "none"

liaoboda

awk  '{for(i=1;i<=NF;i++){if($i~/\<K/){print $1,$i} else {print $1,"none"}}}' Negative.above-Fst-result.txt|wc
这个也会有问题,输出太多

yuloveban

回复 1# liaoboda


   没有看明白问题所在，但是地下的回复写得很精彩，学习了！

yinyuemi

回复 4# liaoboda

1. $ echo 'Ame_R016766        "PPP2R3A; protein phosphatase 2, regulatory subunitB'', alpha; K11583 protein phosphatase 2 (formerly 2A), regulatory subunit B''"
   
2. Ame_R013656        "CDH13; cadherin 13, H-cadherin (heart); K06808 cadherin 13, H-cadherin"
   
3. Ame_R018461        NOTCH4; notch 4; K02599 Notch
   
4. Ame_R018610        rho guanine nucleotide exchange factor 11-like; K12331 Rho guanine nucleotide exchange factor 11
   
5. Ame_R015728        GL19043 gene product from transcript GL19043-RA;nuclear receptor subfamily 2 group A' |\
   
6. awk '{printf $1 FS}{print /; K[0-9]+ /?gensub(/.*; (K[0-9]+).*/,"\\1","g",$0):"none"}'
   
7. Ame_R016766 K11583
   
8. Ame_R013656 K06808
   
9. Ame_R018610 K12331
   
10. Ame_R015728 none

复制代码