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);
RBP AS生信流程(N)RBP表达量提取
最后编辑于 :
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
推荐阅读更多精彩内容
- 成长记录-连载(三十六) ——我的第一篇五千字长文,说了什么,你一定想不到 并不是不想每天写公众号,而是之前思考怎...
- 怎样反思自己的教学——学习杜威《我们怎样思维》有感 2015-01-13 14:56:52 来源:普宁市第三中学网...