| 查看: 1347 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
perl 大神求助 已有1人参与
|
||
|
亲们，我的两个文件是这样的 1文件chr1 - 9901.5 chr1 - 144612.5 chr1 - 227317.5 chr1 - 227321.5 chr1 + 267647.5 chr1 - 267648.5 chr1 - 267649.5 chr2 + 267651.5 chr2 - 267652.5 chr3 - 521268.5 chr3 - 569831.5 chr1 + 713294.5 chr1 - 713307.5 chr1 - 753396.5 chr1 + 762732.5 chr3 - 801094.5 chr1 + 808622.5 chr2 - 808633.5 chr1 - 819067.5 3文件 chr1 - 9802 144712 chr1 - 227222 227421 chr1 - 267552 267751 chr2 + 267552 267753 处理上边两个文件的预期目标是 chr1 - 9802 144712 3 chr1 - 227222 227421 2 chr1 - 267552 267751 2 chr2 + 267552 267753 0 我想在两个文件第一列和第二列相同的情况下，如果1文件的第三列有一个处于3文件的第三列和第四列之间的情况下，$n就加1。有两个就加2，以此类推，$n初始值为0.其中的3,2，2,0是$n. my %hash; my $n=0; open OUTFILE,'>','result1'; open FILE,'<','1'; while (<FILE>  {chomp; ($a,$b,$c)=(split)[0,1,2]; $hash{$a}=1; $hash{$b}=2; open FILE,'<','3'; while (<FILE> {chomp; ($d,$e,$f,$g)=(split)[0,1,2,3]; if(exists $hash{$d} && exists $hash{$e}){ @array=sort{$a <=> $b}($f,$g); next if $c<$array[0] ; next if $c>$array[1]; $n++; }print OUTFILE"$d $e $f $g $n\n"; } } 这是我编写的程序，大家能不能帮我看看是哪里出了问题啊或者帮我重新编写一个，甚是感谢各位大牛了。 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 络合铁相关问题，求大神赐教！
已经有3人回复
- 聚酰亚胺薄膜检测项目中的耐热性标准问题
已经有18人回复
- 想询问一下“某基因”和“某基因-like”的差别? 跪求大神详细解答！
已经有4人回复
- 各位大神求赐教
已经有4人回复
- 编程菜鸟求大神引导入门。
已经有14人回复
- 【求助】请大神帮我翻译一句话
已经有3人回复
- 生物信息学问题，染色体畸变查找
已经有9人回复
- 修改图标问题，请各位大神帮忙！
已经有5人回复
- 编程有点问题，求助
已经有5人回复
- 求大神们雪中送炭
已经有12人回复
- 求助大神帮我看看这段程序是干什么的？（就是定义这个函数算什么）
已经有6人回复
- 求助 一个perl代码修改~~看看吧...
已经有6人回复
- 求大神做一个汇编语言程序
已经有7人回复
- perl程序求助，请大神们赐教~~~
已经有9人回复
- Rg(x) vs time ,Rg(y) vs time,Rg(z) vs time
已经有6人回复
- 请帮忙编写一个perl程序处理文本
已经有26人回复
- 求C语言大神帮忙
已经有6人回复
- 请教perl脚本的一个if语句解释，谢谢
已经有8人回复
- Perl 编程 请大神们帮忙！！！生物信息的牛人请进
已经有11人回复
- perl脚本求助，我想将一万多个文件从fasta格序，改为phylip4格式
已经有8人回复
9楼2014-11-25 09:44:01
wizardfan
至尊木虫 (著名写手)
- BioEPI: 18
- 应助: 599 (博士)
- 贵宾: 1.818
- 金币: 24648.2
- 散金: 197
- 红花: 48
- 沙发: 2
- 帖子: 2254
- 在线: 400.7小时
- 虫号: 1879241
- 注册: 2012-07-05
- 性别: GG
- 专业: 生物信息学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
李靖: 金币+15, ★有帮助 2014-11-25 09:42:42
感谢参与，应助指数 +1
李靖: 金币+15, ★有帮助 2014-11-25 09:42:42
|
首先先让我理解你的算法： 根据你说的 chr1 - 9802 144712 第3，4列是9802和144712，那区间就是9802-144712 在1文件中，只有第二行chr1 - 144612.5 在这个区间，那你的结果为什么是3呢？ 其次，请确定1文件和3文件的分割符是空格，还是tab 接着 看上去第二列的正反链信息好像没有用的 最后 请在你的程序上加上注释，我可以帮你看看，你到底错在哪里： 是算法错误，还是perl实现错误 |
2楼2014-11-10 17:36:47
|
不好意思，昨天把1文件上传错了，如果按一文件的话，预期的目标是 chr1 - 9802 144712 2 chr1 - 227222 227421 2 chr1 - 267552 267751 2 chr2 + 267552 267753 1 我解释一下，我想干什么(上边文件1的最开头还有一行chr1 - 9901.5 )，我另$n的初始值为0，如果文件1的第三列有落在文件三的第三列和第四列之间，那么$n就加1，以此类推， chr1 - 9901.5 chr1 - 144612.5 如上边两行的第三列都在9802 144712 之间，所以$n就为2.同样下边也是这样的也就是生物中数read个数。可以帮我编写一个也是可以的，因为帮别人改程序确实很麻烦的，非常感谢 |
3楼2014-11-11 08:34:54
wizardfan
至尊木虫 (著名写手)
- BioEPI: 18
- 应助: 599 (博士)
- 贵宾: 1.818
- 金币: 24648.2
- 散金: 197
- 红花: 48
- 沙发: 2
- 帖子: 2254
- 在线: 400.7小时
- 虫号: 1879241
- 注册: 2012-07-05
- 性别: GG
- 专业: 生物信息学
4楼2014-11-11 08:39:55
30