免费注册 查看新帖 |

Chinaunix

  平台 论坛 博客 文库
12下一页
最近访问板块 发新帖
查看: 6015 | 回复: 17
打印 上一主题 下一主题

awk下如何对行进行操作~ [复制链接]

论坛徽章:
0
跳转到指定楼层
1 [收藏(0)] [报告]
发表于 2011-10-23 11:06 |只看该作者 |倒序浏览
在选定了RS之后,如何计算每行的字符数?
$2并不能指代第二行~

论坛徽章:
2
射手座
日期:2014-10-10 15:59:4715-16赛季CBA联赛之上海
日期:2016-03-03 10:27:14
2 [报告]
发表于 2011-10-23 11:18 |只看该作者
回复 1# znbznb586


    给点数据看看

论坛徽章:
0
3 [报告]
发表于 2011-10-23 11:37 |只看该作者
ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/
从这里可以下一些数据,主要是想计算每个基因的长度~

论坛徽章:
0
4 [报告]
发表于 2011-10-23 11:38 |只看该作者
回复 2# yinyuemi

ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/
从这里可以下一些数据,主要是想计算每个基因的长度~

论坛徽章:
2
射手座
日期:2014-10-10 15:59:4715-16赛季CBA联赛之上海
日期:2016-03-03 10:27:14
5 [报告]
发表于 2011-10-23 11:53 |只看该作者
回复 4# znbznb586


    运行下,看是你要的不?
awk -vRS='>[^\n]*\n' 'NR>1{printf gensub(/>(gi\|[0-9]+\|).*/,"\\1",1,T)"\t"length(gensub(/\n/,"","g",$0))"\n"}{T=RT}' fa

论坛徽章:
0
6 [报告]
发表于 2011-10-23 11:59 |只看该作者
回复 5# yinyuemi

gensub 没有定义~

论坛徽章:
2
射手座
日期:2014-10-10 15:59:4715-16赛季CBA联赛之上海
日期:2016-03-03 10:27:14
7 [报告]
发表于 2011-10-23 12:03 |只看该作者
回复 6# znbznb586


    不知你是什么环境,用gawk,
   或者试试这个:
awk -vRS='>[^\n]*\n' 'NR>1{gsub("\n","",$0);sub(">","",T);split(T,a,"[|]");printf a[1]"|"a[2]"|\t"length"\n"}{T=RT}' file

论坛徽章:
0
8 [报告]
发表于 2011-10-23 12:12 |只看该作者
出来结果了,太感谢了~
{T=RT}什么意思?去掉后好像没区别

论坛徽章:
0
9 [报告]
发表于 2011-10-23 12:14 |只看该作者
回复 7# yinyuemi


    出来结果了,太感谢了~
{T=RT}什么意思?去掉后好像没区别

论坛徽章:
2
射手座
日期:2014-10-10 15:59:4715-16赛季CBA联赛之上海
日期:2016-03-03 10:27:14
10 [报告]
发表于 2011-10-23 12:15 |只看该作者
回复 8# znbznb586


    是可以出来长度,但是没有gi number了
   

  RT 是RS匹配到的部分
man awk
您需要登录后才可以回帖 登录 | 注册

本版积分规则 发表回复

  

北京盛拓优讯信息技术有限公司. 版权所有 京ICP备16024965号-6 北京市公安局海淀分局网监中心备案编号:11010802020122 niuxiaotong@pcpop.com 17352615567
未成年举报专区
中国互联网协会会员  联系我们:huangweiwei@itpub.net
感谢所有关心和支持过ChinaUnix的朋友们 转载本站内容请注明原作者名及出处

清除 Cookies - ChinaUnix - Archiver - WAP - TOP