Question

我想从@failedamplicons数组中检索包含类似元素（ like ）到@failedamplicons数组的元素。 @amplicon_exon中的每个元素都是唯一的，只能匹配@failedamplicons: example: OCP1_FGFR3_8.87 OCP1_AR_14.89 @amplicon_exon: example: TEST_Focus_ERBB2_2:22:ERBB2:GENE_ID=ERBB2;PURPOSE=CNV,Hotspot;CNV_ID=ERBB2;CNV_HS=1 OCP1_FGFR3_8:intron:FGFR3:GENE_ID=FGFR3;PURPOSE=CNV;CNV_ID=FGFR3;CNV_HS=1 OCP1_CDK6_14:intron:CDK6:GENE_ID=CDK6;PURPOSE=CNV;CNV_ID=CDK6;CNV_HS=1中的单个元素。我尝试了两个for循环，但获得了重复值。有没有更好的方法从两个数组中查找和检索类似的值？

my $i = 0;
my $j = 0;

for ( $i = 0; $i < @amplicon_exon; $i++ ) {

    for ( $j = 0; $j < @failedamplicons; $j++ ) {

        my $fail_amp = ( split /\./, $failedamplicons[$j] )[0];

        #print "the failed amp before match is $fail_amp\n";

        if ( index( $amplicon_exon[$i], $fail_amp ) != -1 ) {

            #print "the amplicon exon that matches $amplicon_exon[$i] and sample is $sample_id\n";
            print "the failed amp that matches $fail_amp and sample is $sample_id\n";

            my @parts = split /:/, $amplicon_exon[$i];
            my $exon_amp = $parts[1];

            next unless $parts[3] =~ /Hotspot/;    #includes only Hotspot amplicons
            my $gene_res   = $parts[2];
            my $depth      = ( split /\./, $failedamplicons[$j] )[1];
            my @total_amps = (
                $run_name, $sample_id, $gene_res, $depth, $fail_amp, $run_date, $matrix_status
            );

            my $lines = join "\t", @total_amps;

            push( @finallines, $lines );
        }
    }
}

这是两个for循环代码：

{{1}}

Answer 1

split和grep是您的朋友，迭代列表的惯用方法也是您的朋友。只需迭代第一个数组，只提取你要匹配的部分（使用split分割.字符上的元素，然后只取第一个条目），然后使用正则表达式， grep表示从元素开头到:的第二个数组中字符串的那一部分：

for my $elem (@failedamplicons){
    my $to_match = (split /\./, $elem)[0];
    if (my ($matched) = grep {$_ =~ /^\Q$to_match:/} @amplicon_exon){
        print "$matched\n";
    }        
}

Perl从两个数组中找到类似的元素

1 个答案: