统计频率

catloverat99

请教各位大神，我有一个文件如下格式，我想统计每个ID下，各个因子的频率。先谢谢了
>g1
V$P53_04 287 20 306 g1_ENSRNOG00000013436_451a Yes
V$P53_04 287 20 306 g1_ENSRNOG00000013436_451a Yes
V$ZFP161_04 288 14 301 g1_ENSRNOG00000013436_451a Yes
V$ZFP161_04 290 14 303 g1_ENSRNOG00000013436_451a Yes
V$CHCH_01 292 6 297 g1_ENSRNOG00000013436_451a Yes
V$ZFP161_04 294 14 307 g1_ENSRNOG00000013436_451a Yes
V$AP2ALPHA_03 304 15 318 g1_ENSRNOG00000013436_451a Yes
V$AP2ALPHA_03 304 15 318 g1_ENSRNOG00000013436_451a Yes
V$FREAC3_01 317 16 332 g1_ENSRNOG00000013436_451a Yes

>g2
V$CREBP1_01 124 8 131 g1_ENSRNOG00000016371_102a Yes
V$CREBP1_01 124 8 131 g1_ENSRNOG00000016371_102a Yes
V$SOX2_Q3_01 144 15 158 g1_ENSRNOG00000016371_102a Yes
V$SOX10_Q3 145 8 152 g1_ENSRNOG00000016371_102a Yes
V$CDX2_01 160 16 175 g1_ENSRNOG00000016371_102a Yes

结果统计每个因子发生的频率，结果应如下

>g1
V$P53_04 2
V$ZFP161_04 2
...
>g2
V$CREBP1_01 2
V$SOX2_Q3_01 1
..

catloverat99

在线等

reyleon

gawk 4.0+

1. awk '/^>/{k=$0;next}NF{a[k][$1]++}END{for(i in a){print i;for(j in a[i])print j,a[i][j]}}' file

复制代码

catloverat99

回复 3# reyleon


运行后，报错，不知道是哪里出了问题？

awk '/^>/{k=$0;next}NF{a[k][$1]++}END{for(i in a){print i;for(j in a)print j,a[j]}}' g2a_tf_yes
awk: syntax error at source line 1
context is
         >>> /^>/{k=$0;next}NF{a[k][ <<<
awk: illegal statement at source line 1
awk: illegal statement at source line 1



   

catloverat99

回复 3# reyleon


gwak要自己安装吗

   

zxy877298415

1. awk -v RS=">" 'NR>1{print ">"$1;for(i=2;i<=NF;i+=6) {a[$i]++} for(i in a) {print i,a[i]} delete a;}' file            
   
2. >g1
   
3. V$CHCH_01 1
   
4. V$ZFP161_04 3
   
5. V$AP2ALPHA_03 2
   
6. V$P53_04 2
   
7. V$FREAC3_01 1
   
8. >g2
   
9. V$CDX2_01 1
   
10. V$SOX2_Q3_01 1
    
11. V$SOX10_Q3 1
    
12. V$CREBP1_01 2

复制代码

回复 1# catloverat99


   

catloverat99

回复 6# zxy877298415


    请问i+＝6是什么意思？我跑了我的例子没问题，可是换了个大点的文件结果就不是我想要的了，请帮忙看看哪里出了问题？ 哎，照猫画虎都不成。。。


文件的样子
>g1_ENSRNOG00000013436_451a  0 -1 g1_ENSRNOG00000013436_451a Yes
V$GFI1_Q6_01 12 10 21 g1_ENSRNOG00000013436_451a Yes
V$SMAD4_Q6_01 29 7 35 g1_ENSRNOG00000013436_451a Yes
V$GKLF_Q4 34 7 40 g1_ENSRNOG00000013436_451a Yes
V$YY1_Q6_03 37 7 43 g1_ENSRNOG00000013436_451a Yes
V$NANOG_01 48 12 59 g1_ENSRNOG00000013436_451a Yes
V$CDX2_01 63 16 78 g1_ENSRNOG00000013436_451a Yes
V$BBX_04 64 17 80 g1_ENSRNOG00000013436_451a Yes
V$BBX_04 65 17 81 g1_ENSRNOG00000013436_451a Yes
V$FREAC3_01 65 16 80 g1_ENSRNOG00000013436_451a Yes
V$AIRE_01 71 26 96 g1_ENSRNOG00000013436_451a Yes
V$TEF1_Q6_04 75 11 85 g1_ENSRNOG00000013436_451a Yes
V$IK_Q5_01 78 7 84 g1_ENSRNOG00000013436_451a Yes
V$YY1_Q6_03 79 7 85 g1_ENSRNOG00000013436_451a Yes
V$IK_Q5_01 83 7 89 g1_ENSRNOG00000013436_451a Yes
V$RFX1_01 107 17 123 g1_ENSRNOG00000013436_451a Yes
V$CHCH_01 111 6 116 g1_ENSRNOG00000013436_451a Yes
V$MYB_05 114 17 130 g1_ENSRNOG00000013436_451a Yes
V$RFX_Q6 115 9 123 g1_ENSRNOG00000013436_451a Yes
V$NANOG_01 116 12 127 g1_ENSRNOG00000013436_451a Yes
V$MYB_Q4 116 12 127 g1_ENSRNOG00000013436_451a Yes
V$CPBP_Q6 120 7 126 g1_ENSRNOG00000013436_451a Yes
V$TBX5_01 120 12 131 g1_ENSRNOG00000013436_451a Yes
V$CHCH_01 121 6 126 g1_ENSRNOG00000013436_451a Yes
V$CDX2_01 122 16 137 g1_ENSRNOG00000013436_451a Yes
V$TATA_01 127 15 141 g1_ENSRNOG00000013436_451a Yes
V$CDX2_Q5_02 128 6 133 g1_ENSRNOG00000013436_451a Yes
V$TBX5_01 129 12 140 g1_ENSRNOG00000013436_451a Yes

跑出来的结果
awk  -v RS=">" 'NR>1{print ">"$1;for(i=2;i<=NF;i+=6) {a[$i]++} for(i in a) {print i,a} delete a;}' g1a_tf_yes
>g1_ENSRNOG00000013436_451a
6 5
26 1
7 15
8 3
9 2
10 11
11 4
12 9
13 1
14 6
15 7
16 4
17 5
19 1
20 4
0 1
>g1_ENSRNOG00000016371_102a
5 2
6 19
7 50
8 47

zxy877298415

回复 7# catloverat99
你的文件的格式都变了啊，这个文本哪个是ID？


   

catloverat99

回复 8# zxy877298415

>g1_ENSRNOG00000013436_451a  0 -1 g1_ENSRNOG00000013436_451a Yes 变成>g1


   

catloverat99

回复 8# zxy877298415


    因为原来的文件很长，我挑了几行，把>id换了个更短的。