免费注册 查看新帖 |

Chinaunix

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

请各位大侠帮忙处理用perl脚本处理文件 [复制链接]

论坛徽章:
0
跳转到指定楼层
1 [收藏(0)] [报告]
发表于 2014-10-23 17:49 |只看该作者 |倒序浏览
附件的文件的第一例为样本的ID数,最后一列是分类,中间的的数字是相对丰度,现在我想将中间的数字都小于0.0001的那一行过滤掉,请各位帮忙,谢谢

QQ图片20141023174817.jpg (97.7 KB, 下载次数: 37)

QQ图片20141023174817.jpg

论坛徽章:
8
技术图书徽章
日期:2013-08-22 11:21:28未羊
日期:2015-01-19 22:22:25巳蛇
日期:2014-08-11 16:53:08子鼠
日期:2014-05-29 09:04:44摩羯座
日期:2014-04-11 14:15:07丑牛
日期:2014-01-24 12:41:28金牛座
日期:2013-11-21 17:38:28射手座
日期:2015-01-21 08:50:32
2 [报告]
发表于 2014-10-23 20:33 |只看该作者
print if /\t([0-9.]+)\t/ and $1 >0.0001;?

论坛徽章:
0
3 [报告]
发表于 2014-10-23 22:28 |只看该作者
您好,请问9后面为什们加.呢?谢谢回复 2# huang6894


   

论坛徽章:
8
技术图书徽章
日期:2013-08-22 11:21:28未羊
日期:2015-01-19 22:22:25巳蛇
日期:2014-08-11 16:53:08子鼠
日期:2014-05-29 09:04:44摩羯座
日期:2014-04-11 14:15:07丑牛
日期:2014-01-24 12:41:28金牛座
日期:2013-11-21 17:38:28射手座
日期:2015-01-21 08:50:32
4 [报告]
发表于 2014-10-24 08:59 |只看该作者
回复 3# huangluhao


    加.表示小数点也考虑在内,我认真看了一下你的需求,发现我的答案是错的

论坛徽章:
32
处女座
日期:2013-11-20 23:41:20双子座
日期:2014-06-11 17:20:43戌狗
日期:2014-06-16 11:05:00处女座
日期:2014-07-22 17:30:47狮子座
日期:2014-07-28 15:38:17金牛座
日期:2014-08-05 16:34:01亥猪
日期:2014-08-18 13:34:25白羊座
日期:2014-09-02 15:03:55金牛座
日期:2014-11-10 10:23:58处女座
日期:2014-12-02 09:17:52程序设计版块每日发帖之星
日期:2015-06-16 22:20:002015亚冠之塔什干火车头
日期:2015-06-20 23:28:22
5 [报告]
发表于 2014-10-24 10:56 |只看该作者
现在我想将中间的数字都小于0.0001的那一行过滤掉

你确定?{:3_183:}

论坛徽章:
0
6 [报告]
发表于 2014-10-24 12:51 |只看该作者
确定的,就是将每行的数字都小于0.0001的滤掉回复 5# yestreenstars


   

论坛徽章:
8
技术图书徽章
日期:2013-08-22 11:21:28未羊
日期:2015-01-19 22:22:25巳蛇
日期:2014-08-11 16:53:08子鼠
日期:2014-05-29 09:04:44摩羯座
日期:2014-04-11 14:15:07丑牛
日期:2014-01-24 12:41:28金牛座
日期:2013-11-21 17:38:28射手座
日期:2015-01-21 08:50:32
7 [报告]
发表于 2014-10-24 13:10 |只看该作者
回复 6# huangluhao


    额。。。还没解决啊?
好吧,菜鸟给你一个先用着吧:
  1. #! /usr/bin/perl -w
  2. while(<DATA>){
  3.         my @all = /\t([0-9.]+)/g;
  4.         pop @all;
  5.         print $_ if grep{$_ >0.001}@all;
  6. }
  7. __DATA__
  8. 1       0       0       0       2       0       1
  9. 2       0.001   0       0       0       0       2
  10. 3       0       0       0.1     0       0       1
  11. 4       0       0.2     0.1     0       0       3
复制代码

论坛徽章:
0
8 [报告]
发表于 2014-10-24 14:53 |只看该作者
How about this one
  1. while(<DATA>)
  2. {
  3.         chomp;
  4.         print $_,"\n" if(eval(join("+",map{$_>=0.01 ? 1 : 0}(split(/\t/))[1..scalar((split(/\t/)))-2])));
  5. }

  6. __DATA__
  7. 1        0        0        0        2        0        1
  8. 2        0.001        0        0        0        0        2
  9. 3        0        0        0.1        0        0        1
  10. 4        0        0.2        0.1        0        0        3
复制代码

论坛徽章:
0
9 [报告]
发表于 2014-10-24 14:53 |只看该作者
How about this one
  1. while(<DATA>)
  2. {
  3.         chomp;
  4.         print $_,"\n" if(eval(join("+",map{$_>=0.01 ? 1 : 0}(split(/\t/))[1..scalar((split(/\t/)))-2])));
  5. }

  6. __DATA__
  7. 1        0        0        0        2        0        1
  8. 2        0.001        0        0        0        0        2
  9. 3        0        0        0.1        0        0        1
  10. 4        0        0.2        0.1        0        0        3
复制代码

论坛徽章:
0
10 [报告]
发表于 2014-10-24 16:28 |只看该作者
谢谢回复 7# huang6894


   
您需要登录后才可以回帖 登录 | 注册

本版积分规则 发表回复

  

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

清除 Cookies - ChinaUnix - Archiver - WAP - TOP