awk 处理两个文件比较复杂的比较不知道可不可以做到

Treywea

回复 30# zxy877298415


    qianwen@epilab:/data/Users/qianwen/bsmap-2.74$ awk 'FNR==NR&&/five_prime_UTR/{a[$1]=a[$1]?a[$1]","$4"#"$54"#"$5;next} FNR>1{split(a[$1],b,",";for(i=1;i<=length(b);i++){split(b,c,"#";if(($2>=c[1])&&($2<=c[2])){meth+=$7;all+=$8;CT+=$6;}}} END{print "meth="meth"\t""all="all"\t""ratio="meth/all"\t""CT="CT"\t";}' five_prime_UTR_Chr1.txt Test11.txt
meth=58 all=17172       ratio=0.00337759        CT=17167.4
额 还算了一下第8列 其实第八列是和第6列一样的 不知道为什么第六列会有小数点啊

还有就是代码中间都看不太懂 可不可以给我稍微解释解释呢 T T

zxy877298415

回复 31# Treywea
meth的结果对吗？你确定文本中的$6和$8是一样的？


   

Treywea

回复 32# zxy877298415


    嗯 结果也少了一倍数字 代码中间怎么运行的我好蒙蒙的。。。

Treywea

回复 32# zxy877298415


    不让我回复你消息。。。积分不够啊。。。