免费注册 查看新帖 |

Chinaunix

  平台 论坛 博客 文库
最近访问板块 发新帖
查看: 11454 | 回复: 7

[练习] 选取所需碱基 [复制链接]

论坛徽章:
7
戌狗
日期:2013-12-15 20:43:38技术图书徽章
日期:2014-03-05 01:33:12技术图书徽章
日期:2014-03-15 20:31:17未羊
日期:2014-03-25 23:48:20丑牛
日期:2014-04-07 22:37:44巳蛇
日期:2014-04-11 21:58:0915-16赛季CBA联赛之青岛
日期:2016-03-17 20:36:13
发表于 2013-11-29 05:38 |显示全部楼层
问题:

一个文件,用tab键隔开。

2 5 ATGCTGCTGTGCACATGTCAT
11 16 AGTCTCGCTCGCAACTAGCG
13 26 AGCTCAGATCGACTAGCATACGATCAGCT
7 15 GTCACATGCATCAGCATC

第一行需要从第2个到第5个碱基,共4个;
第二行需要从第11个到第16个碱基,共6个;
第三行需要从第13个到第26个碱基,共14个;
第四行需要从第7个到第15个碱基,共9个。

output:

TGCT
GCAACT
CTAGCATACGATCA
TGCATCAGC

论坛徽章:
13
双鱼座
日期:2013-10-23 09:30:05数据库技术版块每日发帖之星
日期:2016-04-20 06:20:00程序设计版块每日发帖之星
日期:2016-03-09 06:20:002015亚冠之塔什干火车头
日期:2015-11-02 10:07:452015亚冠之德黑兰石油
日期:2015-08-30 10:07:07数据库技术版块每日发帖之星
日期:2015-08-28 06:20:00数据库技术版块每日发帖之星
日期:2015-08-05 06:20:002015年迎新春徽章
日期:2015-03-04 09:57:09辰龙
日期:2014-12-03 14:45:52酉鸡
日期:2014-07-23 09:46:23亥猪
日期:2014-03-13 08:46:22金牛座
日期:2014-02-11 09:36:21
发表于 2013-11-29 09:22 |显示全部楼层
File.foreach('urfile') {|f|
  f.chomp!
  (a,b,c)=f.split
  puts c.slice(a.to_i-1,b.to_i-1)
}

论坛徽章:
5
丑牛
日期:2014-01-21 08:26:26卯兔
日期:2014-03-11 06:37:43天秤座
日期:2014-03-25 08:52:52寅虎
日期:2014-04-19 11:39:48午马
日期:2014-08-06 03:56:58
发表于 2013-11-29 13:04 |显示全部楼层
{:2_172:}
  1. DATA.each do |l|
  2.   i, j, s = l.split
  3.   i, j    = [i, j].map(&:to_i)
  4.   puts s[i-1..j-1]
  5. end

  6. __END__
  7. 2 5 ATGCTGCTGTGCACATGTCAT
  8. 11 16 AGTCTCGCTCGCAACTAGCG
  9. 13 26 AGCTCAGATCGACTAGCATACGATCAGCT
  10. 7 15 GTCACATGCATCAGCATC
复制代码

论坛徽章:
7
戌狗
日期:2013-12-15 20:43:38技术图书徽章
日期:2014-03-05 01:33:12技术图书徽章
日期:2014-03-15 20:31:17未羊
日期:2014-03-25 23:48:20丑牛
日期:2014-04-07 22:37:44巳蛇
日期:2014-04-11 21:58:0915-16赛季CBA联赛之青岛
日期:2016-03-17 20:36:13
发表于 2013-12-01 04:00 |显示全部楼层
回复 2# bikong0411

大牛,这个
  1. puts c.slice(a.to_i-1,b.to_i-1)
复制代码

论坛徽章:
7
戌狗
日期:2013-12-15 20:43:38技术图书徽章
日期:2014-03-05 01:33:12技术图书徽章
日期:2014-03-15 20:31:17未羊
日期:2014-03-25 23:48:20丑牛
日期:2014-04-07 22:37:44巳蛇
日期:2014-04-11 21:58:0915-16赛季CBA联赛之青岛
日期:2016-03-17 20:36:13
发表于 2013-12-01 04:44 |显示全部楼层
v3:{:3_188:}
  1. #!/usr/bin/ruby
  2. DATA.each do |l|
  3.   i, j, k = l.split
  4.   puts k[i.to_i - 1 .. j.to_i - 1]
  5. end

  6. __END__
  7. 2 5 ATGCTGCTGTGCACATGTCAT
  8. 11 16 AGTCTCGCTCGCAACTAGCG
  9. 13 26 AGCTCAGATCGACTAGCATACGATCAGCT
  10. 7 15 GTCACATGCATCAGCATC
复制代码

论坛徽章:
13
双鱼座
日期:2013-10-23 09:30:05数据库技术版块每日发帖之星
日期:2016-04-20 06:20:00程序设计版块每日发帖之星
日期:2016-03-09 06:20:002015亚冠之塔什干火车头
日期:2015-11-02 10:07:452015亚冠之德黑兰石油
日期:2015-08-30 10:07:07数据库技术版块每日发帖之星
日期:2015-08-28 06:20:00数据库技术版块每日发帖之星
日期:2015-08-05 06:20:002015年迎新春徽章
日期:2015-03-04 09:57:09辰龙
日期:2014-12-03 14:45:52酉鸡
日期:2014-07-23 09:46:23亥猪
日期:2014-03-13 08:46:22金牛座
日期:2014-02-11 09:36:21
发表于 2013-12-01 09:56 |显示全部楼层
回复 4# rubyish


    ruby版本不同,使用"XX"[a..b]这种得到的结果就不一样了,有的得到的是ASCII

论坛徽章:
3
寅虎
日期:2013-11-27 07:53:29申猴
日期:2014-09-12 09:24:152015年迎新春徽章
日期:2015-03-04 09:48:31
发表于 2013-12-02 09:51 |显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽

论坛徽章:
13
双鱼座
日期:2013-10-23 09:30:05数据库技术版块每日发帖之星
日期:2016-04-20 06:20:00程序设计版块每日发帖之星
日期:2016-03-09 06:20:002015亚冠之塔什干火车头
日期:2015-11-02 10:07:452015亚冠之德黑兰石油
日期:2015-08-30 10:07:07数据库技术版块每日发帖之星
日期:2015-08-28 06:20:00数据库技术版块每日发帖之星
日期:2015-08-05 06:20:002015年迎新春徽章
日期:2015-03-04 09:57:09辰龙
日期:2014-12-03 14:45:52酉鸡
日期:2014-07-23 09:46:23亥猪
日期:2014-03-13 08:46:22金牛座
日期:2014-02-11 09:36:21
发表于 2013-12-03 07:48 |显示全部楼层
回复 7# Sevk


    恩,是这样的
您需要登录后才可以回帖 登录 | 注册

本版积分规则 发表回复

  

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

清除 Cookies - ChinaUnix - Archiver - WAP - TOP