20积分悬赏，求个perl程序，新手，写了好几个都无果。

dahe_1984

rubyish 发表于 2013-03-31 07:41

你写的code怎么没有了？ 我还想学习下呢

yanglc2013

回复 51# picbhan

大哥，还有个问题请教你。
如果我想在输出的结果中的结果中，同时加上对应的ref中的数据，比如。
lnc-C1orf227-1:2        chr1        213028819        213031452        -    NM_014053        chr1        213031596        213072705        +        2210

应该怎么修改print。

现在输出的结果形式如下：
lnc-C1orf227-1:2        chr1        213028819        213031452        -


   

picbhan

回复 54# yanglc2013


    你在下面的代码中加入输出的信息（包括@rd和@id），并且把原来的三个赋值语句改为一个表示已经输出过的flag，这样在后面判断是否是non-overlap的时候就只需要看这个flag是不是真值了。

1. 
   
2.                 if ($id[2] >= $rd[2] && $id[3] <= $rd[3]) {
   
3.                     $result{ref_include_inc} = 1;
   
4.                 }
   
5.                 elsif ($rd[2] >= $id[2] && $rd[3] <= $id[3]) {
   
6.                     $result{inc_include_ref} = 1;
   
7.                 }
   
8.                 else {
   
9.                     $result{overlap} = 1;
   
10.                 }
    
11.                 if (%result) {
    
12.                 if ($result{ref_include_inc}) {
    
13.                     print {$ref_include_inc} join("\t", @id), "\n";
    
14.                 }
    
15.                 if ($result{inc_include_ref}) {
    
16.                     print {$inc_include_ref} join("\t", @id), "\n";
    
17.                 }
    
18.                 if ($result{overlap}) {
    
19.                     print {$overlap} join("\t", @id), "\n";
    
20.                 }
    
21.             }
    
22.             else { # if $result is undef, means no overlapped genes
    
23.                 print {$non_include} join("\t", @id), "\n";
    
24.             }

复制代码

yanglc2013

回复 2# picbhan

谢谢啊，你的程序非常好用。还有个难题想请教你，我也写了几个，一塌糊涂啊。

难题.rar (1.89 MB, 下载次数: 9)

2013-04-08 19:16 上传

点击文件名下载附件

   

picbhan

yanglc2013 发表于 2013-04-08 19:17
回复 2# picbhan

谢谢啊，你的程序非常好用。还有个难题想请教你，我也写了几个，一塌糊涂啊。

首先，我没有看到你说的你也写了几个，我只看到你的问题。
其次，其实你的问题很简单，如果你真认真学了Perl并且自己动手写了，肯定是可以完成的，这根本不是什么难题，仅仅只是简单的格式转换而已。
最后，没有人可以一直帮你写代码，而且这种问题都等着别人给你写代码的话只会降低你自己的效率。好好学perl基础编程吧，足以在生物里用了。

1. #!perl
   
2. use List::MoreUtils 'pairwise';
   
3. 
   
4. $fin = 'refGene.txt';
   
5. open my $in, '<', $fin or die "Can't open $fin for reading. $!\n";
   
6. $fout = 'result.txt';
   
7. open my $to, '>', $fout or die "Can't create $fout for writing. $!\n";
   
8. 
   
9. while (<$in>) {
   
10.     chomp;
    
11.     my @t = split;
    
12.     @start = split /,/, $t[1];
    
13.     @end = split /,/, $t[2];
    
14.     pairwise { print {$to} "$t[0]\t$a\t$b\n" } @start, @end;
    
15. }
    
16. close $in;
    
17. close $to;

复制代码