RBP AS生信流程(N)RBP表达量提取

2021-01-24  本文已影响0人  cHarden13
use strict;
use warnings;
my %hash=();
open(RF,"RBPgene.txt") or die $!;
while(my $line=<RF>){
    chomp($line);
    $hash{$line}=1;
}
close(RF);
open(RF,"gene expression.txt") or die $!;
open(WF,">RBP expression.txt") or die $!;
while(my $line=<RF>){
    if($.==1){
        print WF $line;
        next;
    }
    my @arr=split(/\t/,$line);
    my @zeroArr=split(/\|/,$arr[0]);
    if(exists $hash{$zeroArr[0]}){
        print WF $line;
        delete($hash{$zeroArr[0]});
    }
}
close(WF);
close(RF);
上一篇 下一篇

猜你喜欢

热点阅读