这样的程序perl能不能做？

bioinfor

珞水的大叔的程序很快很强大，我顶死你。

liyangole

回复 16# 珞水的大叔


    ok,了解了。多谢。我也对正则感兴趣，到时候跟你交流哈！

liyangole

回复 15# 珞水的大叔


    # for(@A_choose,@B){         
#         $m{$_}++ && $n{$_}++;
#     }      
#     say scalar keys %n;
# }

大叔的这段代码。我在最近帖子中看到有人用过。真是不错，又是并集又是交集，太神奇了。这个用的好，赞！

liyangole

my @A_choose = map{ $A[$_] } keys %rand;

大叔是不是这个就不用调用 rand随机函数了 ？

珞水的大叔

回复 24# liyangole


    %rand里面存放的是随机后的不重复的结果，随机的是这段：

1. while(scalar keys %rand < $yangbengeshu){
   
2.         $rand{int(rand scalar @A)} = 0;
   
3.     }
   

复制代码

liyangole

回复 25# 珞水的大叔


    多谢 ，看到了。

bioinfor

大叔的这个程序在小样本的时候没看出问题，但到了大样本的时候似乎有些不对，如附件给的数据，当 $yangbengeshu = 3时，结果居然都是5或6，理论上应该是＜＝3的数值，不知问题出在哪里呢？

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. use 5.010;
   
5. 
   
6. 
   
7. my $yangbengeshu = 3;   #抽取样本个数
   
8. 
   
9. my @A;
   
10. open A,"<human_gene.txt";
    
11. while(<A>){
    
12.     chomp;
    
13.     s/^\s+$//;
    
14.     s/^\s+//;
    
15.     s/\s+$//;
    
16.     push @A,$_ if $_ ne '';
    
17. }
    
18. close A;
    
19. 
    
20. if($yangbengeshu > scalar @A){
    
21.     die "ERROR: 抽取样本个数大于总个数\n";
    
22. }
    
23. 
    
24. my @B;
    
25. open B,"<targetgene.txt";
    
26. while(<B>){
    
27.     chomp;
    
28.     s/^\s+$//;
    
29.     s/^\s+//;
    
30.     s/\s+$//;
    
31.     push @B,(split /\s*\/\s*/,$_)  if $_ ne '';
    
32. }
    
33. close B;
    
34. 
    
35. for(1 .. 1000){
    
36.     my %rand = ();
    
37.     while(scalar keys %rand < $yangbengeshu){
    
38.         $rand{int(rand scalar @A)} = 0;
    
39.     }
    
40.     my @A_choose = map{ $A[$_] } keys %rand;
    
41.     my (%m,%n);
    
42.     for(@A_choose,@B){         
    
43.         $m{$_}++ && $n{$_}++;
    
44.     }      
    
45.     say scalar keys %n;
    
46. }

复制代码

结果是：
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
6
6
5
6
5
6
5
5
5
5
6
6
5
5
5
5
5
5
6
6
6
5
5
6
5
5
....

珞水的大叔

回复 27# bioinfor


    这是因为targetgene.txt中有重复样本
   程序改了一下，如下：

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. use 5.010;
   
5. 
   
6. 
   
7. my $yangbengeshu = 3;   #抽取样本个数
   
8. 
   
9. my (%A,@A);
   
10. open A,"<human_gene.txt";
    
11. while(<A>){
    
12.     chomp;
    
13.     s/^\s+$//;
    
14.     s/^\s+//;
    
15.     s/\s+$//;
    
16.     $A{$_} = undef if $_ ne '';
    
17. }
    
18. close A;
    
19. @A = keys %A;
    
20. if($yangbengeshu > scalar @A){
    
21.     die "ERROR: 抽取样本个数大于总个数\n";
    
22. }
    
23. 
    
24. my (%B,@B);
    
25. open B,"<targetgene.txt";
    
26. while(<B>){
    
27.     chomp;
    
28.     s/^\s+$//;
    
29.     s/^\s+//;
    
30.     s/\s+$//;
    
31.     map {$B{$_} = undef} (split /\s*\/\s*/,$_)  if $_ ne '';
    
32. }
    
33. close B;
    
34. @B = keys %B;
    
35. 
    
36. for(1 .. 1000){
    
37.     my %rand = ();
    
38.     while(scalar keys %rand < $yangbengeshu){
    
39.         $rand{int(rand scalar @A)} = 0;
    
40.     }
    
41.     my @A_choose = map{ $A[$_] } keys %rand;
    
42.     my (%m,%n);
    
43.     for(@A_choose,@B){         
    
44.         $m{$_}++ && $n{$_}++;
    
45.     }      
    
46.     say scalar keys %n;
    
47. }

复制代码

bioinfor

谢谢大叔！程序很强大！