多个文件求交集并生成矩阵

博主觉得写得还是有点复杂,求再简单一些的方法!

#!perl
use warnings;
use strict;

my %hash;
my %samp;

foreach my $f(@ARGV)
{
	my $id = $f;
	open FF, $f or die $!;
	while(<FF>)
	{
		chomp;
		my @tmp = split;
		$hash{$tmp[0]} = 0;
		$samp{$id}{$tmp[0]} = 0;
	}
	close FF;
}

my @name = sort {$a cmp $b} keys %samp;
print join "\t", "geneid", @name, "\n";

my %mat;
foreach my $k(keys %hash)
{
	foreach my $s(sort {$a cmp $b} keys %samp)
	{
		if(exists $samp{$s}{$k})
		{
			push @{$mat{$k}}, "yes";
		}else{
			push @{$mat{$k}}, "no";
		}
	}
}

foreach my $g(keys %mat)
{
	print join "\t", "$g", @{$mat{$g}}, "\n";
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值